Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,351 T→A 5.5% S64C (AGT→TGT)  alkA ← 3‑methyl‑adenine DNA glycosylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,3510TA 5.5% 261.6 / 3.3 87S64C (AGT→TGT) alkA3‑methyl‑adenine DNA glycosylase II
Reads supporting (aligned to +/- strand):  ref base T (44/38);  new base A (2/3);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

AACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAAC  >  NC_000913/2147208‑2147496
                                                                                                                                                 |                                                                                                                                                 
aaCACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAAt                                                                                                                                              <  2:33414/149‑1 (MQ=255)
aaCACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGAAGATCAAACAGGCGGCTCATTTTCGCAAGAAACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAAt                                                                                                                                              >  2:1491180/1‑149 (MQ=255)
 acacAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAAtt                                                                                                                                             >  1:1347989/1‑149 (MQ=255)
   acaGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTa                                                                                                                                           <  2:1164647/149‑1 (MQ=255)
    caGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTAt                                                                                                                                          <  1:1326846/149‑1 (MQ=255)
         gggTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCACCGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCa                                                                                                                                     >  2:412573/1‑149 (MQ=255)
          ggTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACATAAATTTAAGTGCag                                                                                                                                    <  2:313134/149‑1 (MQ=255)
              aaCGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTa                                                                                                                                >  1:1034927/1‑149 (MQ=255)
                 gCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATgg                                                                                                                             <  1:1021951/149‑1 (MQ=255)
                     tCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCCACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCggg                                                                                                                         >  2:1041868/1‑149 (MQ=255)
                     tCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCggg                                                                                                                         >  2:186888/1‑149 (MQ=255)
                     tCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCggg                                                                                                                         <  2:323868/149‑1 (MQ=255)
                       cGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCt                                                                                                                       >  2:544579/1‑149 (MQ=255)
                        gggCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGGTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGGTCTAAACCGGCACTTAAATTAATCTGCTGACTATGGCGGGCta                                                                                                                      >  1:426796/1‑149 (MQ=255)
                          gccgcgccgcGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGATACATAGCAGATCCAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAGCCTGCACTGAAATTTATGTGCAGAGTATGGCTGGCtata                                                                                                                    >  1:82458/1‑149 (MQ=255)
                           ccgcgccgcgCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTatat                                                                                                                   <  1:465679/149‑1 (MQ=255)
                           ccgcgccgcgCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCCGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTatat                                                                                                                   <  2:105752/149‑1 (MQ=255)
                                ccgcgccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGCCTCTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGa                                                                                                              >  1:819623/1‑149 (MQ=255)
                                 ggggccTACCGTGCCCACGGACCCGTTAACATTTTGTGGGTTCCATTGCTGTTAAACAGGGGGGTTATTTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGaa                                                                                                             <  1:1492968/146‑1 (MQ=255)
                                     ccTAACCTGCCCAACGCACCGATAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGGTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGc                                                                                                         <  1:692600/149‑1 (MQ=255)
                                      cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGaa                                                                                                             >  1:398794/1‑144 (MQ=255)
                                      cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGaa                                                                                                             <  2:398794/144‑1 (MQ=255)
                                      cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCa                                                                                                        <  2:1419911/149‑1 (MQ=255)
                                       tAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAg                                                                                                       <  2:373524/149‑1 (MQ=255)
                                            tGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTcacc                                                                                                  >  1:263604/1‑149 (MQ=255)
                                               ccAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTGTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCAcg                                                                                               >  2:11251/1‑149 (MQ=255)
                                                 aaCGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCAcgcc                                                                                             >  1:472347/1‑149 (MQ=255)
                                                    gCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCcgcg                                                                                          >  2:337437/1‑149 (MQ=255)
                                                     cACCGTTAACAATTTGTGGGTTACATTGCAGATCTAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGATTAGCAGTCACAACGCCGCGa                                                                                         >  1:319612/1‑149 (MQ=255)
                                                     cACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGa                                                                                         <  2:25752/149‑1 (MQ=255)
                                                        cGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGatat                                                                                      <  2:437299/149‑1 (MQ=255)
                                                        cGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGatat                                                                                      <  1:959234/149‑1 (MQ=255)
                                                           tAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTcc                                                                                   >  2:860020/1‑148 (MQ=255)
                                                           tAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCg                                                                                   >  1:717694/1‑149 (MQ=255)
                                                            aaCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGc                                                                                  >  2:708717/1‑149 (MQ=255)
                                                             aCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGcc                                                                                 >  2:514386/1‑149 (MQ=255)
                                                               aaTTTGTGGGTTACATTGCAGATCAAACATGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCa                                                                               >  1:89274/1‑149 (MQ=255)
                                                               aaTTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCa                                                                               >  2:475857/1‑149 (MQ=255)
                                                                aTTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCAc                                                                              <  1:514386/149‑1 (MQ=255)
                                                                  ttGTGGGTTACATTCATATTCAAACCGCCGCCCAATTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGc                                                                           <  1:1218187/149‑1 (MQ=255)
                                                                     tgGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCa                                                                         >  1:1437644/1‑149 (MQ=255)
                                                                       ggTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAgc                                                                       >  2:1502093/1‑148 (MQ=255)
                                                                            cATTGCAGATCATACATGAGGCTCATTTTCGCCAGACACTCTGCGGCAAGAGGTTCTAAACCTGCACATAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACg                                                                  <  1:234983/149‑1 (MQ=255)
                                                                             aTTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACgg                                                                 <  1:412573/149‑1 (MQ=255)
                                                                             aTTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACgg                                                                 <  1:171071/149‑1 (MQ=255)
                                                                                gCAGATCATACAGGCGTCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCAGAGATATTCGCCCACCGCCAGACTACGGGCa                                                              <  2:1275874/149‑1 (MQ=255)
                                                                                 cAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCat                                                             >  2:34442/1‑149 (MQ=255)
                                                                                  aGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCata                                                            >  2:1052636/1‑149 (MQ=255)
                                                                                   gATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCataa                                                           >  1:1022883/1‑149 (MQ=255)
                                                                                       cAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCAtaataa                                                        <  2:1347989/149‑1 (MQ=255)
                                                                                        aaaCAGGCGGCTCATTTTCTCCAGACACTCAGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAAc                                                       <  1:1406532/149‑1 (MQ=255)
                                                                                        aaaCAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAAc                                                       >  1:434269/1‑149 (MQ=255)
                                                                                         aaCAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACt                                                      <  2:263604/149‑1 (MQ=255)
                                                                                         aaCAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACt                                                      >  1:1007327/1‑149 (MQ=255)
                                                                                                ggcTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCcgcg                                               >  2:261757/1‑149 (MQ=255)
                                                                                                  cTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAc                                             >  2:15767/1‑149 (MQ=255)
                                                                                                   tCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAcc                                            <  2:347082/149‑1 (MQ=255)
                                                                                                       ttttCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGtt                                         >  1:805414/1‑149 (MQ=255)
                                                                                                       ttttCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGtt                                         >  1:953451/1‑149 (MQ=255)
                                                                                                        tttCTCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGttt                                        <  1:1491180/149‑1 (MQ=255)
                                                                                                        tttCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGttt                                        >  1:845166/1‑149 (MQ=255)
                                                                                                         ttCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTc                                       <  1:734894/149‑1 (MQ=255)
                                                                                                         ttCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTc                                       <  2:389746/149‑1 (MQ=255)
                                                                                                           cGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTccc                                     >  1:1108319/1‑148 (MQ=255)
                                                                                                           cGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCca                                     <  1:882111/149‑1 (MQ=255)
                                                                                                           cGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCca                                     <  2:446605/149‑1 (MQ=255)
                                                                                                             ccAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCcaca                                   <  1:34442/149‑1 (MQ=255)
                                                                                                                   acTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGGCCGTTTCCTCACTGCTc                             >  1:979274/1‑149 (MQ=255)
                                                                                                                     tctGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCAc                           >  2:1356610/1‑149 (MQ=255)
                                                                                                                         cGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCa                       >  1:1261903/1‑149 (MQ=255)
                                                                                                                           gCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTAACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATACTAACAGTCCGCGACCGTGTCAACACTGCTCACCGCACg                     >  2:566360/1‑149 (MQ=255)
                                                                                                                            cAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCAACGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACgg                    >  1:880256/1‑149 (MQ=255)
                                                                                                                             aaCAGGTTCTTAACCTTCACTTTAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGGGGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACTGTTTCCACACTTCTCACCGCACggg                   <  2:1217363/149‑1 (MQ=255)
                                                                                                                             aaCAGGTTCTAAACCTGAAATTTAATTTATGTGCAGTGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACggg                   <  2:82458/149‑1 (MQ=255)
                                                                                                                                 ggTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggcggc               <  2:1034927/149‑1 (MQ=255)
                                                                                                                                 ggTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGggcggc               <  2:805414/149‑1 (MQ=255)
                                                                                                                                  gTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGgcggcg              >  1:47203/1‑149 (MQ=255)
                                                                                                                                   ttCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCga             >  2:1122819/1‑149 (MQ=255)
                                                                                                                                    tCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCgag            >  1:1056052/1‑149 (MQ=255)
                                                                                                                                      atactcTGCACATAAATTTAAGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGaa          >  1:1264057/6‑149 (MQ=255)
                                                                                                                                       tactcTGCACATAAATTTAAGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGaaa         >  1:746/5‑149 (MQ=255)
                                                                                                                                          tcTGCACATAAATTTAAGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATcc      <  2:1264057/148‑1 (MQ=255)
                                                                                                                                           cTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATccc     <  2:1056052/149‑1 (MQ=255)
                                                                                                                                           cTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATccc     >  2:1234534/1‑149 (MQ=255)
                                                                                                                                           cTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATccc     <  2:1261903/149‑1 (MQ=255)
                                                                                                                                              cACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAAc  <  1:1041868/149‑1 (MQ=255)
                                                                                                                                              cACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAAc  <  1:186888/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAAC  >  NC_000913/2147208‑2147496

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: