Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,100,562:1 +C 7.0% coding (1341/1470 nt) PP_2709 → long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,100,5621.C 7.0% 220.2 / 9.8 57coding (1341/1470 nt)PP_2709long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (28/25);  new base C (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

AAGTGGTGTGATCGCCCAGGCTTTCGTGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGTTGTGAGGTTCACAATGTCCTT  >  NC_002947/3100415‑3100707
                                                                                                                                                    |                                                                                                                                                  
aaGTGGTGTGATCGCCCAGGCTTTCGTGCTTGTCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGTCATCAGAGACGAGGCCATCGATCAGGGCGTACAGCAATGCAACGCGCGATTGCGCGATTACGCCAggg                                                                                                                                                   <  3:35109/149‑1 (MQ=255)
       gtgAGCGCCCCTGGTTTCGGTGATGGCGTAGCCTTGCCCCGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGc                                                                                                                                           <  8:52766/149‑1 (MQ=255)
          aTCGCCCAGGCTTTCGTGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGCGTGCAccctg                                                                                                                                         >  7:322479/1‑145 (MQ=255)
              cccAGGCTTTCTTGCATGGCGAAGCCTTGCCACGCAACCATGCCTTGATCTGGCCGCTGGATCCCGGCATCAGCGACGAGCCCATCGAACAGGCCGTACAGCAATGCAACGCGCGAATGCCCGATTACGCCAGGGTGCACGCCTGgcgc                                                                                                                                    >  1:137775/1‑149 (MQ=255)
              cccAGGCTTTCGTGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGgcgc                                                                                                                                    >  8:86175/1‑149 (MQ=255)
              cccAGGCTTTCGTGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGgcgc                                                                                                                                    >  7:75185/1‑149 (MQ=255)
                 aGGCTTTCGTGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGcgccgc                                                                                                                                 >  8:183880/1‑149 (MQ=255)
                   gCTTTCGTGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGCGTGCACCCTGgcgtaatcg                                                                                                                               <  8:322479/149‑7 (MQ=255)
                          tGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGCGTGCACCCTGgcg                                                                                                                                     <  2:138728/136‑1 (MQ=255)
                          tGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGCGTGCACCCTGgcg                                                                                                                                     >  1:138728/1‑136 (MQ=255)
                            cATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGc                                                                                                                       <  3:180943/148‑1 (MQ=255)
                                 cGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCctct                                                                                                                   >  1:93996/1‑147 (MQ=255)
                                 cGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTa                                                                                                                 <  3:185745/149‑1 (MQ=255)
                                       cTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGc                                                                                                           <  8:75185/149‑1 (MQ=255)
                                        ttGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCt                                                                                                          >  5:243213/1‑149 (MQ=255)
                                          gCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGTCAGTGCCTCTTAGCACGCtt                                                                                                         <  4:192021/148‑1 (MQ=255)
                                             aCGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGc                                                                                                                       >  6:106835/1‑131 (MQ=255)
                                             aCGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGc                                                                                                                       <  5:106835/131‑1 (MQ=255)
                                                 aaCCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAAc                                                                                                 >  7:217527/1‑149 (MQ=255)
                                                   ccttgccttgCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAAccc                                                                                               >  3:112138/1‑149 (MQ=255)
                                                     ttgccttgcTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGaaa                                                                                                  <  3:303808/144‑1 (MQ=255)
                                                     ttgccttgcTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGaaa                                                                                                  >  4:303808/1‑144 (MQ=255)
                                                        ccttgcTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGAc                                                                                           <  4:20909/148‑1 (MQ=255)
                                                           tgcTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACTCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGc                                                                                        >  3:319683/1‑148 (MQ=255)
                                                            gcTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGc                                                                                        <  7:183880/147‑1 (MQ=255)
                                                             ctctGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCaa                                                                                     >  7:136965/1‑149 (MQ=255)
                                                                  gCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGgc                                                                                 <  3:135091/148‑1 (MQ=255)
                                                                    cGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGgccgc                                                                              >  7:159302/1‑149 (MQ=255)
                                                                     gCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGgccgcc                                                                             >  8:216451/1‑149 (MQ=255)
                                                                             ccGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGc                                                                                        >  5:265016/1‑130 (MQ=255)
                                                                             ccGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGc                                                                                        <  6:265016/130‑1 (MQ=255)
                                                                             ccGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGt                                                                     >  2:224949/1‑149 (MQ=255)
                                                                              cGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGt                                                                     >  7:272504/1‑148 (MQ=255)
                                                                                gCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAg                                                                  <  7:216451/149‑1 (MQ=255)
                                                                                    cAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGa                                                              >  5:3738/1‑149 (MQ=255)
                                                                                             ggCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGc                                                     <  6:243213/149‑1 (MQ=255)
                                                                                              gCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACTCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCg                                                    <  4:319683/149‑1 (MQ=255)
                                                                                               ccATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGt                                                   <  1:224949/149‑1 (MQ=255)
                                                                                                cATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGt                                                   >  8:105133/1‑148 (MQ=255)
                                                                                                cATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGt                                                   >  7:268519/1‑148 (MQ=255)
                                                                                                cATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATAGCCGCCCGCGACGTGAGGCGATCTTGCAGCGt                                                   >  2:31291/1‑148 (MQ=255)
                                                                                                 aTCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTa                                                 >  3:140634/1‑149 (MQ=255)
                                                                                                 aTCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTa                                                 >  7:85901/1‑149 (MQ=255)
                                                                                                      aCAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACcaca                                            <  4:112138/149‑1 (MQ=255)
                                                                                                       cAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCacac                                           >  6:280606/1‑149 (MQ=255)
                                                                                                           ccGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCtg                                       <  8:268519/149‑1 (MQ=255)
                                                                                                                   cAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGa                               <  6:3738/149‑1 (MQ=255)
                                                                                                                     aTGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACa                             <  4:140634/149‑1 (MQ=255)
                                                                                                                     aTGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACa                             <  8:217527/149‑1 (MQ=255)
                                                                                                                       gCAACGCGCGATTGCCCGATTACGCCAGGGTGCAGCGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGGGGCGATCTTg                                                         >  1:232086/1‑120 (MQ=255)
                                                                                                                       gCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGGGGCGATCTTg                                                         <  2:232086/119‑1 (MQ=255)
                                                                                                                            gcgcgATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGtt                      <  5:280606/149‑1 (MQ=255)
                                                                                                                            gcgcgATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTCCCTCTTAGCTCGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGtt                      <  2:93996/149‑1 (MQ=255)
                                                                                                                                 aTTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGTTGTGAg                 <  7:86175/149‑1 (MQ=255)
                                                                                                                                              gCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGTTGTGAGGTTCACAATGTcc    >  7:34460/1‑149 (MQ=255)
                                                                                                                                               ccAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGTTGTGAGGTTCACAATGTc     <  1:117591/147‑1 (MQ=255)
                                                                                                                                               ccAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGTTGTGAGGTTCACAATGTc     >  2:117591/1‑147 (MQ=255)
                                                                                                                                                cAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGTTGTGAGGTTCACAATGTCCtt  >  3:179260/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGTGGTGTGATCGCCCAGGCTTTCGTGCATGGCGAAGCCTTGCCACGCAACCTTGCCTTGCTCTGGCCGCTGGATCCCGGCATCAGCGACGAGGCCATCGAACAGGCCGTACAGCAATGCAACGCGCGATTGCCCGATTACGCCAGGGTGCACGCCTGGCGCCGCTTGCCAGTGCCTCTTAGCACGCTTGACGAAACCCTGACTGCCAATGGCCGCCCGCGACGTGAGGCGATCTTGCAGCGTTACCACACCCTGTTGTCCGACATCTCGTTGTGAGGTTCACAATGTCCTT  >  NC_002947/3100415‑3100707

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