Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,698 G→T 9.7% intergenic (+244/+93) codA → / ← cynR cytosine/isoguanine deaminase/DNA‑binding transcriptional dual regulator CynR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,6980GT 9.7% 325.4 / 17.6 108intergenic (+244/+93)codA/cynRcytosine/isoguanine deaminase/DNA‑binding transcriptional dual regulator CynR
Reads supporting (aligned to +/- strand):  ref base G (48/47);  new base T (9/2);  total (59/49)
Fisher's exact test for biased strand distribution p-value = 5.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

Only 100 of 111 total aligned reads displayed.

GTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA  >  NC_000913/357549‑357847
                                                                                                                                                     |                                                                                                                                                      
gTAGCTCATCAGGCAATTCCATCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAAACGCGCATTCCTCCAACGCCATGTTCAAAAATGCCTTCTCTGCGGCTGTCTGCCATCTGTAGTTTATGCCGGATGCGGCGTg                                                                                                                                                        <  2:1977272/150‑1 (MQ=18)
gTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                                                                                                                                                        <  1:2313775/150‑1 (MQ=255)
 tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCTGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCTCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGa                                                                                                                                                       <  2:44288/150‑1 (MQ=255)
  aGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGaa                                                                                                                                                      <  2:79111/150‑1 (MQ=255)
   gcgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCATCAACCGCGCATTCCTCCAACTCCATGTGATAAAATGCCTTCGCAGCGGCGGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAAc                                                                                                                                                     <  2:2200226/150‑1 (MQ=18)
    cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACg                                                                                                                                                    <  1:690059/150‑1 (MQ=255)
     gcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                                   >  1:1932942/1‑150 (MQ=255)
     gcATCAGGCAATTCCAGCCGCTGATCTGTGGCAGCGGCTACCTTGATTCATTCCCGCCACCAACCGCGCATTCCTCCCACGCCATATGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                                   <  1:1406400/150‑1 (MQ=255)
      cATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                                                  <  1:814668/150‑1 (MQ=255)
            gCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                            >  1:1969041/1‑150 (MQ=255)
             cAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                           >  2:524048/1‑150 (MQ=255)
              aaTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          <  1:1662242/150‑1 (MQ=255)
              aaTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          <  1:2014996/150‑1 (MQ=255)
               aTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAAAAACCAAGAAATACTAGAACGCCAAGAGAAAAAATACCTTAACATCTGCTTTCTCCCAGATTTAGTTTATAAAAAATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                         >  2:1135854/1‑150 (MQ=12)
                ttCCAGCCGCTTATCTGTGTCAGCGGATACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCAAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                                                        <  2:1051561/150‑2 (MQ=18)
                ttCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGTCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                        <  1:1515580/150‑1 (MQ=255)
                   cAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTATTTAATCCCCGAGGAGGCGTGACAGCCTTAGcctgggtgc                                                                                                                                     >  1:75593/1‑143 (MQ=17)
                   cAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                     >  1:449988/1‑150 (MQ=255)
                    aGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                    <  2:1685165/150‑1 (MQ=255)
                     gCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                                                                                                                   >  2:1412376/1‑150 (MQ=255)
                      ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                                                                                                                  <  2:1672260/150‑1 (MQ=255)
                          tGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTTCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGc                                                                                                                              <  2:476424/150‑1 (MQ=255)
                            aTCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCAc                                                                                                                            <  1:56588/150‑1 (MQ=255)
                              cTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTg                                                                                                                          >  1:2334759/1‑150 (MQ=255)
                               tgtgTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGa                                                                                                                         >  1:211068/1‑150 (MQ=255)
                                 tgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                                       <  2:810438/150‑1 (MQ=255)
                                 tgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                                       <  2:609964/150‑1 (MQ=255)
                                     aGCTGCTAACTTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCCCAGCTGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGt                                                                                                                  <  1:1639340/150‑1 (MQ=37)
                                                  attcattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGc                                                                                                     >  1:2016489/1‑150 (MQ=255)
                                                  attcattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTTCAGTGCTCTACGTAGGCCGGATAAgg                                                                                                     <  1:330237/150‑2 (MQ=21)
                                                      attcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTACAAAAATGCCTTCGCAGCGGCTTTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCCGATAAg                                                                                                      <  2:2187739/145‑1 (MQ=255)
                                                      attcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTACAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCCGATAAg                                                                                                      >  1:2187739/1‑145 (MQ=255)
                                                        tcCCGCCACCAACCGCGCATTCCTCCCACGCCATGTGCAAAAATGCCTTCGCAGCGGCGGTCTTCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcg                                                                                               <  1:1250388/150‑1 (MQ=255)
                                                         cccGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGAAGCGGCGTCCGGTCCATAACCGGCGTCCGGAGAGAGCTGACCGAGTAGCCCTCATAAGCGGAccgc                                                                                              >  1:42078/1‑150 (MQ=12)
                                                          actacaataaCCGTGGATTATTCTAACGCCATGGTCTAAAATGCATTCTCTGCGGCTGTATTCCTTCTGTAGTTTATGCCGGATGCGGAGTTAACGCGTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCa                                                                                             <  2:1317955/146‑1 (MQ=9)
                                                             ccAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCTCAGCTGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCa                                                                                          <  2:19357/150‑1 (MQ=255)
                                                                 caacCGCGCATTCCTCACACGCCATGTGCAAAAATGGCTTCTCTTCGGCTTTATTGCTTCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAGTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                      <  2:2222256/150‑1 (MQ=17)
                                                                      gcgcATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTcc                                                                                 <  2:1932942/150‑1 (MQ=255)
                                                                           ttCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCTGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCg                                                                            <  2:200099/150‑1 (MQ=255)
                                                                           ttCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATACGGCCTACGTAGCGCACTGAACTCGTAGGCCTGATAAGAGTAGCGGATCAGGCAATTCCAgaag                                                                            >  1:284498/1‑147 (MQ=21)
                                                                                 cAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGAcc                                                                      <  1:524048/150‑1 (MQ=255)
                                                                                     gCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATCCCGGATGCGGCGTGAAATCATTATCCGGCATACGTAGAGCACTGAAATCGTAGGCCTGATAAGAGTAGAGAATCAGGAAATTCCAGCCGAAGACAtgtg                                                                  >  2:1426/1‑150 (MQ=37)
                                                                                      ccATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGACTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTgtgt                                                                 >  1:1232660/1‑150 (MQ=255)
                                                                                           tgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCg                                                            >  1:78571/1‑150 (MQ=255)
                                                                                           ggCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCg                                                            <  2:1969041/149‑1 (MQ=255)
                                                                                             cAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACTCCTTATACGGGATACGTAGAGAACTGAACTCGTAGGCCTGATAAGCGTAGCGCGGCAGGCGATTCACGCCGCAGACGTGTGTGAGCggg                                                          >  2:1584422/1‑149 (MQ=18)
                                                                                             cAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGc                                                          <  1:2215987/150‑1 (MQ=255)
                                                                                                aaaTGCCTTCGCAGCGGCTGTCTGCCATCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTAc                                                       <  1:2411146/150‑1 (MQ=255)
                                                                                                    gCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                   <  1:1403639/150‑1 (MQ=255)
                                                                                                    gCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                   <  1:1412376/150‑1 (MQ=255)
                                                                                                    gCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                   <  1:398291/150‑1 (MQ=255)
                                                                                                    gCCTTCGCAGCGGCTGTCTGACAGCTGTAGATTATGCAGGATGCGGAGTGAACGAATTATACGGCCTAAGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                   >  2:305965/1‑150 (MQ=17)
                                                                                                     ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGa                                                  <  1:810781/150‑1 (MQ=255)
                                                                                                      cTTCGCAGCGGCTGTCTGCCAGCTTTAGTTTATGCCGGAAGAGGCGTGAAAGCATTATCCTGCATACGTAGAGCACTGAAATCGTAGGCCTGATAAGCGTAGCGCAGCAGGCAATTCAAGCAGAAGACCTGTGTAAGAGGCTACCGTGAt                                                 >  2:1325951/1‑150 (MQ=17)
                                                                                                      cTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAt                                                 >  1:1393835/1‑150 (MQ=255)
                                                                                                      cTTCGCAGCGCCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAt                                                 >  1:2077756/1‑150 (MQ=255)
                                                                                                         cGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTTAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCa                                              <  2:397151/150‑1 (MQ=255)
                                                                                                          gCAGCGGCTGTCTGCCAGCTGTTTTTAATCCCGGAAACAGCGTCCAATCATTAAACGTCCTACGTAGAGCAATGAACTAGTAGGACTGATAAGCGTAGCGAATCAGGCAATTCCAGCCGCATACATATGTAAGCGGCTACCGTGATTCAt                                             >  2:1696689/1‑150 (MQ=255)
                                                                                                              cGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCCGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTcc                                         >  1:1585785/1‑150 (MQ=255)
                                                                                                               ggCTGTCTTCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCg                                        <  2:2044331/150‑1 (MQ=255)
                                                                                                                gCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGc                                       <  2:2016489/150‑1 (MQ=255)
                                                                                                                 cggTCTGCCAGCTGTAGTTTATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCATGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGcc                                      <  1:1135854/148‑1 (MQ=255)
                                                                                                                 cTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGcc                                      >  2:1609347/1‑150 (MQ=255)
                                                                                                                       gCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCaacaac                                >  1:1790795/1‑150 (MQ=255)
                                                                                                                       gCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCaacaac                                >  1:1154578/1‑150 (MQ=255)
                                                                                                                       gCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCaacaac                                >  2:2451690/1‑150 (MQ=255)
                                                                                                                           gCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGt                                                    <  2:543314/126‑1 (MQ=39)
                                                                                                                           gCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcgcg                            >  2:77738/1‑150 (MQ=40)
                                                                                                                            ccgTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTcc                                                                                 <  1:1764867/94‑1 (MQ=14)
                                                                                                                             tGTAGTTTATGACGGATTCGGCGTGAAAGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCa                          >  2:1619909/1‑150 (MQ=255)
                                                                                                                                aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAttt                       >  1:2215098/1‑150 (MQ=38)
                                                                                                                                 gTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcgcg                            <  1:1123203/144‑1 (MQ=37)
                                                                                                                                 gTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcgcg                            >  2:1123203/1‑144 (MQ=37)
                                                                                                                                 gTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTa                      >  2:1029475/1‑150 (MQ=38)
                                                                                                                                   ttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATc                    <  1:1149545/150‑1 (MQ=40)
                                                                                                                                    tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATcc                   <  2:1393835/150‑1 (MQ=40)
                                                                                                                                     aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCa                  >  1:1124842/3‑150 (MQ=14)
                                                                                                                                     aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCa                  >  1:2231815/1‑150 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCaa                 >  2:1394841/1‑150 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCaa                 >  1:1221105/1‑150 (MQ=255)
                                                                                                                                      ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCaa                 >  2:1586789/2‑150 (MQ=14)
                                                                                                                                       gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCCTTTATCCAAc                >  2:2449680/1‑150 (MQ=14)
                                                                                                                                       gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAAc                >  1:162516/1‑150 (MQ=14)
                                                                                                                                       gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAAc                >  1:2295117/1‑150 (MQ=14)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg               >  1:1252539/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg               <  1:355096/150‑1 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg               >  1:226622/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATAAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg               >  2:2176827/1‑150 (MQ=255)
                                                                                                                                        ccGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACg               >  1:615885/1‑150 (MQ=14)
                                                                                                                                           gataaggcgTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCa            >  1:23113/6‑150 (MQ=14)
                                                                                                                                           gataaggcgTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCa            >  1:851773/6‑150 (MQ=14)
                                                                                                                                             tGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTCCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtg          >  2:821292/1‑150 (MQ=255)
                                                                                                                                               cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGCTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtgtg        >  1:864709/1‑150 (MQ=255)
                                                                                                                                               cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtgtg        >  1:397534/1‑150 (MQ=255)
                                                                                                                                               cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtgtg        <  2:1760636/150‑1 (MQ=255)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                <  1:638420/109‑1 (MQ=35)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                >  2:638420/1‑109 (MQ=35)
                                                                                                                                                ggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGc       >  2:1858390/1‑150 (MQ=255)
                                                                                                                                                  cgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGCCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaa     >  1:1897529/1‑150 (MQ=255)
                                                                                                                                                     gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaaa  >  1:1006156/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA  >  NC_000913/357549‑357847

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