Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,161,134 A→T 5.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,161,1340AT 5.0% 1335.6 / 36.0 388L217* (TTA→TAA) PP_3659transcriptional regulator, AraC family
Reads supporting (aligned to +/- strand):  ref base A (174/193);  new base T (11/9);  total (186/202)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 414 total aligned reads displayed.

GGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACG  >  NC_002947/4160986‑4161282
                                                                                                                                                    |                                                                                                                                                    
 gTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaa                                                                                                                                                    >  1:4562156/1‑150 (MQ=255)
  tGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaa                                                                                                                                                   <  2:5883310/150‑1 (MQ=255)
   ggAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaaa                                                                                                                                                  <  2:4706/150‑1 (MQ=255)
    gAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAg                                                                                                                                                 >  2:6543074/1‑150 (MQ=255)
     aaCTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAgg                                                                                                                                                >  2:518424/1‑150 (MQ=255)
       cTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCa                                                                                                                                              <  2:3872419/150‑1 (MQ=255)
       cTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCa                                                                                                                                              <  1:5790188/150‑1 (MQ=255)
       cTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCa                                                                                                                                              <  1:2186680/150‑1 (MQ=255)
        tGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAg                                                                                                                                             <  1:4753601/150‑1 (MQ=255)
           aCATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGaa                                                                                                                                          <  2:8570656/150‑1 (MQ=255)
           aCATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGaa                                                                                                                                          >  2:7249018/1‑150 (MQ=255)
            cATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAag                                                                                                                                         <  1:7039902/150‑1 (MQ=255)
             aTGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAaga                                                                                                                                        >  2:2341633/1‑150 (MQ=255)
                 tCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGaaa                                                                                                                                    >  1:2846987/1‑150 (MQ=255)
                    aGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTCTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACgg                                                                                                                                 <  2:6030742/150‑1 (MQ=255)
                      tCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGc                                                                                                                               <  1:2043001/150‑1 (MQ=255)
                      tCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGc                                                                                                                               <  2:5985561/150‑1 (MQ=255)
                         cgGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGtt                                                                                                                            <  2:1518049/150‑1 (MQ=255)
                          gggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                           >  1:6384660/1‑150 (MQ=255)
                          gggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                           <  1:5558743/150‑1 (MQ=255)
                           ggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTg                                                                                                                          <  2:2085302/150‑1 (MQ=255)
                           ggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTg                                                                                                                          <  2:3591127/150‑1 (MQ=255)
                            gCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTgg                                                                                                                         <  1:4394023/150‑1 (MQ=255)
                            gCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTgg                                                                                                                         <  2:2534033/150‑1 (MQ=255)
                            gCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTgg                                                                                                                         <  2:5620008/150‑1 (MQ=255)
                                gCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCaaa                                                                                                                     >  2:8304912/1‑150 (MQ=255)
                                    aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCATGAAGAGAAACGGGCGTTTGGCAAACCtt                                                                                                                 >  1:2984585/1‑150 (MQ=255)
                                    aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCtt                                                                                                                 >  2:2194522/1‑150 (MQ=255)
                                    aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCtt                                                                                                                 <  1:6705111/150‑1 (MQ=255)
                                    aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCtt                                                                                                                 <  1:6575163/150‑1 (MQ=255)
                                     tCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTg                                                                                                                >  1:8566749/1‑150 (MQ=255)
                                       gAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTa                                                                                                              <  2:6729685/150‑1 (MQ=255)
                                          ttCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGc                                                                                                           >  1:8261140/1‑150 (MQ=255)
                                           tCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGcc                                                                                                          <  1:8802043/150‑1 (MQ=255)
                                            cAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCa                                                                                                         <  1:4448173/150‑1 (MQ=255)
                                             aGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAg                                                                                                        >  1:1179396/1‑150 (MQ=255)
                                             aGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAg                                                                                                        >  2:7013921/1‑150 (MQ=255)
                                               gCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGaa                                                                                                      >  1:6060372/1‑150 (MQ=255)
                                                ccACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGACGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCATTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGaaa                                                                                                     <  2:436834/150‑1 (MQ=255)
                                                  aCGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGc                                                                                                   >  1:7804114/1‑150 (MQ=255)
                                                           gagagaGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGa                                                                                          <  2:5421941/150‑1 (MQ=255)
                                                            agagagCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGaa                                                                                         <  1:6404182/150‑1 (MQ=255)
                                                            agagagCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGaa                                                                                         <  1:8328415/150‑1 (MQ=255)
                                                              agagCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAaca                                                                                       <  1:516570/150‑1 (MQ=255)
                                                                 gCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAg                                                                                    >  2:1327473/1‑150 (MQ=255)
                                                                  cAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGc                                                                                   >  2:8431219/1‑150 (MQ=255)
                                                                    gTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCgg                                                                                 >  1:4904701/1‑150 (MQ=255)
                                                                     ttttAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGa                                                                                >  2:2691666/1‑150 (MQ=255)
                                                                     ttttAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGa                                                                                <  1:337327/150‑1 (MQ=255)
                                                                      tttAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAg                                                                               <  1:1389292/150‑1 (MQ=255)
                                                                         aaTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGtt                                                                            >  1:7378171/1‑150 (MQ=255)
                                                                         aaTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGtt                                                                            >  1:5996441/1‑150 (MQ=255)
                                                                           tCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCt                                                                          >  1:3887884/1‑150 (MQ=255)
                                                                             gccgccGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTAt                                                                        <  2:1643341/150‑1 (MQ=255)
                                                                             gccgccGGATCATAATGCTCGTACTGTCGTCGTTACGGTATATCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTAt                                                                        >  2:4517690/1‑150 (MQ=255)
                                                                                   ggATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGccc                                                                  >  1:8679390/1‑150 (MQ=255)
                                                                                   ggATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGccc                                                                  >  1:7525215/1‑150 (MQ=255)
                                                                                     aTCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTc                                                                <  2:8797368/150‑1 (MQ=255)
                                                                                      tCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTcc                                                               >  2:2969922/1‑150 (MQ=255)
                                                                                      tCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCa                                                              >  2:4637802/1‑150 (MQ=255)
                                                                                       cATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCa                                                              <  1:7885262/150‑1 (MQ=255)
                                                                                       cATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCa                                                              >  2:3370466/1‑150 (MQ=255)
                                                                                        aTAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCAc                                                             >  1:4195560/1‑150 (MQ=255)
                                                                                               tCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGAttt                                                      >  2:7781924/1‑150 (MQ=255)
                                                                                                 gTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTcc                                                    <  1:6839213/150‑1 (MQ=255)
                                                                                                       tcgtcgTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACaa                                              >  2:3282483/1‑150 (MQ=255)
                                                                                                          tcgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAg                                           <  2:4148219/150‑1 (MQ=255)
                                                                                                           cgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAgc                                          <  1:1012055/150‑1 (MQ=255)
                                                                                                           cgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAgc                                          <  2:7473895/150‑1 (MQ=255)
                                                                                                              tACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCa                                       <  1:5670450/150‑1 (MQ=255)
                                                                                                              tACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCa                                       >  1:8525375/1‑150 (MQ=255)
                                                                                                                 ggTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTc                                    >  2:4955917/1‑150 (MQ=255)
                                                                                                                  gTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCaaa                                                                                                                     <  2:1926267/68‑1 (MQ=255)
                                                                                                                  gTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCaaa                                                                                                                     >  1:1926267/1‑68 (MQ=255)
                                                                                                                   tATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGa                                  <  2:4894651/150‑1 (MQ=255)
                                                                                                                   tATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGa                                  >  1:5375151/1‑150 (MQ=255)
                                                                                                                     tttCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATa                                <  1:8110633/150‑1 (MQ=255)
                                                                                                                        cacaAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTc                             >  1:2692139/1‑150 (MQ=255)
                                                                                                                         acaAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCg                            >  2:810963/1‑150 (MQ=255)
                                                                                                                          caAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGc                           <  1:1219629/150‑1 (MQ=255)
                                                                                                                           aaGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCt                          >  1:3590139/1‑150 (MQ=255)
                                                                                                                           aaGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCt                          <  2:5797998/150‑1 (MQ=255)
                                                                                                                              gAAGTTACTTGTCTTGCATATTAAAAAGTCAGGAAGAGAAACGGGCGTTTGTCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTgtc                       <  2:7869239/150‑1 (MQ=255)
                                                                                                                              gAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTgtc                       <  1:2497439/150‑1 (MQ=255)
                                                                                                                                      ttGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCg               >  2:6366842/1‑150 (MQ=255)
                                                                                                                                      tcttcctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCg               >  1:3179460/7‑150 (MQ=255)
                                                                                                                                        ggCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGc             >  1:6252129/1‑150 (MQ=255)
                                                                                                                                           gTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga          <  2:6686858/150‑1 (MQ=255)
                                                                                                                                           gTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga          <  1:2628904/150‑1 (MQ=255)
                                                                                                                                           ctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga          >  1:734328/2‑150 (MQ=255)
                                                                                                                                           ctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga          >  2:1224358/2‑150 (MQ=255)
                                                                                                                                            tGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgag         >  2:4194167/1‑150 (MQ=255)
                                                                                                                                            tGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgag         <  1:5911365/150‑1 (MQ=255)
                                                                                                                                               cTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgaggag      <  2:5346383/150‑1 (MQ=255)
                                                                                                                                               cTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgaggag      <  2:1836817/150‑1 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGaggagg     >  2:3064131/1‑150 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGaggagg     >  1:3951458/1‑150 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGaggagg     >  1:1686381/1‑150 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGaggagg     >  2:4284430/1‑150 (MQ=255)
                                                                                                                                                  ttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGAc   <  1:462429/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACG  >  NC_002947/4160986‑4161282

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