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

AACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAAC  >  NC_000913/2147208‑2147496
                                                                                                                                                 |                                                                                                                                                 
aaCACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAAt                                                                                                                                              <  2:33414/149‑1 (MQ=255)
aaCACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGAAGA‑TCAAACAGGCGGCTCA‑TTTTCGCAAGAAACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAAt                                                                                                                                              >  2:1491180/1‑149 (MQ=255)
 acacAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAAtt                                                                                                                                             >  1:1347989/1‑149 (MQ=255)
   acaGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTa                                                                                                                                           <  2:1164647/149‑1 (MQ=255)
    caGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTAt                                                                                                                                          <  1:1326846/149‑1 (MQ=255)
         gggTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCACCGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCa                                                                                                                                     >  2:412573/1‑149 (MQ=255)
          ggTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACATAAATTTAAGTGCag                                                                                                                                    <  2:313134/149‑1 (MQ=255)
              aaCGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTa                                                                                                                                >  1:1034927/1‑149 (MQ=255)
                 gCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATgg                                                                                                                             <  1:1021951/149‑1 (MQ=255)
                     tCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCCACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCggg                                                                                                                         >  2:1041868/1‑149 (MQ=255)
                     tCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCggg                                                                                                                         >  2:186888/1‑149 (MQ=255)
                     tCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCggg                                                                                                                         <  2:323868/149‑1 (MQ=255)
                       cGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCt                                                                                                                       >  2:544579/1‑149 (MQ=255)
                        gggCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGGTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGGTCTAAACCGGCACTTAAATTAATCTGCTGACTATGGCGGGCta                                                                                                                      >  1:426796/1‑149 (MQ=255)
                          gccgcgccgcGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGATACATAGCAGA‑TCCAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAGCCTGCACTGAAATTTATGTGCAGAGTATGGCTGGCtata                                                                                                                    >  1:82458/1‑149 (MQ=255)
                           ccgcgccgcgCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTatat                                                                                                                   <  1:465679/149‑1 (MQ=255)
                           ccgcgccgcgCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCCGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTatat                                                                                                                   <  2:105752/149‑1 (MQ=255)
                                ccgcgccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGCCTCTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGa                                                                                                              >  1:819623/1‑149 (MQ=255)
                                 ggggccTACCGTGCCCACGGACCCGTTAACATTTTGTGGGTTCCATTGCTG‑‑TTAAACAGGGGGGTTATTTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGaa                                                                                                             <  1:1492968/146‑1 (MQ=255)
                                     ccTAACCTGCCCAACGCACCGATAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGGTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGc                                                                                                         <  1:692600/149‑1 (MQ=255)
                                      cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGaa                                                                                                             >  1:398794/1‑144 (MQ=255)
                                      cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGaa                                                                                                             <  2:398794/144‑1 (MQ=255)
                                      cTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCa                                                                                                        <  2:1419911/149‑1 (MQ=255)
                                       tAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAg                                                                                                       <  2:373524/149‑1 (MQ=255)
                                            tGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTcacc                                                                                                  >  1:263604/1‑149 (MQ=255)
                                               ccAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTGTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCAcg                                                                                               >  2:11251/1‑149 (MQ=255)
                                                 aaCGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCAcgcc                                                                                             >  1:472347/1‑149 (MQ=255)
                                                    gCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCcgcg                                                                                          >  2:337437/1‑149 (MQ=255)
                                                     cACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCTAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGATTAGCAGTCACAACGCCGCGa                                                                                         >  1:319612/1‑149 (MQ=255)
                                                     cACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGa                                                                                         <  2:25752/149‑1 (MQ=255)
                                                        cGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGatat                                                                                      <  2:437299/149‑1 (MQ=255)
                                                        cGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGatat                                                                                      <  1:959234/149‑1 (MQ=255)
                                                           tAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTcc                                                                                   >  2:860020/1‑148 (MQ=255)
                                                           tAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCg                                                                                   >  1:717694/1‑149 (MQ=255)
                                                            aaCAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGc                                                                                  >  2:708717/1‑149 (MQ=255)
                                                             aCAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGcc                                                                                 >  2:514386/1‑149 (MQ=255)
                                                               aaTTTGTGGGTTACATTGCAGA‑TCAAACATGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCa                                                                               >  1:89274/1‑149 (MQ=255)
                                                               aaTTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCa                                                                               >  2:475857/1‑149 (MQ=255)
                                                                aTTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCAc                                                                              <  1:514386/149‑1 (MQ=255)
                                                                  ttGTGGGTTACATT‑CATATTCAAAC‑CGCCGCCCA‑ATTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGc                                                                           <  1:1218187/149‑1 (MQ=255)
                                                                     tgGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCa                                                                         >  1:1437644/1‑149 (MQ=255)
                                                                       ggTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAgc                                                                       >  2:1502093/1‑148 (MQ=255)
                                                                            cATTGCAGA‑TCATACATGAGGCTCA‑TTTTCGCCAGACACTCTGCGGCAAGAGGTTCTAAACCTGCACATAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACg                                                                  <  1:234983/149‑1 (MQ=255)
                                                                             aTTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACgg                                                                 <  1:412573/149‑1 (MQ=255)
                                                                             aTTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACgg                                                                 <  1:171071/149‑1 (MQ=255)
                                                                                gCAGA‑TCATACAGGCGTCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCAGAGATATTCGCCCACCGCCAGACTACGGGCa                                                              <  2:1275874/149‑1 (MQ=255)
                                                                                 cAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCat                                                             >  2:34442/1‑149 (MQ=255)
                                                                                  aGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCata                                                            >  2:1052636/1‑149 (MQ=255)
                                                                                   gA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCataa                                                           >  1:1022883/1‑149 (MQ=255)
                                                                                       cAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCAtaataa                                                        <  2:1347989/149‑1 (MQ=255)
                                                                                        aaaCAGGCGGCTCA‑TTTTCTCCAGACACTCAGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAAc                                                       <  1:1406532/149‑1 (MQ=255)
                                                                                        aaaCAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAAc                                                       >  1:434269/1‑149 (MQ=255)
                                                                                         aaCAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACt                                                      <  2:263604/149‑1 (MQ=255)
                                                                                         aaCAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACt                                                      >  1:1007327/1‑149 (MQ=255)
                                                                                                ggcTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCcgcg                                               >  2:261757/1‑149 (MQ=255)
                                                                                                  cTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAc                                             >  2:15767/1‑149 (MQ=255)
                                                                                                   tCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAcc                                            <  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)
                                                                                                                                                 |                                                                                                                                                 
AACACAGCCGGGTAAACGCAATCCGGGCCGCGCCGCGCCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGA‑TCAAACAGGCGGCTCA‑TTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAAC  >  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: ‑