Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,159,020 T→G 30.6% K237T (AAG→ACG)  proP_4 ← Proline/betaine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,159,0200TG30.6% 67.3 / 7.0 58K237T (AAG→ACG) proP_4Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base T (14/25);  new base G (17/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCA  >  ACB143/4158885‑4159154
                                                                                                                                        |                                                                                                                                      
gACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCt                                                                                                                                        <  6:39832/136‑1 (MQ=255)
 aCAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCtt                                                                                                                                       >  4:256894/1‑136 (MQ=255)
     gtaggtCGGCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACTCCCATGCCGATTGCCAGACTGCGCCGGTGCTTGCCCAGCAGTTCTATAAGCTGCGTCGTGGCCATGCCGGCCTTGGtt                                                                                                                                   <  4:272825/136‑1 (MQ=255)
        ggtCGGCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCTTCGTGGCCCTGCCCGCggtggtgtac                                                                                                                                >  7:247647/1‑132 (MQ=255)
        ggtCGGCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATTCCGGCCTTGGTTTcc                                                                                                                                >  6:287863/1‑136 (MQ=255)
        ggtCGGCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCCGCCGTGGttta                                                                                                                                 >  1:277921/1‑134 (MQ=255)
             gCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATTGCGGCCGTGGTTTCCAAgcg                                                                                                                           >  2:225921/1‑136 (MQ=255)
               aTGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGGCCATCCGGCCGTGGGTTACAGgcggg                                                                                                                         >  3:153998/1‑134 (MQ=255)
               aTGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCgct                                                                                                                         <  8:195373/136‑1 (MQ=255)
                  tAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGTTGACACCCAAGCCGATTTCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGAAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCgctgc                                                                                                                       <  5:75801/135‑1 (MQ=255)
                  tAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGGCATGCCGGCCGTTGTTTGCACGCgctgcg                                                                                                                      >  4:77875/1‑135 (MQ=255)
                  tAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCTGCCGGCCGTGGTTTCCAGGCgctgct                                                                                                                      >  8:80529/1‑136 (MQ=255)
                   aaCTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTc                                                                                                                     <  3:36094/136‑1 (MQ=255)
                    aCTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTc                                                                                                                     <  2:229617/135‑1 (MQ=255)
                        aGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCAtc                                                                                                                 <  1:225921/135‑1 (MQ=255)
                           aaCAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGGCAAGGCGGCCTTTGTTTCCAGGGGCTGCTCCATCTCCa                                                                                                             >  6:164182/1‑136 (MQ=255)
                           aaCAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCTGCCGGCCGTGGTTTCCCGGCGCTGCTCCAtctc                                                                                                               >  3:288385/1‑134 (MQ=255)
                           aaCAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTTGTTTCCAGGCGCTGCTCCAtatc                                                                                                               >  2:211579/1‑134 (MQ=255)
                           aaCAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCAtctc                                                                                                               >  1:180210/1‑134 (MQ=255)
                           aaCAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCa                                                                                                             >  1:157847/1‑136 (MQ=255)
                             cAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCAGCCGGCCTTGGGTTGCACGCGCGGCTCCAACTCCCgg                                                                                                           >  3:102241/1‑136 (MQ=255)
                                     aCGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCGTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTggg                                                                                                   >  3:291934/1‑136 (MQ=255)
                                           ccGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCTGCCGGCCGTGGTTTCCCGGCGCTGCGGCCTCTCCCGGAACTTGGGCggggc                                                                                             >  1:8248/1‑132 (MQ=255)
                                            cGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGGCGTGGTTTCCAGGCGCGGCGCCATCGCCCGGAACTTGGGCGTGGcc                                                                                            >  1:225613/1‑136 (MQ=255)
                                              cATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCTGCCGGCCTTGGTTTCCCGGCGCTGCTCCATCTCCAGGAACTTTGGCGTGTCCTg                                                                                          >  2:85391/1‑136 (MQ=255)
                                               aTTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGc                                                                                         <  3:285983/136‑1 (MQ=255)
                                                   aggcaggTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGc                                                                                     <  8:170226/136‑1 (MQ=255)
                                                   aggcaggTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGc                                                                                     <  8:146713/136‑1 (MQ=255)
                                                      caggTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCa                                                                                  <  1:145552/136‑1 (MQ=255)
                                                           gACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCGTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCgg                                                                             >  7:9637/1‑136 (MQ=255)
                                                            acacCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCAGCCGGCCGTGGTTTCCCGGCGCGGGGCCCTCTGCCCGAACGTGGGCGGGGCCGGGAAGCTCCTGCGGa                                                                            >  3:262357/1‑136 (MQ=255)
                                                            acacCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGa                                                                            <  7:57343/136‑1 (MQ=255)
                                                                ccATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGAtata                                                                        <  1:146652/136‑1 (MQ=255)
                                                                      cGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCAGCCGGCCTTGGTTTCCAGGCGCTGCAACATCTCCAGGAACTTGGGCGTGTCCAGCAAACTCAAGCCGATATAACGcc                                                                   >  2:225442/1‑135 (MQ=255)
                                                                      cGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCCTCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATATCGcc                                                                   >  2:80863/1‑135 (MQ=255)
                                                                      cGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCGGCTCCATATCCAAGAACTTTGGCGTGTCCATCCAGCTCATGCGGCTTTAGCGccc                                                                  >  8:64493/1‑136 (MQ=255)
                                                                        aTTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGGGTCCTGCAAGCTCCTGCGGATATAGCGCCCTa                                                                >  4:130303/1‑136 (MQ=255)
                                                                           gCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCTGCCGGCCGTGGTTTCCAGGCGCTGGTCCCTCTCCAGGAACTTGGGCGGGTCCTGCAAGCTCAAGCGGATATAGCGCCCTACCa                                                             >  2:45799/1‑136 (MQ=255)
                                                                                    cgcCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACg                                                    <  3:76109/136‑1 (MQ=255)
                                                                                     gcCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCGTGGTTTCCAGGCGCTGCTCCCTCTCCAGGAACTTGGGCGTGGGCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGGACgg                                                   >  1:206517/1‑136 (MQ=255)
                                                                                       cGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCCTGCCGGGCGTGGTTTCCCGGCGCGGCGGCCTGTCCAGGAACTTGGGCGTGTCCTGCCAGCGCATGCGGATATAGCGCCCGACCAGGCCGCACGgcg                                                 >  2:164771/1‑136 (MQ=255)
                                                                                            cTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCgccgcc                                            <  5:171075/136‑1 (MQ=255)
                                                                                               gCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCGTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCcagc                                         <  8:178523/136‑1 (MQ=255)
                                                                                                 ccAGCAGTTCGCGCAGCGGCGTCGTGGCCCTGCCGGCCGTGGTTTGCCGGCGCGGCGCCCTCTCCAGGCACTTGGGCGTGGGCGGCCAGCTCATGGGGGTATAGCGCCCGACCCGGCCGCACCGCGCCGCccgcag                                       >  2:258861/1‑136 (MQ=255)
                                                                                                       gTTCGCGCAGCGGCGTCGTGGGCAAGCCGGCCGTGGTTTCCAGGCGCTGCGCCCTATACATGTAAATGGGCGTGTGCTGGATTGATGTGGGGTTATTGCTGCCGACCACGCCGGACGGCGCCGCCAGCGGGGACgg                                 >  8:35118/1‑136 (MQ=255)
                                                                                                           gcgcAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAgg                             <  8:79170/136‑1 (MQ=255)
                                                                                                            cgcAGCGGCGTCGTGGCCCTGCCGGCCGTGGTTTCCAGGCGCGGCGCCCTCTCCAGGGACTTTGGGGGGTGCGGCCAGCTCATGCGGATTTTGCGCCCGACCCGGCCCCCCGGCGCCGCCCGCCGGCAcggccggc                            >  4:54332/1‑136 (MQ=255)
                                                                                                             gcAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTGCTTCTCCAGGAACTTGGGCGTGTCCTGCCAGCTCATGCGGATATAGCGCCCTTCCAGGCCGAACGGCGCCGCCCGCAGGAACGGCAGgcg                           >  3:148638/1‑136 (MQ=255)
                                                                                                               aGCGGCGGCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGcc                         <  5:112389/136‑1 (MQ=255)
                                                                                                                 cggcgTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGTAGGCGCCAg                       <  3:213860/136‑1 (MQ=255)
                                                                                                                    cgtcgtGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGccc                    <  7:52902/136‑1 (MQ=255)
                                                                                                                         tGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTc               <  2:277921/136‑1 (MQ=255)
                                                                                                                          gggCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTc               <  4:260071/135‑1 (MQ=255)
                                                                                                                             ccATGCCGGCCGTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATg            <  6:267321/136‑1 (MQ=255)
                                                                                                                                  ccGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGTTCTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAggg       <  2:176252/136‑1 (MQ=255)
                                                                                                                                  ccGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAggg       <  5:86392/136‑1 (MQ=255)
                                                                                                                                    ggCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAgc                                                                      <  1:182055/71‑1 (MQ=255)
                                                                                                                                    ggCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAgc                                                                      >  2:182055/1‑71 (MQ=255)
                                                                                                                                    ggCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATACCGCCCAACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCACGCGCCAGCCCCACTCATGCAGGGCt     >  7:24817/1‑136 (MQ=255)
                                                                                                                                       cTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCa  <  3:85375/136‑1 (MQ=255)
                                                                                                                                       cTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCa  <  2:157847/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GACAGGTAGGTCGGCATGTAACTGAGCAACAGGTAGAACGCAACCGCATTCAGGCAGGTGACACCCATGCCGATTGCCAGACTGCGCCGGTGCTGGCCCAGCAGTTCGCGCAGCGGCGTCGTGGCCATGCCGGCCTTGGTTTCCAGGCGCTGCTCCATCTCCAGGAACTTGGGCGTGTCCTGCAAGCTCATGCGGATATAGCGCCCTACCAGGCCGAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCA  >  ACB143/4158885‑4159154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: