Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 852,857 A→C 5.0% intergenic (‑260/+12) 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,8570AC 5.0% 377.2 / 6.1 130intergenic (‑260/+12)opgE/mntSOPG biosynthetic transmembrane phosphoethanolamine transferase/Mn(2)‑response protein, MntR‑repressed
Reads supporting (aligned to +/- strand):  ref base A (67/56);  new base C (3/4);  total (70/60)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 132 total aligned reads displayed.

GCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCC  >  NC_000913/852712‑853006
                                                                                                                                                 |                                                                                                                                                     
gCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACcgc                                                                                                                                                   >  2:847096/1‑150 (MQ=255)
gCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACcgc                                                                                                                                                   >  1:1924007/1‑150 (MQ=255)
 cTGCAATAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGGTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACcgcg                                                                                                                                                  <  2:1294463/150‑1 (MQ=255)
 cTGCAAGAAGATAGAGCGCACCGCCGTTCAGGTCAACTACTGACGCACAAATCTTTTGCCAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTATTCAGCACAAACCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCTGCGTACcgcg                                                                                                                                                  >  2:1390352/1‑150 (MQ=255)
    cAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACTAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcg                                                                                                                                               >  2:1908686/1‑150 (MQ=255)
    cAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcg                                                                                                                                               <  2:1479258/150‑1 (MQ=255)
    cAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAACCGGCGTACcgcgccc                                                                                                                                               >  2:1262989/1‑147 (MQ=255)
      agaagaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAAGCCGACAAAAAGTGGGGTTTTTCAGACCAAACCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcgcc                                                                                                                                             >  1:855042/1‑150 (MQ=255)
       gaagaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcgccg                                                                                                                                            >  2:1083341/1‑150 (MQ=255)
        aagaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcgccgc                                                                                                                                           <  2:1404690/150‑1 (MQ=255)
        aagaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGAcgccgc                                                                                                                                           <  1:507762/150‑1 (MQ=255)
         agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCt                                                                                                                                          <  1:1663115/150‑1 (MQ=255)
         agaTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCt                                                                                                                                          >  1:2048380/1‑150 (MQ=255)
           aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttc                                                                                                                                        <  1:1044883/150‑3 (MQ=255)
           aTAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAt                                                                                                                                        <  1:929598/150‑1 (MQ=255)
            tAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttca                                                                                                                                       <  2:290505/150‑4 (MQ=255)
            tAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAtt                                                                                                                                       <  2:1195918/150‑1 (MQ=255)
             agagCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTAttt                                                                                                                                      >  1:1930654/1‑150 (MQ=255)
               agcgagcCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCGCGGTACGCCGCttcactc                                                                                                                                    <  2:888345/150‑7 (MQ=255)
               agcgagcCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTAt                                                                                                                                    >  2:1240614/1‑150 (MQ=255)
                    gccgccGTTCAGGTCACCTCCTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAATAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCAGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGACGCAATTTATCGGaa                                                                                                                               >  1:1321424/1‑150 (MQ=255)
                    gccgccGTTCAGGTCAAATACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGAGGGGTTTTTCAGACCAAAACGCGTTATTCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGaa                                                                                                                               >  1:461499/1‑150 (MQ=255)
                             cAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTCCCGCGACGCCGCTATTTATCGGAAGGTTTATCt                                                                                                                      >  2:1976239/1‑150 (MQ=255)
                             cAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCt                                                                                                                      >  1:1379485/1‑150 (MQ=255)
                             cAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGATTATCt                                                                                                                      <  2:822729/150‑1 (MQ=255)
                             cAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGCAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCt                                                                                                                      >  1:58764/1‑150 (MQ=255)
                              aGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCtt                                                                                                                     >  1:1254477/1‑150 (MQ=255)
                               ggTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtg                                                                                                                    >  2:760846/1‑150 (MQ=255)
                                gTCAACTACTTACGAACAAATCTTTTTCGAAGGGGATGGCAAATCCGACAAATAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGATAGAAGCGGCGTACCGCGACGCCTCTATTTATCGGAAGGTTTATCTtgc                                                                                                                   <  1:102966/150‑1 (MQ=255)
                                gTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtgc                                                                                                                   <  1:1009317/150‑1 (MQ=255)
                                  cAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTtgctg                                                                                                                 >  2:284045/1‑150 (MQ=255)
                                    actactGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACTCCGCTATTTATCGGAAGGGTTATCTTGATGCg                                                                                                               >  2:390720/1‑150 (MQ=255)
                                    actactGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCg                                                                                                               >  2:390718/1‑150 (MQ=255)
                                     ctactGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGTCGCTATTTATCGGAAGGTTTATCTTGCTGCgg                                                                                                              >  2:451702/1‑150 (MQ=255)
                                       actGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGtt                                                                                                            <  1:1795952/150‑1 (MQ=255)
                                          gACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgt                                                                                                         >  1:2035679/1‑150 (MQ=255)
                                          gACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgt                                                                                                         >  2:2073295/1‑150 (MQ=255)
                                           aCGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGGGGGGTGTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgtt                                                                                                        <  1:176921/150‑1 (MQ=255)
                                            cGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGATGCGGTttgttg                                                                                                       >  1:839819/1‑150 (MQ=255)
                                               aCAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                    <  1:1262708/150‑1 (MQ=255)
                                               aCAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgttg                                                                                                    <  2:1963223/150‑1 (MQ=255)
                                                  aaTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGGGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGGTTATCTTGCTGCGGTTTGTTGTTGAcc                                                                                                 <  1:1182920/150‑1 (MQ=255)
                                                   atattttGCGAATGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCa                                                                                                <  2:1832591/147‑1 (MQ=255)
                                                   aTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGTGGTTTTTCATACCAAAAAGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATAGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCa                                                                                                <  1:1391315/150‑1 (MQ=255)
                                                   aTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGATTGTTGTTGACCa                                                                                                >  1:106883/1‑150 (MQ=255)
                                                     cTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCata                                                                                              >  1:954703/1‑150 (MQ=255)
                                                       tttGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACGAAAACGCGTTATGCCTGCATCAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGGTATTTAGCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATc                                                                                            >  1:1451045/1‑150 (MQ=255)
                                                       tttGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGCCGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATc                                                                                            >  1:607579/1‑150 (MQ=255)
                                                         tGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGc                                                                                          >  1:1001535/1‑150 (MQ=255)
                                                           cGAAGGGGAGTGCAAATCCGACAAAATGTGGGGTTTTTCAGACCAAAAGGCGTTATGCCTGCTTGATGGCAGGAGTGAAGCGGCCTACCGCGACGCCGCAAGTTATCGGAACGTTTATCTTTCTGCGGTTGGTTGTTGACATTATCGcac                                                                                        >  2:1280687/1‑150 (MQ=255)
                                                              aGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCacaac                                                                                     <  1:1943269/150‑1 (MQ=255)
                                                               ggggATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCacaaca                                                                                    <  1:391497/150‑1 (MQ=255)
                                                                    ttGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagag                                                                               >  1:1876207/1‑150 (MQ=255)
                                                                     tGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagaga                                                                              <  2:2054048/150‑1 (MQ=255)
                                                                      gCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagagcg                                                                             >  2:850375/1‑148 (MQ=255)
                                                                      gCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagagag                                                                             >  1:1661604/1‑150 (MQ=255)
                                                                       cAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGagc                                                                            <  1:173144/150‑1 (MQ=255)
                                                                         aaTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGagcag                                                                          >  1:1451009/1‑150 (MQ=255)
                                                                           tCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAgcagca                                                                        >  2:1334230/1‑150 (MQ=255)
                                                                             cGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCAtt                                                                      >  2:882576/1‑150 (MQ=255)
                                                                              gACAAAAAGCGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTa                                                                     <  2:103207/150‑1 (MQ=255)
                                                                                      gTGGGGTGTTTCAGACCAAAACGCTTTATGCCTGCATGAAAGCAGGAGTGACGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACt                                                             <  2:855042/150‑1 (MQ=255)
                                                                                           gTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCAGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTaaa                                                        >  1:1876081/1‑150 (MQ=255)
                                                                                            tttttCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAg                                                       <  2:175689/150‑1 (MQ=255)
                                                                                              tttCAGACCAAAACGCGTTATGCCTGCATGAAATCAGGAGGGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGGTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAggg                                                     <  2:801796/150‑1 (MQ=255)
                                                                                                tCAGACCAATACGCGTTATGCCTGCATGAAAGCTGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgag                                                   <  2:128861/150‑1 (MQ=255)
                                                                                                tCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgag                                                   >  1:759112/1‑150 (MQ=255)
                                                                                                 cAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgaga                                                  <  1:983824/150‑1 (MQ=255)
                                                                                                   gACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGGAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAAt                                                >  1:725490/1‑150 (MQ=255)
                                                                                                    aCCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCTTATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATg                                               <  2:1063556/150‑1 (MQ=255)
                                                                                                    aCCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATg                                               <  2:1254477/150‑1 (MQ=255)
                                                                                                    aCCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATg                                               >  1:471778/1‑150 (MQ=255)
                                                                                                     ccAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGa                                              >  1:1891461/1‑150 (MQ=255)
                                                                                                     ccAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGa                                              >  1:1478169/1‑150 (MQ=255)
                                                                                                               gtATGCCTGCTTGAACGCCGGAGTGACGCGGCGTACCGCGACGCCGCTATTTATCGTAAGGTTTATCTTGCGGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgc                                    <  2:1399185/149‑1 (MQ=255)
                                                                                                                tATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgcg                                   >  2:1158791/1‑150 (MQ=255)
                                                                                                                 aTGCCTGCATGAATGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTTTTGTAGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACgcgc                                  <  1:390720/150‑1 (MQ=255)
                                                                                                                 aTGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgcgg                                  <  1:390718/150‑2 (MQ=255)
                                                                                                                    ccTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATa                               >  2:470005/1‑150 (MQ=255)
                                                                                                                    ccTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATa                               >  2:107095/1‑150 (MQ=255)
                                                                                                                    ccTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATa                               >  1:1633761/1‑150 (MQ=255)
                                                                                                                      tGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATaca                             >  1:1399855/1‑150 (MQ=255)
                                                                                                                      tGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATaca                             >  2:600037/1‑150 (MQ=255)
                                                                                                                           gAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTACACACGCGCATACACctct                        >  1:583844/1‑150 (MQ=255)
                                                                                                                           gAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACctct                        >  1:767317/1‑150 (MQ=255)
                                                                                                                           gAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACctct                        <  2:1379485/150‑1 (MQ=255)
                                                                                                                             aaGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTg                      <  2:284850/150‑1 (MQ=255)
                                                                                                                             aaGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTg                      >  1:1863354/1‑150 (MQ=255)
                                                                                                                              aGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGa                     <  1:2030110/150‑1 (MQ=255)
                                                                                                                               gCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGGTTGTTGTTGACCATATCGCACGACATAGAGAGCAGCATTAACCGTAATGGAAAGGGAGAATGAATAAACACGCGCATACACCTCTGGaa                    >  1:991139/1‑150 (MQ=255)
                                                                                                                                cAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAAc                   >  2:476208/1‑150 (MQ=255)
                                                                                                                                 aGGAGTGAAGCGGCGGACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACt                  <  2:1539085/150‑1 (MQ=255)
                                                                                                                                      tGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcattc             <  1:516164/150‑1 (MQ=255)
                                                                                                                                      aataGCGGCGTCGCGGTACGCCGCTATTTATCGGATGCTTTATCTTGCTGCGGTTTGTTGTTGACCATATGGCACAACATAGAGAGCACCATTCACCGTACTTTAAAGGTATAATGACTAAACACGCGCATACCCCACTTGAAGtcattc             >  2:429027/4‑150 (MQ=255)
                                                                                                                                      aataGCGGCGTCGCGGTACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcattc             >  1:522259/4‑150 (MQ=255)
                                                                                                                                       gAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcattca            <  1:847096/150‑1 (MQ=255)
                                                                                                                                       gAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcattca            <  2:2093750/150‑1 (MQ=255)
                                                                                                                                       gAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcattca            >  2:1741395/1‑150 (MQ=255)
                                                                                                                                              cgTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcc     <  1:1083341/150‑1 (MQ=255)
                                                                                                                                                 aCCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGAcctcc  >  1:479589/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCTGCAAGAAGATAGAGCGAGCCGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAAGACCTCC  >  NC_000913/852712‑853006

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