Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,247,728 G→A P269S (CCG→TCG)  lysP ← lysine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,247,7280GA99.0% 342.0 / ‑4.8 98P269S (CCG→TCG) lysPlysine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (40/57);  total (40/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GGAGTTACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGGATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGAT  >  NC_000913/2247591‑2247863
                                                                                                                                         |                                                                                                                                       
ggAGTTACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAAt                                                                                                                                       <  2:1103059/140‑1 (MQ=255)
 gAGTTCCCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATc                                                                                                                                      <  2:631174/140‑1 (MQ=255)
 gAGTTACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATc                                                                                                                                      >  1:1222527/1‑140 (MQ=255)
 gAGTTACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATc                                                                                                                                      <  2:1201955/140‑1 (MQ=255)
 gAGTTACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATc                                                                                                                                      >  1:199694/1‑140 (MQ=255)
  aGTTACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCg                                                                                                                                     >  1:464500/1‑140 (MQ=255)
    ttACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgt                                                                                                                                   <  2:1085113/140‑1 (MQ=255)
    ttACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgt                                                                                                                                   <  1:460658/140‑1 (MQ=255)
     tACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtg                                                                                                                                  >  1:962675/1‑140 (MQ=255)
      ccccGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 <  1:143464/139‑1 (MQ=255)
      aCCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 >  1:928124/1‑140 (MQ=255)
      aCCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 <  1:51500/140‑1 (MQ=255)
      aCCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 >  2:805846/1‑140 (MQ=255)
      aCCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 <  2:563155/140‑1 (MQ=255)
       cccGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 >  1:608010/1‑139 (MQ=255)
       cccGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 <  2:608010/139‑1 (MQ=255)
       cccGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTa                                                                                                                                <  2:63789/140‑1 (MQ=255)
              aCAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGaata                                                                                                                         >  1:362495/1‑140 (MQ=255)
                aGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 <  2:1304035/130‑1 (MQ=255)
                aGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGgtgt                                                                                                                                 >  1:1304035/1‑130 (MQ=255)
                aGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGaataat                                                                                                                       <  1:805846/140‑1 (MQ=255)
                   cccgccgTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAg                                                                                                                    <  1:24120/139‑1 (MQ=255)
                    ccgccgTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAgg                                                                                                                   <  2:409427/140‑1 (MQ=255)
                     cgccgTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGc                                                                                                                  >  2:1106789/1‑140 (MQ=255)
                      gccgTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCt                                                                                                                 >  1:815423/1‑140 (MQ=255)
                       ccgTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTg                                                                                                                <  1:752510/140‑1 (MQ=255)
                       ccgTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTg                                                                                                                <  2:915729/140‑1 (MQ=255)
                       ccgTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTg                                                                                                                <  1:1156299/140‑1 (MQ=255)
                            aGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAAt                                                                                                           <  2:1222527/140‑1 (MQ=255)
                              aaTAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCa                                                                                                         <  2:1092910/140‑1 (MQ=255)
                                  aCTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTGCGGAATAATCAGGCTGATAATCAGGAt                                                                                                     >  1:84133/1‑140 (MQ=255)
                                      cGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATcgcg                                                                                                 <  2:978443/140‑1 (MQ=255)
                                       gTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGa                                                                                                <  1:1034930/140‑1 (MQ=255)
                                       gTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGa                                                                                                <  2:164837/140‑1 (MQ=255)
                                       gTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGa                                                                                                <  2:217072/140‑1 (MQ=255)
                                       gTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGa                                                                                                <  1:823267/140‑1 (MQ=255)
                                        ttCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGaa                                                                                               >  2:272392/1‑140 (MQ=255)
                                          catcaCCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAaca                                                                                             <  1:699304/140‑1 (MQ=255)
                                          catcaCCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAaca                                                                                             <  2:288802/140‑1 (MQ=255)
                                          catcaCCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAaca                                                                                             <  2:199694/140‑1 (MQ=255)
                                           atcaCCGCCGCCGCAGAGGGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAacac                                                                                            >  2:1311555/1‑140 (MQ=255)
                                           atcaCCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGaa                                                                                               >  1:271964/1‑137 (MQ=255)
                                           atcaCCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGaa                                                                                               <  2:271964/137‑1 (MQ=255)
                                             caCCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACAt                                                                                          <  2:464500/140‑1 (MQ=255)
                                               ccgccgccgCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAg                                                                                        <  2:627137/140‑1 (MQ=255)
                                                       gcAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAg                                                                                >  2:164588/1‑140 (MQ=255)
                                                         agagagCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGa                                                                              <  1:232358/140‑1 (MQ=255)
                                                         agagagCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGa                                                                              <  1:272392/140‑1 (MQ=255)
                                                         agagagCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGa                                                                              <  2:929239/140‑1 (MQ=255)
                                                         agagagCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGa                                                                              <  2:284449/140‑1 (MQ=255)
                                                           agagCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGAtt                                                                            <  1:601208/140‑1 (MQ=255)
                                                            gggCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTc                                                                           <  1:1311555/138‑1 (MQ=255)
                                                             agCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCg                                                                          >  1:1259590/1‑140 (MQ=255)
                                                               cAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGcc                                                                        >  1:928457/1‑140 (MQ=255)
                                                                    ccGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGTTAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGaaa                                                                   >  2:288599/1‑139 (MQ=255)
                                                                      gcgTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAacac                                                                 <  1:536451/140‑1 (MQ=255)
                                                                       cgTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAg                                                                      <  1:530998/134‑1 (MQ=255)
                                                                       cgTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAg                                                                      >  2:530998/1‑134 (MQ=255)
                                                                       cgTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACAcc                                                                >  2:588001/1‑140 (MQ=255)
                                                                       cgTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACAcc                                                                >  2:1097664/1‑140 (MQ=255)
                                                                                acacCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTAcc                                                       >  1:1330823/1‑140 (MQ=255)
                                                                                 cacCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACcg                                                      <  1:19477/140‑1 (MQ=255)
                                                                                  acCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACcgc                                                     <  1:1192956/140‑1 (MQ=255)
                                                                                     aGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTTGCCAGAACACCTGACGTACcgcgcg                                                  <  1:425729/140‑1 (MQ=255)
                                                                                      gggTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGt                                                 >  2:665378/1‑140 (MQ=255)
                                                                                           aaCGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACcgcgcg                                                  >  2:561675/1‑134 (MQ=255)
                                                                                           aaCGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACcgcgcg                                                  <  1:561675/134‑1 (MQ=255)
                                                                                            aCGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATg                                           >  2:1180464/1‑140 (MQ=255)
                                                                                            aCGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATg                                           >  1:625578/1‑140 (MQ=255)
                                                                                               gACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTGCGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGttt                                        <  2:84133/140‑1 (MQ=255)
                                                                                                     cGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGcc                                  >  2:381651/1‑140 (MQ=255)
                                                                                                           tGTCTTTAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCt                            <  2:1020444/140‑1 (MQ=255)
                                                                                                               tttAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCgg                        >  1:918332/1‑140 (MQ=255)
                                                                                                                 tAACATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGAc                      <  1:54804/140‑1 (MQ=255)
                                                                                                                  aaCATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACt                     <  2:644111/140‑1 (MQ=255)
                                                                                                                  aaCATCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACt                     <  1:1106789/140‑1 (MQ=255)
                                                                                                                      tCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGcc                 <  2:362495/140‑1 (MQ=255)
                                                                                                                      tCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGcc                 <  1:502468/140‑1 (MQ=255)
                                                                                                                      tCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGcc                 <  1:1042568/140‑1 (MQ=255)
                                                                                                                      tCGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGcc                 <  2:185064/140‑1 (MQ=255)
                                                                                                                       cGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCt                >  2:513071/1‑140 (MQ=255)
                                                                                                                       cGTTACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCt                >  2:678070/1‑140 (MQ=255)
                                                                                                                         ttACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTgc              >  2:1022733/1‑140 (MQ=255)
                                                                                                                         ttACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTgc              >  2:640924/1‑140 (MQ=255)
                                                                                                                         ttACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTgc              <  1:640924/140‑1 (MQ=255)
                                                                                                                         ttACGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTgc              >  1:362464/1‑140 (MQ=255)
                                                                                                                           aCGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTgcag            <  2:1389584/140‑1 (MQ=255)
                                                                                                                           aCGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTgcag            <  1:678070/140‑1 (MQ=255)
                                                                                                                            cGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTg               >  1:802423/1‑136 (MQ=255)
                                                                                                                            cGCAGCAGGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTg               <  2:802423/136‑1 (MQ=255)
                                                                                                                              cagcagGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCaa         >  2:533436/1‑140 (MQ=255)
                                                                                                                                gcagGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATa       <  2:656488/140‑1 (MQ=255)
                                                                                                                                 cagGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATAc      >  2:413295/1‑140 (MQ=255)
                                                                                                                                 cagGCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATAc      >  1:599069/1‑140 (MQ=255)
                                                                                                                                   ggCTCGGATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATCCCg    <  1:1039309/140‑1 (MQ=255)
                                                                                                                                    gCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGa   >  2:1154227/1‑140 (MQ=255)
                                                                                                                                    gCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGa   >  2:1123129/1‑140 (MQ=255)
                                                                                                                                    gCTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGa   >  1:1270429/1‑140 (MQ=255)
                                                                                                                                     cTCGAATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGAt  <  1:44580/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GGAGTTACCCGCTGACAGCACCGCCGTCAGAATAACTGCGTTCATCACCGCCGCCGCAGAGAGCAGACCCGCGTGCTGGAACACCAGGGTGAACGGACTAACGCTGATGTCTTTAACATCGTTACGCAGCAGGCTCGGATCGGTGTACGGAATAATCAGGCTGATAATCAGGATCGCGAACACATAGAACAACAGGATTCGCCAGAACACCTGACGTACCGCGCGTGGAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGAT  >  NC_000913/2247591‑2247863

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