Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,156,141 A→T 11.0% L217* (TTA→TAA)  PP_3659 ← transcriptional regulator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,156,1410AT11.0% 189.9 / 13.5 80L217* (TTA→TAA) PP_3659transcriptional regulator, AraC family
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base T (6/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

GGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGA  >  NC_002947/4155996‑4156284
                                                                                                                                                 |                                                                                                                                               
ggAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTtaa                                                                                                                                               >  1:241896/1‑148 (MQ=255)
 gAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaaa                                                                                                                                             <  4:173623/149‑1 (MQ=255)
       gACATGCTCAAGTCTCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTTAAAAGGCac                                                                                                                                        <  8:41416/148‑2 (MQ=255)
       gACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCa                                                                                                                                         <  8:241686/147‑1 (MQ=255)
       gACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAg                                                                                                                                        >  7:256702/1‑148 (MQ=255)
       gACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAg                                                                                                                                        >  5:241109/1‑148 (MQ=255)
         cATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACTTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGCCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGaa                                                                                                                                     <  3:8835/149‑1 (MQ=255)
         cATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGaa                                                                                                                                     <  4:20201/149‑1 (MQ=255)
           tGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTTAAAAGGCacgccaa                                                                                                                                   <  2:241896/149‑7 (MQ=255)
           tGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCCTGCCTTTTTAAAAGGCacgccaa                                                                                                                                   <  3:337568/149‑7 (MQ=255)
              tCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAgaga                                                                                                                                 >  4:7624/1‑148 (MQ=255)
                aaGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGaaa                                                                                                                               <  3:7624/148‑1 (MQ=255)
                  gTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGCATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACgg                                                                                                                            >  6:109724/1‑149 (MQ=255)
                     gcgGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGTAGAGAGCAGTTTTAATCGCAGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTGTTAAAAAGGCAGGAAGAGAAACGGGCg                                                                                                                         >  8:223760/1‑149 (MQ=255)
                       gggCGAGCTCATCGAGTTCAGGCCACGGGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGtt                                                                                                                       <  3:7579/149‑1 (MQ=255)
                        ggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                      <  2:209817/149‑1 (MQ=255)
                         gCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTg                                                                                                                     >  5:175995/1‑149 (MQ=255)
                             gCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCa                                                                                                                  >  5:112730/1‑148 (MQ=255)
                               tCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCTTTTGGCaaa                                                                                                                >  4:49103/1‑148 (MQ=255)
                               tCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCaaa                                                                                                                >  3:18124/1‑148 (MQ=255)
                                 aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCt                                                                                                             >  3:253907/1‑149 (MQ=255)
                                 aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCt                                                                                                             >  1:277022/1‑149 (MQ=255)
                                  tCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTAATTGTCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGTAGTTTGGCAAAACtt                                                                                                            >  5:191940/1‑149 (MQ=255)
                                          aGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGc                                                                                                      >  2:291137/1‑147 (MQ=255)
                                          aGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCa                                                                                                    <  4:253907/149‑1 (MQ=255)
                                           ggCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAg                                                                                                   <  3:327063/149‑1 (MQ=255)
                                            gCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAg                                                                                                   >  4:34703/1‑148 (MQ=255)
                                                       ggAGAGAGAAGTTTTAATCGCCGCAGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGCGCGTTTGGCAAAACTTGTACGCCAGAAACCtcttc                                                                                         >  7:140281/1‑147 (MQ=255)
                                                        gagagaGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTtctc                                                                                       >  1:159362/1‑148 (MQ=255)
                                                            gagCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACGTTGTACGCCAGAAAGCTCTTCTCGAaca                                                                                  <  5:374066/149‑1 (MQ=255)
                                                               cAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAg                                                                               <  8:106576/149‑1 (MQ=255)
                                                                 gTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCg                                                                             <  2:159362/149‑1 (MQ=255)
                                                                  ttttAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAgag                                                                             >  7:278573/1‑146 (MQ=255)
                                                                  ttttAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCg                                                                             >  6:286934/1‑148 (MQ=255)
                                                                  ttttAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCAAAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACAAAgag                                                                             >  6:127064/1‑146 (MQ=255)
                                                                       aTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCGGAAAGCTCTTCTCGAACACAGCGGAGGtt                                                                       >  3:287230/1‑149 (MQ=255)
                                                                       aTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGt                                                                        >  8:29960/1‑148 (MQ=255)
                                                                           ccgccgGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTAt                                                                   <  8:217454/149‑1 (MQ=255)
                                                                              ccgGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCgg                                                                >  7:247148/1‑149 (MQ=255)
                                                                              ccgCATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCgg                                                                <  5:109724/149‑1 (MQ=255)
                                                                                  aTCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGcc                                                              <  6:241109/147‑1 (MQ=255)
                                                                                   tCATAATGCTGGTACTGTCGTCGTTACGGTATTTCACAAGTAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGTAAGCTCTTCTCGAACACAGCGGAGGTTCGATCGGCCCTc                                                           <  2:5805/149‑1 (MQ=255)
                                                                                       aaTGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACt                                                       >  1:198849/1‑149 (MQ=255)
                                                                                       aaTGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACt                                                       >  5:85614/1‑149 (MQ=255)
                                                                                           cTCGTACTGTCGTCGTTTCGGTATTTCACAAGGATGTTACTTGGCGTGCGTTTTAAAAAGGCAGGATGAGTAACGGGCTTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACAGAGCGGAGGTTCTATCGGCCCTCCACTCGAt                                                   <  8:140281/149‑1 (MQ=255)
                                                                                           cTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGa                                                    <  2:198849/148‑1 (MQ=255)
                                                                                            tCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGAtt                                                  >  3:30608/1‑149 (MQ=255)
                                                                                            tCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGAtt                                                  >  7:8457/1‑149 (MQ=255)
                                                                                            tCGTACTGTCGTCGTTACGCTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGAtt                                                  >  2:6547/1‑149 (MQ=255)
                                                                                               tACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCGGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTcc                                               <  4:287230/149‑1 (MQ=255)
                                                                                                    tcgtcgTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTAc                                                                                                        <  6:137847/87‑1 (MQ=255)
                                                                                                    tcgtcgTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTAc                                                                                                        >  5:137847/1‑87 (MQ=255)
                                                                                                    tcgtcgTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTcac                                           >  1:235868/1‑148 (MQ=255)
                                                                                                    tcgtcgTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTca                                            <  6:85614/147‑1 (MQ=255)
                                                                                                       tcgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCtt                                                                                                            >  6:150201/1‑80 (MQ=255)
                                                                                                          ttACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAgcg                                    >  8:214418/1‑149 (MQ=255)
                                                                                                            aCGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCg                                                                 <  4:350704/118‑1 (MQ=255)
                                                                                                            aCGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCa                                  >  8:261146/1‑149 (MQ=255)
                                                                                                             cGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCa                                  <  3:49103/148‑1 (MQ=255)
                                                                                                              ggTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCAtt                                >  4:258291/1‑149 (MQ=255)
                                                                                                               gTATTTCACAATGTAGTTACTTGGGGTGCCTTTTAAAAAGGCAGTAAGAGAATCGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAAAACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTc                               <  8:256702/149‑1 (MQ=255)
                                                                                                                    tCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGATGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATa                           <  7:223760/148‑1 (MQ=255)
                                                                                                                      acaAGGAAGTTACTTGGCGTGCCTTTTAAAAAGTCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTc                        <  6:112730/149‑1 (MQ=255)
                                                                                                                        aaGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGc                      >  7:46138/1‑149 (MQ=255)
                                                                                                                         aGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCt                     <  4:18124/149‑1 (MQ=255)
                                                                                                                              gTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTgtcgt                >  5:45954/1‑149 (MQ=255)
                                                                                                                                   ttGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCAcc           >  7:54103/1‑149 (MQ=255)
                                                                                                                                     ttcctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCGTCCACTCGATTTCCGCACAAAAGCGCATTCGATAGTCGCTGTCGTCCAGCgg         >  2:164114/5‑149 (MQ=255)
                                                                                                                                     ttcctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCgg         >  4:77185/5‑149 (MQ=255)
                                                                                                                                      gCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTc                                             >  7:274581/1‑112 (MQ=255)
                                                                                                                                      gCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTc                                             <  8:274581/112‑1 (MQ=255)
                                                                                                                                      gCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCgg         <  2:277022/148‑1 (MQ=255)
                                                                                                                                       cctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTACACTCGATTTCCTCCCAAAAGCGCATTCGATAGTCTCTGTCGTCCACCGGCt       >  8:155762/3‑149 (MQ=255)
                                                                                                                                       cGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCt       <  7:261146/149‑1 (MQ=255)
                                                                                                                                       agTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACATTGTACGCAAGAAAGCTCTTCTCGAACACAGAGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCAAAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCt       <  8:278573/148‑1 (MQ=255)
                                                                                                                                        ctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACGCAGCGGAGGTTCTATCGGCCCTCGACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTg      >  8:240101/2‑149 (MQ=255)
                                                                                                                                        ctGCCTTTTTAAAAGGCAGGAAGAAAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCt       >  8:134703/2‑148 (MQ=255)
                                                                                                                                         ctccTTTTAATAACGCTGGAAGAGATACGTGGGTTTGGCAAACGTTGTACGCCAGATATTTCTTATCTATCACTGCGGATTTTCTATCGGCCCTCCACTCGATTTCGTGTCAATAGCCTATTGGATAGTCGCTGTGGTCGACCTGCTga     <  4:311511/147‑1 (MQ=255)
                                                                                                                                          gCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgag    >  7:99422/1‑149 (MQ=255)
                                                                                                                                            cTTTTAAAAAGTCATGAAGAGAAACGGGCGTTTGGCAAAACTTGTACGCCAGATAGCTCTTCTCGAACACATCGGAGGTTCTATCTGCCCTACACTCGATTTCCTGACAAAAGCTCATTGGATAGTCGCTGTCGTCCACCGGCTgagg   <  6:191940/148‑1 (MQ=255)
                                                                                                                                            cTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgagga  <  7:263286/149‑1 (MQ=255)
                                                                                                                                            cTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgagga  <  1:291137/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGA  >  NC_002947/4155996‑4156284

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: