Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,902,194 G→T 7.6% intergenic (‑56/+88) PP_mr24 ← / ← rdgC mini‑ykkC regulatory region/recombination‑associated protein RdgC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,902,1940GT 7.6% 821.6 / 42.0 249intergenic (‑56/+88)PP_mr24/rdgCmini‑ykkC regulatory region/recombination‑associated protein RdgC
Reads supporting (aligned to +/- strand):  ref base G (122/108);  new base T (7/12);  total (129/120)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 263 total aligned reads displayed.

AACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAATTCG  >  NC_002947/1902045‑1902342
                                                                                                                                                        |                                                                                                                                                    
aaCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAgg                                                                                                                                                      <  1:1777259/150‑1 (MQ=255)
 aCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                     <  2:1444699/150‑1 (MQ=255)
 aCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                     <  1:3225421/150‑1 (MQ=255)
 aCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                     <  1:6129617/150‑1 (MQ=255)
 aCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                     <  2:5564421/150‑1 (MQ=255)
 aCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                     <  2:2887629/150‑1 (MQ=255)
 aCAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                     <  2:3907332/150‑1 (MQ=255)
  cAGGATGATCCAGGACATGATTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAgggg                                                                                                                                                    <  1:3967932/150‑1 (MQ=255)
    ggATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGtt                                                                                                                                                  <  1:6365438/150‑1 (MQ=255)
     gatgatCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTa                                                                                                                                                 >  1:5662023/1‑150 (MQ=255)
          tCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACt                                                                                                                                            >  2:4791702/1‑150 (MQ=255)
          tCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACt                                                                                                                                            >  1:198357/1‑150 (MQ=255)
           ccAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACtt                                                                                                                                           >  1:5028075/1‑150 (MQ=255)
           ccAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACtt                                                                                                                                           <  1:501463/150‑1 (MQ=255)
           ccAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACtt                                                                                                                                           <  2:5783582/150‑1 (MQ=255)
             aGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCCAGACCGCGCCTGCACGGGGTTACCACTGGc                                                                                                                                         >  2:6080135/1‑150 (MQ=255)
             aGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGTGGTAACCCCTTGc                                                                                                                                         <  2:84530/150‑1 (MQ=255)
             aGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGc                                                                                                                                         >  2:1744210/1‑150 (MQ=255)
              ggACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCa                                                                                                                                        <  1:2641864/150‑1 (MQ=255)
               gACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTCTTCGCGGGTAAACCCGCGCCTGCAAGTGGTAACCCCTTg                                                                                                                                          <  1:6756822/150‑1 (MQ=37)
               gACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAggg                                                                                                                                                     >  1:1752786/1‑136 (MQ=255)
               gACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGTGGTAACCCCTTGCAg                                                                                                                                       <  1:633/150‑1 (MQ=255)
                  aTGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGgcg                                                                                                                                    >  2:1761230/1‑150 (MQ=255)
                   tGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTCGCGGGTAAACCCGCGCCTGCAAGTGGTAACCCCTTGCAGgcgc                                                                                                                                   <  2:5712574/150‑1 (MQ=38)
                     aGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCgg                                                                                                                                 <  1:4468718/150‑1 (MQ=255)
                                  cGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGaa                                                                                                                    <  2:4877768/150‑1 (MQ=255)
                                    tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCCCCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCAGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGcgacga                                                                                                                  >  1:4511154/1‑150 (MQ=255)
                                    tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAATGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAaga                                                                                                                  >  1:3709511/1‑150 (MQ=255)
                                    tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAaga                                                                                                                  >  2:3125027/1‑150 (MQ=255)
                                    tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAaga                                                                                                                  <  1:4298103/150‑1 (MQ=255)
                                      tGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAgg                                                                                                                >  1:1047308/1‑150 (MQ=255)
                                       gAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGc                                                                                                               <  2:2802305/150‑1 (MQ=255)
                                        aaTGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGcc                                                                                                              >  2:6297480/1‑150 (MQ=255)
                                             tgGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGc                                                                                                         >  1:2828851/1‑150 (MQ=255)
                                             tgGGGCCGGCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGc                                                                                                         <  2:3577352/150‑1 (MQ=255)
                                               gggCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTg                                                                                                       >  1:597370/1‑150 (MQ=255)
                                                ggCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTgg                                                                                                      <  2:1020481/150‑1 (MQ=255)
                                                       cccGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGaaaa                                                                                               <  1:5131844/150‑1 (MQ=255)
                                                       cccGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGaaaa                                                                                               >  2:3210056/1‑150 (MQ=255)
                                                         cGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACg                                                                                             >  2:6374565/1‑150 (MQ=255)
                                                          ggAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGa                                                                                            <  1:1848014/150‑1 (MQ=255)
                                                                 tCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCAGGGTGCCTGTCTGCGGC‑‑CTGGCGGGCAAGACCGCGCCTGCACGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGc                                                                                     <  1:161496/150‑1 (MQ=255)
                                                                   aCGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑ATTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCa                                                                                   <  2:5447137/150‑1 (MQ=255)
                                                                     gcgGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGc                                                                                 >  2:6129729/1‑150 (MQ=255)
                                                                       ggTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGTAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCAc                                                                               <  2:3795788/150‑1 (MQ=255)
                                                                       ggTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGTGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCAc                                                                               >  2:1114199/1‑150 (MQ=255)
                                                                       ggTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGAGTCAGGCCAGCAc                                                                               >  2:4558899/1‑150 (MQ=255)
                                                                           cgtcgtGAGGTCGTCCTCCCCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGACAACGATTCAGGCCAGCACCTGc                                                                           <  1:149302/150‑1 (MQ=255)
                                                                            acggtGAGGACGACCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCa                                                                          >  1:2540441/4‑150 (MQ=255)
                                                                              cgtGAGGGCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCAtt                                                                        <  1:132032/150‑1 (MQ=255)
                                                                                tGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTc                                                                      >  2:647649/1‑150 (MQ=255)
                                                                                 gAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCa                                                                     >  2:1867720/1‑150 (MQ=255)
                                                                                    gtcgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATAc                                                                  >  1:1228823/1‑150 (MQ=255)
                                                                                    gtcgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CGTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATAc                                                                  <  1:1516844/150‑1 (MQ=255)
                                                                                     tcgtcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATAcc                                                                 >  2:2366361/1‑150 (MQ=255)
                                                                                        tcCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAAt                                                              <  1:1002675/150‑1 (MQ=255)
                                                                                        gccTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAAt                                                              <  1:5939490/149‑1 (MQ=255)
                                                                                          cTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCa                                                            <  1:4791702/150‑1 (MQ=255)
                                                                                          cTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCa                                                            <  2:1900209/150‑1 (MQ=255)
                                                                                           tCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAg                                                           >  2:6755980/1‑150 (MQ=255)
                                                                                           tCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAg                                                           >  1:5541904/1‑150 (MQ=255)
                                                                                           tCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTGTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAg                                                           >  2:6777427/1‑150 (MQ=255)
                                                                                                 tGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCt                                                     >  1:6196651/1‑150 (MQ=255)
                                                                                                   ggCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTg                                                   >  2:1411519/1‑150 (MQ=255)
                                                                                                    gCACCATGGTGCCCGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGc                                                  >  2:4296835/1‑150 (MQ=255)
                                                                                                       ccATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGa                                               >  2:3420526/1‑150 (MQ=255)
                                                                                                        cATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGaa                                              <  1:3390533/150‑1 (MQ=255)
                                                                                                         aTGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAAt                                             <  2:3298844/150‑1 (MQ=255)
                                                                                                              gCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTACATACCAATCAGACCCCTTGCGGAATCTCtt                                        >  2:1281799/1‑150 (MQ=255)
                                                                                                               ccTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTc                                       <  2:450791/150‑1 (MQ=255)
                                                                                                                cTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTCTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTca                                      >  1:6598542/1‑150 (MQ=255)
                                                                                                                 tGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTcac                                     >  1:6849291/1‑150 (MQ=255)
                                                                                                                  gTCTGCGGC‑‑CTTGCGGGCAATACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTcacc                                    >  2:5272431/1‑150 (MQ=255)
                                                                                                                     tGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCaccacc                                 >  1:679096/1‑150 (MQ=255)
                                                                                                                              ttGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGcc                          >  1:4229847/1‑150 (MQ=255)
                                                                                                                                 gCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCTGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTc                        <  1:189298/150‑1 (MQ=255)
                                                                                                                                 gCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCTGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTc                        <  2:5159820/150‑1 (MQ=255)
                                                                                                                                 gCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTc                        <  1:4378603/150‑1 (MQ=255)
                                                                                                                                 gCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTc                        <  1:5749592/150‑1 (MQ=255)
                                                                                                                                      cAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACa                   >  2:1607048/1‑150 (MQ=255)
                                                                                                                                       aaGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAg                  >  1:5243190/1‑150 (MQ=255)
                                                                                                                                       aaGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAg                  >  1:3334560/1‑150 (MQ=255)
                                                                                                                                        aGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGc                 <  2:597370/150‑1 (MQ=255)
                                                                                                                                        aGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGATTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGc                 <  2:150482/150‑1 (MQ=255)
                                                                                                                                         gACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCATCACCACCCAGCGCCTCGAACAGCa                >  1:4500947/1‑150 (MQ=255)
                                                                                                                                           ccGCGCCTGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCAcc              >  1:3096054/1‑150 (MQ=255)
                                                                                                                                           ccGCGCCTGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCAcc              >  2:662231/1‑150 (MQ=255)
                                                                                                                                             gcgcCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCgg            >  1:1656493/1‑150 (MQ=255)
                                                                                                                                              cgcCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGc           <  2:1190381/150‑1 (MQ=255)
                                                                                                                                                 cTGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTc                                                                      >  2:535599/1‑88 (MQ=37)
                                                                                                                                                 cTGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTc                                                                      <  1:535599/88‑1 (MQ=37)
                                                                                                                                                 cTGCAAGTGGTAACCCCTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAgg        >  2:3931974/1‑150 (MQ=255)
                                                                                                                                                  tGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGa       >  1:9247/1‑150 (MQ=255)
                                                                                                                                                   gCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGaa      <  1:2366361/150‑1 (MQ=255)
                                                                                                                                                    cAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCACTCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGACCAGCACCGGCAGGact     >  2:6093991/1‑148 (MQ=255)
                                                                                                                                                    cAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAAt     >  1:5878140/1‑150 (MQ=255)
                                                                                                                                                     aaGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAAtctc                                          >  1:1298729/1‑112 (MQ=255)
                                                                                                                                                     aaGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAAtctc                                          <  2:1298732/112‑1 (MQ=255)
                                                                                                                                                     aaGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAAtt    >  1:1999490/1‑150 (MQ=255)
                                                                                                                                                      aGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCCATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAATTc   <  1:927906/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
AACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGC‑‑CTTGCGGGCAAGACCGCGCCTGCAAGGGGTTACCACTTGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAGGCCAGCACCTGCATTTCATACCAATCAGACCCCTTGCGGAATCTCTTCACCACCCAGCGCCTCGAACAGCACCGGCAGGAATTCG  >  NC_002947/1902045‑1902342

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