Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 852,861 C→G 7.4% 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,8610CG 7.4% 381.5 / 17.9 122intergenic (‑264/+8)opgE/mntSOPG biosynthetic transmembrane phosphoethanolamine transferase/Mn(2)‑response protein, MntR‑repressed
Reads supporting (aligned to +/- strand):  ref base C (50/63);  new base G (4/5);  total (54/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

Only 100 of 126 total aligned reads displayed.

GGCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCTGAC  >  NC_000913/852711‑853010
                                                                                                                                                      |                                                                                                                                                     
ggCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACcgc                                                                                                                                                       >  2:580285/1‑150 (MQ=255)
 gCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTGAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGAGTACcgc                                                                                                                                                       >  1:92265/1‑150 (MQ=255)
 gCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACcgc                                                                                                                                                       <  1:501329/150‑1 (MQ=255)
     cAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcg                                                                                                                                                   >  2:1325880/1‑150 (MQ=255)
          agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCt                                                                                                                                              <  1:373878/150‑1 (MQ=255)
          agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCt                                                                                                                                              >  1:1704997/1‑150 (MQ=255)
           gaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCtt                                                                                                                                             <  1:1498649/150‑2 (MQ=255)
            aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttc                                                                                                                                            <  2:294637/150‑3 (MQ=255)
            aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCttc                                                                                                                                            <  1:588470/150‑3 (MQ=255)
            aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAt                                                                                                                                            >  2:201010/1‑150 (MQ=255)
            aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAt                                                                                                                                            >  1:1388319/1‑150 (MQ=255)
               gagcgagcCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTa                                                                                                                                         <  1:1006968/150‑1 (MQ=255)
               gagcgagcCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTa                                                                                                                                         >  1:516861/1‑150 (MQ=255)
                 gcgagcCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATc                                                                                                                                       <  2:450840/150‑1 (MQ=255)
                 ccgagcCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGACGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATc                                                                                                                                       <  1:694020/149‑1 (MQ=255)
                  cgagcCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCg                                                                                                                                      <  1:1360057/150‑1 (MQ=255)
                      ccgccgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAg                                                                                                                                  >  2:95203/1‑150 (MQ=255)
                        gccgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGt                                                                                                                                >  2:782139/1‑150 (MQ=255)
                         ccgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGtt                                                                                                                               <  1:69619/150‑1 (MQ=255)
                          cgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGttt                                                                                                                              <  2:1704997/150‑1 (MQ=255)
                          cgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGttt                                                                                                                              <  2:651286/150‑1 (MQ=255)
                            ttCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTAt                                                                                                                            <  1:1064457/150‑1 (MQ=255)
                             tCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATc                                                                                                                           <  2:1348487/150‑1 (MQ=255)
                                ggacAACTACTGACGAACAAATCATTTGCGAAGGGGATTGCAAATACGACAAAAAGTGTGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtg                                                                                                                        <  2:1290394/147‑1 (MQ=255)
                                 gTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtgc                                                                                                                       <  1:1142972/150‑1 (MQ=255)
                                 gTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtga                                                                                                                       >  1:1385907/1‑149 (MQ=255)
                                 gTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAATTATCGGAAGGTTTATCTtgc                                                                                                                       >  1:252943/1‑150 (MQ=255)
                                  tCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtgct                                                                                                                      <  2:1721753/150‑1 (MQ=255)
                                    aaCTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtgctgc                                                                                                                    <  2:499904/150‑1 (MQ=255)
                                     actactGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCg                                                                                                                   <  2:99417/150‑1 (MQ=255)
                                           gACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCtt                                                                                                                                             >  2:1113295/1‑117 (MQ=255)
                                           gACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgt                                                                                                             <  2:1140970/150‑1 (MQ=255)
                                           gACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgt                                                                                                             >  1:1537080/1‑150 (MQ=255)
                                            aCGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgtt                                                                                                            <  1:1666446/150‑1 (MQ=255)
                                             cGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACTCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttg                                                                                                           <  1:1568734/150‑1 (MQ=255)
                                             cGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttg                                                                                                           <  2:788234/150‑1 (MQ=255)
                                              gAACAAATCTTTTGCGAAGGGGATTGCATATCCGACAAAAAGTGGGGTTTTTCAGACCTAAACGCGTTATTCCTGCATGAAAGAAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgt                                                                                                          <  1:277588/150‑1 (MQ=255)
                                               aaCAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgtt                                                                                                         >  1:466789/1‑150 (MQ=255)
                                               aaCAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgtt                                                                                                         <  2:1668195/150‑1 (MQ=255)
                                                 cAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGa                                                                                                       >  1:1656980/1‑150 (MQ=255)
                                                  aaaTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGAc                                                                                                      >  1:318446/1‑150 (MQ=255)
                                                  aaaTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGAc                                                                                                      >  1:1278167/1‑150 (MQ=255)
                                                         ttGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCg                                                                                               <  1:238365/150‑1 (MQ=255)
                                                            cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcac                                                                                            <  1:1157735/150‑1 (MQ=255)
                                                            cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcac                                                                                            <  2:529287/150‑1 (MQ=255)
                                                             gAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcaca                                                                                           >  2:830711/1‑150 (MQ=255)
                                                             gAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcaca                                                                                           >  2:1067308/1‑150 (MQ=255)
                                                             gAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcaca                                                                                           <  1:735114/150‑1 (MQ=255)
                                                                ggggATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAATGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCacaaca                                                                                        <  2:609241/150‑1 (MQ=255)
                                                                   gATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATag                                                                                     <  1:743527/150‑1 (MQ=255)
                                                                   gATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATag                                                                                     >  2:548472/1‑150 (MQ=255)
                                                                   gATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATag                                                                                     <  2:140197/150‑1 (MQ=255)
                                                                       gCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACagagagag                                                                                 >  1:947403/1‑150 (MQ=255)
                                                                       gCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagagag                                                                                 >  2:1687749/1‑150 (MQ=255)
                                                                           aTCCGACAAAACGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGagcagc                                                                             <  1:1443938/150‑1 (MQ=255)
                                                                           aTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGagcagc                                                                             <  1:1450694/150‑1 (MQ=255)
                                                                             ccGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCAt                                                                           <  1:943792/150‑1 (MQ=255)
                                                                               gACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTa                                                                         >  2:503551/1‑150 (MQ=255)
                                                                               gACAAAAACTGGGGGGTTTCAGCCCAAAACGCGTTATGCCCGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTa                                                                         <  2:261425/150‑1 (MQ=255)
                                                                                aCAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTaa                                                                        >  2:598134/1‑150 (MQ=255)
                                                                                   aaaaGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCg                                                                     >  1:347453/1‑150 (MQ=255)
                                                                                   aaaaGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCg                                                                     >  1:1326621/1‑150 (MQ=255)
                                                                                   aaaaGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCg                                                                     <  1:179249/150‑1 (MQ=255)
                                                                                   aaaaGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCg                                                                     >  1:1562993/1‑150 (MQ=255)
                                                                                   aaaaGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCg                                                                     <  1:201010/150‑1 (MQ=255)
                                                                                        gggggTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGGGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACtt                                                                <  2:1406545/149‑1 (MQ=255)
                                                                                         ggggTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAATCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACttt                                                               <  2:1537080/150‑1 (MQ=255)
                                                                                                 tCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGAGCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgat                                                       >  2:28264/1‑149 (MQ=255)
                                                                                                 tCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgag                                                       <  1:580285/150‑1 (MQ=255)
                                                                                                    gACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAAt                                                    >  1:1430985/1‑150 (MQ=255)
                                                                                                    gACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAAt                                                    <  2:466789/150‑1 (MQ=255)
                                                                                                      ccAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGa                                                  >  1:671735/1‑150 (MQ=255)
                                                                                                        aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                                >  2:1290721/1‑150 (MQ=255)
                                                                                                          aaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTaa                                              >  1:1658772/1‑150 (MQ=255)
                                                                                                          aaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTaa                                              >  2:528026/1‑150 (MQ=255)
                                                                                                            cgcgTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGACGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAac                                            >  2:141562/1‑150 (MQ=255)
                                                                                                                 tATGCCTGCATGAAATCAGGAGTGAAGCGGCGTACCGCGTCGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgcg                                       <  1:357520/150‑1 (MQ=255)
                                                                                                                 tATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgcg                                       >  2:1615465/1‑150 (MQ=255)
                                                                                                                   tGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCa                                     >  1:731352/1‑150 (MQ=255)
                                                                                                                   tGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATATATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGAGCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCa                                     <  2:670341/150‑1 (MQ=255)
                                                                                                                      cTGCATGAAAGCAGGAGTGAAGCGGGGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCGGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATac                                  <  2:92265/150‑1 (MQ=255)
                                                                                                                         cATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACAcc                               <  1:1226318/150‑1 (MQ=255)
                                                                                                                         cATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACAcc                               <  1:1619119/150‑1 (MQ=255)
                                                                                                                          aTGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACct                              <  2:1359786/150‑1 (MQ=255)
                                                                                                                            gAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACctct                            >  2:1359045/1‑150 (MQ=255)
                                                                                                                              aaGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTg                          <  2:1656980/150‑1 (MQ=255)
                                                                                                                                 cAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAAc                       <  1:1467664/150‑1 (MQ=255)
                                                                                                                                     taaataGCGGCGTCGCGGTACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCTCAACATAGAGAGCAGCAGTAACCGTACTTTTAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcat                   >  1:720589/6‑150 (MQ=255)
                                                                                                                                      gTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATTCACCTCTTGAACtcatt                  >  1:1005791/1‑150 (MQ=255)
                                                                                                                                      gTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcatt                  >  1:687426/1‑150 (MQ=255)
                                                                                                                                      gTGAAGCGGCGTACCGCGACGCCGCGATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCAGAACATAGAGAGCCGCATTACCCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGATCtcatt                  >  1:1080160/1‑150 (MQ=255)
                                                                                                                                       tGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACCCCTCTTGAACtcattc                 >  2:1274404/1‑150 (MQ=255)
                                                                                                                                        gAAGCGGCGTACCGCGACGCCGCTATTTATCGTAATGTTTACCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcattca                <  2:653058/150‑1 (MQ=255)
                                                                                                                                        gAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcattca                <  2:318446/150‑1 (MQ=255)
                                                                                                                                        ataGCGGCGTCGCGGTACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcattca                >  2:1156244/3‑150 (MQ=255)
                                                                                                                                          aGCGGCGTCGCGGTACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATa              >  1:512097/1‑150 (MQ=255)
                                                                                                                                          aGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATa              <  1:598134/150‑1 (MQ=255)
                                                                                                                                              gcgTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAc          <  1:1724392/150‑1 (MQ=255)
                                                                                                                                               cgTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcc         <  1:95203/150‑1 (MQ=255)
                                                                                                                                                      cgACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCTGAc  <  2:1178047/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCTGAC  >  NC_000913/852711‑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: