Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 852,862 G→T 10.0% intergenic (‑265/+7) 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,8620GT10.0% 303.1 / 20.0 110intergenic (‑265/+7)opgE/mntSOPG biosynthetic transmembrane phosphoethanolamine transferase/Mn(2)‑response protein, MntR‑repressed
Reads supporting (aligned to +/- strand):  ref base G (58/41);  new base T (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 = 1.11e-01

Only 100 of 111 total aligned reads displayed.

TGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCTGAC  >  NC_000913/852714‑853010
                                                                                                                                                    |                                                                                                                                                    
    agaagaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGCCCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcgcc                                                                                                                                                 >  2:1522346/1‑150 (MQ=255)
       agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCt                                                                                                                                              <  2:503167/150‑1 (MQ=255)
       agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCt                                                                                                                                              <  2:1171213/150‑1 (MQ=255)
       agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCt                                                                                                                                              <  1:629394/150‑1 (MQ=255)
        gaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCtt                                                                                                                                             <  2:713234/150‑2 (MQ=255)
         aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttc                                                                                                                                            <  2:1487149/150‑3 (MQ=255)
         aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttc                                                                                                                                            <  2:569056/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                                                                                                                                           <  1:1309127/140‑1 (MQ=255)
                    cgccgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAtt                                                                                                                                           >  2:1309127/1‑140 (MQ=255)
                          tCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTAta                                                                                                                           >  2:1090561/1‑149 (MQ=255)
                           cAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCt                                                                                                                          >  2:532662/1‑150 (MQ=255)
                            aGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCtt                                                                                                                         >  1:1441888/1‑150 (MQ=255)
                             ggTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtg                                                                                                                        >  1:185960/1‑150 (MQ=255)
                             ggTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtg                                                                                                                        >  2:943838/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:570084/1‑150 (MQ=255)
                                     actGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGtt                                                                                                                <  2:185289/150‑1 (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)
                                                aaTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGAcc                                                                                                     <  1:827846/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)
                                                              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)
                                                                    gCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCATACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagagag                                                                                 >  1:1478448/1‑150 (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)
                                                                                   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)
                                                                                               cAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgaga                                                      >  2:1481465/1‑150 (MQ=255)
                                                                                                 gACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCACATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAAt                                                    >  1:802373/1‑150 (MQ=255)
                                                                                                   ccAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATg                                                   >  1:499958/1‑149 (MQ=255)
                                                                                                   ccAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGa                                                  >  1:1677657/1‑150 (MQ=255)
                                                                                                     aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                                >  1:1598770/1‑150 (MQ=255)
                                                                                                     aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                                <  1:532662/150‑1 (MQ=255)
                                                                                                     aaaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACt                                                <  1:562551/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)
                                                                                                                  ccTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATa                                   >  2:922348/1‑150 (MQ=255)
                                                                                                                    tGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATaca                                 >  2:1032458/1‑150 (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)
                                                                                                                             gCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGaa                        >  1:1250479/1‑150 (MQ=255)
                                                                                                                             gCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGaa                        >  2:101924/1‑150 (MQ=255)
                                                                                                                              cAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAAc                       <  1:570084/150‑1 (MQ=255)
                                                                                                                              cAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAAc                       <  1:388380/150‑1 (MQ=255)
                                                                                                                                     gAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGTGAATGACTAAACACGCGCATACACCTCTTGAACtcattca                >  1:66860/1‑150 (MQ=255)
                                                                                                                                       aGCGGCGTCGCGGTACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATa              >  2:710707/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)
                                                                                                                                         cggcgTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAg            <  2:1250479/150‑1 (MQ=255)
                                                                                                                                          ggcgTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGa           >  2:1066956/1‑150 (MQ=255)
                                                                                                                                                ccGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcctcct     >  1:385461/1‑150 (MQ=255)
                                                                                                                                                ccGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcctcct     <  2:1221091/150‑1 (MQ=255)
                                                                                                                                                   cgACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTAATAAGACCTCCTGAc  >  2:1593581/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCCTGAC  >  NC_000913/852714‑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: