Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,287,771 C→A 26.9% intergenic (‑171/+58) PP_2010 ← / ← PP_2011 cytochrome b561/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,287,7710CA26.9% 330.9 / 147.1 193intergenic (‑171/+58)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (70/71);  new base A (32/20);  total (102/91)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

Only 100 of 195 total aligned reads displayed.

CTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGT  >  NC_002947/2287622‑2287919
                                                                                                                                                     |                                                                                                                                                     
cTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGTCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGc                                                                                                                                                       <  2:2678482/150‑1 (MQ=255)
   aaTGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtg                                                                                                                                                    <  2:3357281/150‑1 (MQ=255)
   aaTGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtg                                                                                                                                                    >  2:2091004/1‑150 (MQ=255)
     tGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg                                                                                                                                                  <  1:6276585/150‑1 (MQ=255)
     tGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACtgtg                                                                                                                                                  <  1:464641/150‑1 (MQ=255)
      ggCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGc                                                                                                                                                 <  1:2579002/150‑1 (MQ=37)
       gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa                                                                                                                                                <  1:1245148/150‑1 (MQ=255)
       gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa                                                                                                                                                <  1:6403358/150‑1 (MQ=255)
       gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa                                                                                                                                                <  1:4366038/150‑1 (MQ=255)
       gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                                                                                                <  1:1904936/150‑1 (MQ=255)
       gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCGCCCACAGTCCAACTGCACAGTCTTGGAACCCGGTGcc                                                                                                                                                >  2:864641/1‑149 (MQ=255)
       gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCa                                                                                                                                                <  2:5572969/150‑1 (MQ=37)
        cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                                                                                               <  1:1080657/150‑1 (MQ=255)
        cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                                                                                               <  1:3570959/150‑1 (MQ=255)
          ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt                                                                                                                                             <  2:2182444/150‑1 (MQ=255)
          ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt                                                                                                                                             <  2:6286536/150‑1 (MQ=255)
          ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt                                                                                                                                             <  1:2743520/150‑1 (MQ=255)
           tctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTg                                                                                                                                            <  1:5434216/150‑1 (MQ=255)
           actTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCATGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTg                                                                                                                                            <  1:3798740/149‑1 (MQ=255)
              tGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCGCCCACAGTCCAACTGCACAGTCGTGGGCGCGTGGGCCGTTGGAc                                                                                                                                         >  2:1984696/1‑150 (MQ=255)
                cTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtg                                                                                                                                       <  1:2028204/150‑1 (MQ=255)
                   gcgcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAATTGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTgg                                                                                                                                    <  2:1899493/150‑1 (MQ=255)
                     gcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGa                                                                                                                                  <  1:5974385/150‑1 (MQ=255)
                       tCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGc                                                                                                                                <  2:4650758/150‑1 (MQ=255)
                             gTCATTTTAATGGCTACCGGAAACCGCTCTTGGCGCCAGGTCTGTGAAAAATGTGTGGCCCGCCACCGGCCTCTGCGCGGGCACGCCCGCTCCCCCAGCCCACCGGCCCAGTCTTGGAAGCCTGGGCTGTTGGACTGTGGGAGCGGGCGt                                                                                                                          <  1:780961/150‑1 (MQ=25)
                             gTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGt                                                                                                                          <  2:4108538/150‑1 (MQ=255)
                              tCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTACCACAGGCCAAAACCAAAGTATTGGAATCATGTGCAGTTGGTCTGTGGGAAAGGGCGTg                                                                                                                         >  1:3360836/1‑150 (MQ=255)
                                aTTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACAGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGCGCGGGCGTGcc                                                                                                                       >  2:5986143/1‑150 (MQ=255)
                                      ttGGCTACGGTAAACAGATCTTTGCGCCTGGTCTGTTAAAATTGTGTTGCCAGCCCCGGGCCTCTTCGGGGGCACGCCCGCTCCCACCGTCCAACTGCACAGTCTTGGACGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                 <  1:5287261/149‑1 (MQ=21)
                                              cGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCGCCCAACGTCCAAAGGCCAACAGTTGGGAGCCGGTGGAGGTGGGCTGGGGGAGCGGGCGGGCCCCCGGACGGGCCgg                                                                                                         >  1:6447670/1‑150 (MQ=255)
                                               ggAAACCGTTGTTTGCGCCAGGTCTGTTAAAAATGTGTTGCCAGCCACCGGCCTCTTGGCTGGCACGCCCGCTCCCAAAGTCCAACTGCACAGTCTTGTAAGCGTGTGCTTTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGa                                                                                                       <  1:6586507/150‑1 (MQ=25)
                                               ggAAACAGATCTTTGCGCCTGGTCGGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGACGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGa                                                                                                       <  2:1373178/150‑1 (MQ=255)
                                                   aCAGATCGTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCGGTGCAGTTGGACGGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAg                                                                                                   >  2:2348479/1‑150 (MQ=255)
                                                     aGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCGCCCAAAGTCCAACTGCACAGTCGTGGCAGCCAGTGCAGTTTGACTGTGGGGGGGGGCGTGCCCCCGAAGGGGCCGGAGCAGGc                                                                                                 >  1:1562686/1‑150 (MQ=255)
                                                        tCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCt                                                                                                                                                     >  1:442901/1‑96 (MQ=255)
                                                            tGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCa                                                                                          <  2:2820442/150‑1 (MQ=255)
                                                              cgcCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGt                                                                                        <  2:471832/150‑1 (MQ=255)
                                                              cgcCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCCCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGt                                                                                        <  1:5948465/150‑1 (MQ=255)
                                                                  aGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTAcc                                                                                    <  2:4523985/150‑1 (MQ=255)
                                                                   ggTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGACGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCg                                                                                   <  1:1508276/150‑1 (MQ=255)
                                                                    gTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCGTGGAAGCCTGTGCAGTTTGACTGTGGGGGCGGGCGTGCCCGCGGAGGGGCCGGCGCAGGCCACAACAGTTACCgg                                                                                  >  2:4215205/1‑150 (MQ=25)
                                                                             aatggtgtTGCCACCCACCGGCCTCTTCGCTGGCCCGCCCGCCCCCACAGTCCAACTGCACAGTCTTGGAAGCGTGTGCAGTTGGAGTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAgg                                                                         <  1:4613025/146‑1 (MQ=25)
                                                                              attgtgtgGCCCGCCACCGGCCTCTTCGCTGGCACGCCCGCACCCCCAGTCCAACTGCACAGTCTTGGAAGCCTGTGCTGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAggg                                                                        <  1:382583/148‑1 (MQ=25)
                                                                              aaTGTGTGGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCCCCCACAGTCCAACGCAACCGTCTTGTAAGCCTGTGCTGTTGGACTGTGGGTGGGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAggg                                                                        <  1:591576/150‑1 (MQ=255)
                                                                                        ccgccaCGGGCCTCTGCGCTGGCACGCCCGCACCCACCGTCCACCGGCCCCGTCTTGTAAGCCTGTGCAGTTGGCGTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGtcatca                                                              <  1:5775558/148‑1 (MQ=255)
                                                                                               cGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACAGCACACTCGTGGAAGCCTGTGCAGTTGGACGGGGGGAGCGGGCGTGCCCGCGACGGGGCCGGAGCAGGCGACAACAGTTACCGGGCCGATCGGGCCAGGCATCACccacca                                                       >  2:3617274/1‑150 (MQ=25)
                                                                                                 gCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACAACAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGa                                                     >  2:186331/1‑150 (MQ=255)
                                                                                                   ctctTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCGTGGAAGCCTGTGCAGGTGGACTGGGGGAGCGGGCGGGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGAAGATAGGGCCAGTCCTCACCCACCAGaaa                                                   >  2:2051115/1‑150 (MQ=37)
                                                                                                    cctTCGCTGGCACGCCCGCTCCCCCAGTCCAACGGCACAGTCTTGGAAGCCTGTGCTGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAac                                                  <  1:4695071/149‑1 (MQ=255)
                                                                                                      ttCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAacac                                                <  2:3992508/150‑1 (MQ=255)
                                                                                                      ttCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAaaac                                                >  2:1437625/1‑150 (MQ=255)
                                                                                                       tCGCGGGCACGCCCGCTCCCACAGTCCAAATGCACAGGCTTCCAAGACTGTGCAGTTGGACAGTGGGAGCGGGCGTGCCAGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCACTCAACACCCACCACAAACAcc                                               >  1:3585934/1‑150 (MQ=17)
                                                                                                        cGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACAccc                                              <  2:5777813/150‑1 (MQ=255)
                                                                                                           tGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGtt                                           <  1:2090954/150‑1 (MQ=255)
                                                                                                             gCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACCGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATGAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCt                                         <  1:4528254/150‑1 (MQ=255)
                                                                                                                    cGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcctcct                                  <  2:2491005/150‑1 (MQ=255)
                                                                                                                       tCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGcc                               >  1:3296564/1‑150 (MQ=21)
                                                                                                                          cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCCGGTGGACGGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAgg                            >  1:17699/1‑150 (MQ=255)
                                                                                                                          cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAgg                            >  1:3105770/1‑150 (MQ=255)
                                                                                                                           accgTCCAACTGCACCGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGc                           <  1:1975850/147‑1 (MQ=255)
                                                                                                                           acaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGc                           >  1:5215328/1‑150 (MQ=255)
                                                                                                                             aGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCAc                         >  1:662978/1‑150 (MQ=255)
                                                                                                                                  aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGc                                                                                                 <  2:2571111/73‑1 (MQ=255)
                                                                                                                                  aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                    >  1:2492621/1‑150 (MQ=255)
                                                                                                                                  aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                    >  2:73317/1‑150 (MQ=255)
                                                                                                                                  aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                    >  1:6008126/1‑150 (MQ=255)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGTAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGATACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACAATCg                    >  1:5613615/1‑150 (MQ=37)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                    >  2:3487208/1‑150 (MQ=255)
                                                                                                                                  aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                    >  2:5620822/1‑150 (MQ=255)
                                                                                                                                    cTGCACAGGCTTCCAAGACTGTGCCGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTCCCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGc                  >  2:4400161/1‑150 (MQ=37)
                                                                                                                                     tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                 >  1:2394022/1‑150 (MQ=255)
                                                                                                                                     tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                 >  1:3906613/1‑150 (MQ=255)
                                                                                                                                     tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                 >  1:4040584/1‑150 (MQ=255)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAg                                                                                         <  1:6136160/78‑1 (MQ=14)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAg                                                                                         >  2:6136160/1‑78 (MQ=14)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                 >  1:537802/1‑150 (MQ=255)
                                                                                                                                     tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa                 >  2:810600/1‑150 (MQ=255)
                                                                                                                                      gCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt                >  2:945736/1‑150 (MQ=255)
                                                                                                                                      gCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt                >  2:2078616/1‑150 (MQ=255)
                                                                                                                                       cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc               >  1:6361874/1‑150 (MQ=255)
                                                                                                                                       cacaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGAc                                                                                              <  1:2995979/71‑1 (MQ=11)
                                                                                                                                       cacaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc               >  1:1306815/1‑150 (MQ=255)
                                                                                                                                       cacaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCAt                >  1:1982657/1‑150 (MQ=21)
                                                                                                                                        acaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCa              <  1:2569641/150‑1 (MQ=255)
                                                                                                                                        acaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCa              >  2:3313047/1‑150 (MQ=255)
                                                                                                                                        acaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCa              >  2:521323/1‑150 (MQ=255)
                                                                                                                                          cgTCTGGGACGCCGGTGCTGTTGGCGTGGGGGCGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAgg            <  1:1168800/149‑1 (MQ=255)
                                                                                                                                          aGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCCTCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCCTCAgg            >  2:3791581/1‑150 (MQ=255)
                                                                                                                                           gTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGc           >  1:5623726/1‑150 (MQ=255)
                                                                                                                                            tCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCa          >  2:6679816/1‑150 (MQ=255)
                                                                                                                                             cTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAc         <  1:4120364/150‑1 (MQ=255)
                                                                                                                                              ttGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAcc        >  2:455454/1‑150 (MQ=255)
                                                                                                                                              ttGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAcc        >  2:1862647/1‑150 (MQ=255)
                                                                                                                                              atGGAATACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGAGAAGGCGACGGAGCAGGCGACATCACATACCGGGCAGATAGGGCCAGTCATCAACCACCAGAAACACCCGTACTGCCTCCTGCAAGGAACCATCGCCATCAGGCAca        >  2:2416555/2‑149 (MQ=255)
                                                                                                                                                ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAg      >  2:1158190/1‑150 (MQ=255)
                                                                                                                                                ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAg      >  1:214367/1‑150 (MQ=255)
                                                                                                                                                 gAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACAATCGGCATCAGGAACCact     >  1:699134/1‑148 (MQ=255)
                                                                                                                                                  aaGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTc    >  2:1530348/1‑150 (MQ=255)
                                                                                                                                                    gCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGt  >  1:4921229/1‑150 (MQ=255)
                                                                                                                                                    gCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGGCATCACCCACCAGAAACACCCGGTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGt  >  2:5310513/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAGTCGT  >  NC_002947/2287622‑2287919

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