Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 852,861 C→G 10.0% intergenic (‑264/+8) opgE ← / ← mntS OPG biosynthetic transmembrane phosphoethanolamine transferase/Mn(2)‑response protein, MntR‑repressed

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913852,8610CG10.0% 325.4 / 24.9 110intergenic (‑264/+8)opgE/mntSOPG biosynthetic transmembrane phosphoethanolamine transferase/Mn(2)‑response protein, MntR‑repressed
Reads supporting (aligned to +/- strand):  ref base C (57/41);  new base G (2/9);  total (60/50)
Fisher's exact test for biased strand distribution p-value = 2.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

Only 100 of 112 total aligned reads displayed.

GCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCTGAC  >  NC_000913/852712‑853010
                                                                                                                                                     |                                                                                                                                                     
gCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACcgc                                                                                                                                                       >  2:1301375/1‑150 (MQ=255)
  tGCAAGAAGATAGCGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCCCGCGa                                                                                                                                                     >  2:1154856/1‑150 (MQ=255)
      agaagaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGCCCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcgcc                                                                                                                                                 >  2:1522346/1‑150 (MQ=255)
         agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCt                                                                                                                                              <  2:1171213/150‑1 (MQ=255)
         agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCt                                                                                                                                              <  2:503167/150‑1 (MQ=255)
          gaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCtt                                                                                                                                             <  2:713234/150‑2 (MQ=255)
           aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttc                                                                                                                                            <  2:569056/150‑3 (MQ=255)
           aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttc                                                                                                                                            <  2:1487149/150‑3 (MQ=255)
           aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTGCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttc                                                                                                                                            <  1:346536/150‑3 (MQ=255)
            tAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttca                                                                                                                                           <  2:1270690/150‑4 (MQ=255)
            tAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttca                                                                                                                                           <  1:1021256/150‑4 (MQ=255)
            tAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAtt                                                                                                                                           <  1:956400/150‑1 (MQ=255)
                 cgagcCGCCGTACACGTCAACTACTGACGAACAACTCATTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGCAGGGAAGCGGCGTACCGCGACGCCGCTATTTATCg                                                                                                                                      <  1:1504959/150‑1 (MQ=255)
                  tagcCGCCGTTCAGGTCATCTACTGTCGAACAAATCTTTTGCGAAGGGGTTTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAACCGCTTTATGACTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATctg                                                                                                                                     <  2:1283117/149‑3 (MQ=255)
                    gccgccGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGaa                                                                                                                                   >  2:334974/1‑150 (MQ=255)
                    gccgccGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGaa                                                                                                                                   >  1:889266/1‑150 (MQ=255)
                      cgccgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAtt                                                                                                                                           >  2:1309127/1‑140 (MQ=255)
                      cgccgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAtt                                                                                                                                           <  1:1309127/140‑1 (MQ=255)
                        ccgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGtt                                                                                                                               >  2:1720202/1‑150 (MQ=255)
                            tCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTAta                                                                                                                           >  2:1090561/1‑149 (MQ=255)
                             cAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCt                                                                                                                          >  2:532662/1‑150 (MQ=255)
                              aGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCtt                                                                                                                         >  1:1441888/1‑150 (MQ=255)
                               ggTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtg                                                                                                                        >  2:943838/1‑150 (MQ=255)
                               ggTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtg                                                                                                                        >  1:185960/1‑150 (MQ=255)
                                gTCAACTACTGACGAACAAATCTTTAGCGAAGGGGATTGCAAATACGACAAAAAGTGGGGTTATTCAGACCAAAACGCGTTATGCCAGCATGAAAGCAGGAGTGAAGCAGAGAACCGAGACGACGCTATTTATAGGAAGGATTATCTtga                                                                                                                       >  1:695830/1‑149 (MQ=255)
                                  cAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtgctg                                                                                                                     <  1:1463785/150‑1 (MQ=255)
                                  cAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTAGCGGAAGGTTTATCTtgctg                                                                                                                     >  2:1160014/1‑150 (MQ=255)
                                       actGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGtt                                                                                                                <  2:185289/150‑1 (MQ=255)
                                       actGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGtt                                                                                                                >  2:570084/1‑150 (MQ=255)
                                        ctGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGttt                                                                                                               >  1:1627416/1‑150 (MQ=255)
                                             gAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAACACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGCCGCCGCTATTTATCGGAGGGTTTATCTTGCTGCGGTttgttgt                                                                                                          >  2:310621/1‑150 (MQ=255)
                                             gAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgt                                                                                                          >  2:1506240/1‑150 (MQ=255)
                                               aCAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                        <  2:1534711/150‑1 (MQ=255)
                                               aCAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                        <  2:858529/150‑1 (MQ=255)
                                                   aTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCa                                                                                                    <  1:989725/150‑1 (MQ=255)
                                                           cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcac                                                                                            <  2:631887/150‑1 (MQ=255)
                                                           cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGc                                                                                              <  2:670189/148‑1 (MQ=255)
                                                           cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGc                                                                                              >  1:670189/1‑148 (MQ=255)
                                                            gAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTTCTGCGGTTTGTTGTTGACCATATCGaaca                                                                                           >  2:216961/1‑150 (MQ=255)
                                                            gAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcaca                                                                                           >  2:681615/1‑150 (MQ=255)
                                                                gggATTGCAAATCCGACAAAAAGGGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACAt                                                                                       <  2:1396721/150‑1 (MQ=255)
                                                                 ggATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATa                                                                                      <  2:875053/150‑1 (MQ=255)
                                                                    ttGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTTTGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagag                                                                                   <  1:334974/150‑1 (MQ=255)
                                                                    ttGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGCCCATATCGCACAACATagag                                                                                   >  1:859815/1‑150 (MQ=255)
                                                                    ttGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagag                                                                                   >  1:628205/1‑150 (MQ=255)
                                                                    ttGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagag                                                                                   <  1:1720202/150‑1 (MQ=255)
                                                                       cAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGagc                                                                                >  2:818962/1‑150 (MQ=255)
                                                                        aaaTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                        >  2:359565/1‑125 (MQ=255)
                                                                        aaaTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                        <  1:359565/125‑1 (MQ=255)
                                                                              gACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTa                                                                         >  1:1542086/1‑150 (MQ=255)
                                                                              gACAAAAAGTGGGGGTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTa                                                                         <  1:1310788/150‑1 (MQ=255)
                                                                               aCAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTaa                                                                        >  2:598232/1‑150 (MQ=255)
                                                                                cAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACTTAGAGAGAAGCATTaaa                                                                       >  2:1095515/1‑149 (MQ=255)
                                                                                cAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAAc                                                                       >  1:1267598/1‑150 (MQ=255)
                                                                                cAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAAc                                                                       >  1:616165/1‑150 (MQ=255)
                                                                                 aaaaaGTGGGGTGTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAAcc                                                                      <  1:1496595/150‑1 (MQ=255)
                                                                                     aGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTAc                                                                  >  1:196614/1‑150 (MQ=255)
                                                                                      gTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACt                                                                 >  1:1221091/1‑150 (MQ=255)
                                                                                         gggTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTa                                                              <  2:191187/150‑1 (MQ=255)
                                                                                         gggTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTa                                                              >  1:777481/1‑150 (MQ=255)
                                                                                          ggTTTTTCAGACAAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTTTTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTaa                                                             >  1:1532279/1‑150 (MQ=255)
                                                                                           gTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTaaa                                                            >  2:388380/1‑150 (MQ=255)
                                                                                            tttttCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAg                                                           <  1:1301375/150‑1 (MQ=255)
                                                                                            gtttgCAGACCATAACCCGTTATGGCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTTTTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATTTCGCACAACATAGAGAGCAGCATTAACCGGACTTTAAAg                                                           <  1:842762/149‑1 (MQ=255)
                                                                                              tttCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAggg                                                         >  1:704591/1‑150 (MQ=255)
                                                                                                 cAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACAATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgaga                                                      >  1:929718/1‑150 (MQ=255)
                                                                                                   gACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCACATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAAt                                                    >  1:802373/1‑150 (MQ=255)
                                                                                                     ccAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATg                                                   <  2:499958/149‑1 (MQ=255)
                                                                                                     ccAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATg                                                   >  1:499958/1‑149 (MQ=255)
                                                                                                       aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                                <  1:562551/150‑1 (MQ=255)
                                                                                                       aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                                >  1:1598770/1‑150 (MQ=255)
                                                                                                       aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                                <  1:532662/150‑1 (MQ=255)
                                                                                                        aaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTa                                               >  2:378652/1‑150 (MQ=255)
                                                                                                          aCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTaaa                                             <  1:411667/150‑1 (MQ=255)
                                                                                                           cgcgTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAac                                            <  2:663005/150‑1 (MQ=255)
                                                                                                           cgcgTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGATATTTATCGGAAGGTATATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCTTACTTTAAAGGGAGAATGACTAAac                                            >  2:1141745/1‑150 (MQ=255)
                                                                                                            gcgTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAACAGCATTAACCGTACTTTAAAGGGAGAATGACTAAaca                                           >  2:1250299/1‑150 (MQ=255)
                                                                                                                tATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgcg                                       <  1:716431/150‑1 (MQ=255)
                                                                                                                tATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgcg                                       >  2:770768/1‑150 (MQ=255)
                                                                                                                 aTGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACgcgc                                      <  1:1522346/150‑1 (MQ=255)
                                                                                                                   gCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCAt                                    <  1:475423/150‑1 (MQ=255)
                                                                                                                    ccTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATa                                   >  2:922348/1‑150 (MQ=255)
                                                                                                                      tGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATaca                                 >  2:1032458/1‑150 (MQ=255)
                                                                                                                      tGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATaca                                 <  1:681615/150‑1 (MQ=255)
                                                                                                                          tGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACctc                             >  1:425533/1‑150 (MQ=255)
                                                                                                                          tGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACctc                             <  2:1441888/150‑1 (MQ=255)
                                                                                                                          tGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCAGTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACctc                             >  1:170728/1‑150 (MQ=255)
                                                                                                                             aaGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTg                          <  2:1542086/150‑1 (MQ=255)
                                                                                                                              aGCAGGATTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATATAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGa                         <  2:889266/150‑1 (MQ=255)
                                                                                                                               gCAGGAGTGAAGCGTCGTACCGCGACGCCGCTATTTATCTGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGATAATGACTAAACACGCGCATACACCACTTGaa                        >  2:1314472/1‑150 (MQ=255)
                                                                                                                               gCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGaa                        >  1:1250479/1‑150 (MQ=255)
                                                                                                                               gCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGaa                        >  2:101924/1‑150 (MQ=255)
                                                                                                                                cAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAAc                       <  1:570084/150‑1 (MQ=255)
                                                                                                                                       gAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGTGAATGACTAAACACGCGCATACACCTCTTGAACtcattca                >  1:66860/1‑150 (MQ=255)
                                                                                                                                         aGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATa              <  2:1267598/150‑1 (MQ=255)
                                                                                                                                         aGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATa              >  2:1145264/1‑150 (MQ=255)
                                                                                                                                          gcggcgTCGCGGTACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATaa             >  1:1632159/1‑150 (MQ=255)
                                                                                                                                                  ccGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcctcct     <  2:1221091/150‑1 (MQ=255)
                                                                                                                                                  ccGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcctcct     >  1:385461/1‑150 (MQ=255)
                                                                                                                                                     cgACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTAATAAGACCTCCTGAc  >  2:1593581/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCTGAC  >  NC_000913/852712‑853010

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