Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 852,858 C→G 10.0% intergenic (‑261/+11) 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,8580CG10.0% 339.2 / 25.0 110intergenic (‑261/+11)opgE/mntSOPG biosynthetic transmembrane phosphoethanolamine transferase/Mn(2)‑response protein, MntR‑repressed
Reads supporting (aligned to +/- strand):  ref base C (58/41);  new base G (2/9);  total (60/50)
Fisher's exact test for biased strand distribution p-value = 2.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-02

Only 100 of 112 total aligned reads displayed.

GGCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCT  >  NC_000913/852711‑853007
                                                                                                                                                   |                                                                                                                                                     
ggCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACcg                                                                                                                                                     <  2:595862/150‑1 (MQ=255)
 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                                                                                                                                        <  1:1021256/150‑4 (MQ=255)
             tAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttca                                                                                                                                        <  2:1270690/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)
                   gagcCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCgg                                                                                                                                  <  2:522988/150‑1 (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)
                               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)
                                          tGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttg                                                                                                           <  2:185960/150‑1 (MQ=255)
                                              gAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAACACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGCCGCCGCTATTTATCGGAGGGTTTATCTTGCTGCGGTttgttgt                                                                                                       >  2:310621/1‑150 (MQ=255)
                                              gAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgt                                                                                                       >  2:1506240/1‑150 (MQ=255)
                                                aCAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                     <  2:858529/150‑1 (MQ=255)
                                                aCAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                     <  2:1534711/150‑1 (MQ=255)
                                                    aTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCa                                                                                                 <  1:989725/150‑1 (MQ=255)
                                                            cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcac                                                                                         <  2:631887/150‑1 (MQ=255)
                                                            cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGc                                                                                           >  1:670189/1‑148 (MQ=255)
                                                            cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGc                                                                                           <  2:670189/148‑1 (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                                                                                                     <  1:359565/125‑1 (MQ=255)
                                                                         aaaTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                     >  2:359565/1‑125 (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:616165/1‑150 (MQ=255)
                                                                                 cAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAAc                                                                    >  1:1267598/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                                                           >  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:532662/150‑1 (MQ=255)
                                                                                                        aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                             <  1:562551/150‑1 (MQ=255)
                                                                                                        aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                             >  1:1598770/1‑150 (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                                    >  2:770768/1‑150 (MQ=255)
                                                                                                                 tATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgcg                                    <  1:716431/150‑1 (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                          <  2:1441888/150‑1 (MQ=255)
                                                                                                                           tGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACctc                          >  1:425533/1‑150 (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                     >  2:101924/1‑150 (MQ=255)
                                                                                                                                gCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGaa                     >  1:1250479/1‑150 (MQ=255)
                                                                                                                                 cAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAAc                    <  1:570084/150‑1 (MQ=255)
                                                                                                                                        gAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGTGAATGACTAAACACGCGCATACACCTCTTGAACtcattca             >  1:66860/1‑150 (MQ=255)
                                                                                                                                          aGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATa           >  2:1145264/1‑150 (MQ=255)
                                                                                                                                          aGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATa           <  2:1267598/150‑1 (MQ=255)
                                                                                                                                           gcggcgTCGCGGTACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATaa          >  1:1632159/1‑150 (MQ=255)
                                                                                                                                                   ccGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcctcct  <  2:1221091/150‑1 (MQ=255)
                                                                                                                                                   ccGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcctcct  >  1:385461/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCT  >  NC_000913/852711‑853007

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: