Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 691,390 G→C 8.1% A132G (GCC→GGC)  ybeX ← CorC‑HlyC family protein YbeX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913691,3900GC 8.1% 283.7 / 13.3 86A132G (GCC→GGC) ybeXCorC‑HlyC family protein YbeX
Reads supporting (aligned to +/- strand):  ref base G (36/43);  new base C (3/4);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

CAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTC  >  NC_000913/691241‑691539
                                                                                                                                                     |                                                                                                                                                     
cAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAgg                                                                                                                                                       >  2:1193334/1‑150 (MQ=255)
cAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAgg                                                                                                                                                       >  2:738504/1‑150 (MQ=255)
 aaTGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGc                                                                                                                                                      <  2:723960/150‑1 (MQ=255)
   tGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCtt                                                                                                                                                    >  1:776298/1‑150 (MQ=255)
    ggTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTAAGCATACGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTc                                                                                                                                                   >  2:759699/1‑150 (MQ=255)
     gTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                  <  1:1291789/150‑1 (MQ=255)
     gTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                  >  2:356416/1‑150 (MQ=255)
      tCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAg                                                                                                                                                 <  2:41164/150‑1 (MQ=255)
       cACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGc                                                                                                                                                >  2:1074190/1‑150 (MQ=255)
         ccAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCAt                                                                                                                                              >  1:64679/1‑150 (MQ=255)
          cAGACCGGAAACCCCACCGAATTCGTCATTAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTc                                                                                                                                                   <  2:581426/144‑1 (MQ=255)
          cAGACCGGAAACCCCACCGAATTCGTCATTAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTc                                                                                                                                                   >  1:581426/1‑144 (MQ=255)
           aGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgg                                                                                                                                            <  2:485547/150‑3 (MQ=255)
           aGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgg                                                                                                                                            <  1:770415/150‑3 (MQ=255)
            gACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgga                                                                                                                                           <  1:264875/150‑4 (MQ=255)
             aCCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtggac                                                                                                                                          <  1:914500/150‑5 (MQ=255)
              ccGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTg                                                                                                                                         <  2:322992/150‑1 (MQ=255)
                ggAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcg                                                                                                                                       >  1:273378/1‑150 (MQ=255)
                 gAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcgc                                                                                                                                      <  1:1245238/150‑1 (MQ=255)
                  aaaCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCa                                                                                                                                     <  1:1596397/150‑1 (MQ=255)
                  aaaCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCa                                                                                                                                     <  2:663165/150‑1 (MQ=255)
                      cccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATaaa                                                                                                                                 <  1:1650894/150‑1 (MQ=255)
                       ccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAAc                                                                                                                                >  1:611438/1‑150 (MQ=255)
                         aCCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACgg                                                                                                                              <  1:1074190/150‑1 (MQ=255)
                          ccGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgc                                                                                                                             <  1:1193334/150‑1 (MQ=255)
                            gAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCCCTGCGCATAAAAGgcag                                                                                                                           >  1:1062577/1‑150 (MQ=255)
                                 cGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGt                                                                                                                      <  2:1442410/150‑1 (MQ=255)
                                 cGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGt                                                                                                                      >  1:1475425/1‑150 (MQ=255)
                                    cAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCtt                                                                                                                   <  2:225347/150‑1 (MQ=255)
                                       tAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCg                                                                                                                >  1:1445328/1‑150 (MQ=255)
                                        aaCGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGACCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGc                                                                                                               <  2:1250301/150‑1 (MQ=255)
                                                ccATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGaa                                                                                                       >  1:378917/1‑150 (MQ=255)
                                                     tgGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCt                                                                                                  >  1:1003512/1‑150 (MQ=255)
                                                       gTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTc                                                                                                <  2:64679/150‑1 (MQ=255)
                                                          aCGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAt                                                                                             >  1:1260918/1‑150 (MQ=255)
                                                           cGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtg                                                                                            >  1:207498/1‑150 (MQ=255)
                                                                agagCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATc                                                                                       <  2:1035986/150‑1 (MQ=255)
                                                                   gCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                    >  1:549930/1‑150 (MQ=255)
                                                                   gCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                    >  1:555581/1‑150 (MQ=255)
                                                                   gCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                    >  1:1326736/1‑150 (MQ=255)
                                                                    cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                                                                   <  2:370162/150‑1 (MQ=255)
                                                                    cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                                                                   <  1:1149433/150‑1 (MQ=255)
                                                                      aaaCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTc                                                                                 <  2:469089/150‑1 (MQ=255)
                                                                           ctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                            >  2:730435/1‑150 (MQ=255)
                                                                           ctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                            <  2:1634046/150‑1 (MQ=255)
                                                                              tCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTaa                                                                         >  1:542014/1‑150 (MQ=255)
                                                                                 gCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCa                                                                      >  2:1375406/1‑150 (MQ=255)
                                                                                   aTCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTTATCAcc                                                                    >  2:630924/1‑150 (MQ=255)
                                                                                    tCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCg                                                                   <  2:273378/150‑1 (MQ=255)
                                                                                         tCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCATCAATGTGATCTTTGTCTTCGCTAATCACCGGGaaa                                                              >  2:1243710/1‑150 (MQ=255)
                                                                                           tACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACg                                                            <  2:805508/150‑1 (MQ=255)
                                                                                             cgcgCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg                                                          <  1:1146217/150‑1 (MQ=255)
                                                                                             cgcgCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg                                                          <  2:776298/150‑1 (MQ=255)
                                                                                             cgcgCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg                                                          <  1:738504/150‑1 (MQ=255)
                                                                                             cgcgCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg                                                          <  2:515536/150‑1 (MQ=255)
                                                                                             cgcgCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg                                                          <  2:1483156/150‑1 (MQ=255)
                                                                                             cgcgCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg                                                          >  1:572508/1‑150 (MQ=255)
                                                                                              gcgcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGa                                                         >  2:1600819/1‑150 (MQ=255)
                                                                                               cgcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAg                                                        <  1:759699/150‑1 (MQ=255)
                                                                                               cgcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAg                                                        >  1:485387/1‑150 (MQ=255)
                                                                                                 cTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTg                                                      <  1:294795/150‑1 (MQ=255)
                                                                                                    aCTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCTGGAAACGTGAGTGGGc                                                   <  1:431295/150‑1 (MQ=255)
                                                                                                    aCTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGc                                                   >  2:1006830/1‑150 (MQ=255)
                                                                                                     cTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCg                                                  >  1:1633153/1‑150 (MQ=255)
                                                                                                       ttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGa                                                >  2:450184/1‑150 (MQ=255)
                                                                                                        tCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTTATCACCGGGAAACGTGAGTGGGCGGAc                                               <  1:630924/150‑1 (MQ=255)
                                                                                                         cAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACt                                              <  2:378917/150‑1 (MQ=255)
                                                                                                           ggAACGACAACCGCCTGACGTAACACTTTGTCCGTGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCg                                            >  1:839224/1‑150 (MQ=255)
                                                                                                              aCGACAACGGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTGAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatt                                         <  2:587962/150‑2 (MQ=255)
                                                                                                                 aCAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGatgatg                                      <  2:711045/150‑1 (MQ=255)
                                                                                                                 aCAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCAAATCACCGGGAAACGTGAGTGGGCGGACTCGatgatg                                      >  2:748760/1‑150 (MQ=255)
                                                                                                                    aCCGCCTGTCGTTACACTTTGTCCATGCTGAATGCTTCAGCATCGCTGCGCATAAAAGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGTTCTTTGTCTTCGCTAATCACGGGGAAACGTGAGTGGGCTGACTCGATGATGACa                                   <  2:1062577/150‑1 (MQ=255)
                                                                                                                    aCCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACa                                   <  1:22297/150‑1 (MQ=255)
                                                                                                                    aCCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACa                                   <  1:450184/150‑1 (MQ=255)
                                                                                                                      cGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACATGGAAACGTGAGTGGGCGGACTCGATGATGACATc                                 <  1:356416/150‑1 (MQ=255)
                                                                                                                       gCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTggg                                                    <  2:1611440/130‑1 (MQ=255)
                                                                                                                       gCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTggg                                                    >  1:1611440/1‑130 (MQ=255)
                                                                                                                       gCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCa                                >  1:530384/1‑150 (MQ=255)
                                                                                                                          tGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATTATGACATCAAGa                             <  2:1475425/150‑1 (MQ=255)
                                                                                                                             cGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACAt                          <  1:1334177/150‑1 (MQ=255)
                                                                                                                                   acTTTGTCCATGCTGAAGGCTTCATCATCGCTTCTCATAAACGGCAGCAAGTCTTTCGCCATCATAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTcc                    <  1:527619/150‑1 (MQ=255)
                                                                                                                                       tGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAGTCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCg                <  1:374060/150‑1 (MQ=255)
                                                                                                                                        gTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGt               >  1:615996/1‑150 (MQ=255)
                                                                                                                                         gcgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTc              >  1:179846/4‑150 (MQ=255)
                                                                                                                                          cgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCt             >  1:670736/3‑150 (MQ=255)
                                                                                                                                          cgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCt             >  1:832976/3‑150 (MQ=255)
                                                                                                                                             tGCTGAAGGCTTCAGAATCGCTGCGCATAAACGCCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGt          >  2:272796/1‑150 (MQ=255)
                                                                                                                                               cTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTg        >  1:1579590/1‑150 (MQ=255)
                                                                                                                                                     gCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCCTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTc  <  1:1243710/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTC  >  NC_000913/691241‑691539

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