Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 496,099 T→A 8.5% intergenic (‑63/‑33) PP_0408 ← / → PP_0409 hypothetical protein/two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947496,0990TA 8.5% 320.2 / 14.8 105intergenic (‑63/‑33)PP_0408/PP_0409hypothetical protein/two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (50/46);  new base A (6/3);  total (56/49)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 115 total aligned reads displayed.

TTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCTG  >  NC_002947/495950‑496248
                                                                                                                                                     |                                                                                                                                                     
tttttCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgct                                                                                                                                                       <  1:1095988/150‑1 (MQ=255)
  tttCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgctgc                                                                                                                                                     <  1:1098060/150‑1 (MQ=255)
  tttCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgctgc                                                                                                                                                     <  1:1098081/150‑1 (MQ=255)
  tttCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgctgc                                                                                                                                                     <  1:2358512/150‑1 (MQ=255)
    tgggCATCTGCGGCCAGGGCGCCGAGTTGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCAt                                                                                                                                                   <  2:947608/148‑1 (MQ=255)
    tCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCAt                                                                                                                                                   <  1:1754430/150‑1 (MQ=255)
     tccactCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATg                                                                                                                                                  >  1:1530037/6‑150 (MQ=255)
     cGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATg                                                                                                                                                  >  1:30554/1‑150 (MQ=255)
       cactCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCCTg                                                                                                                                                  >  2:1087429/4‑148 (MQ=255)
       cactCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATg                                                                                                                                                  <  1:1087429/145‑1 (MQ=255)
        actCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCTGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCGGGATGTGAGTGGCCCGTGACGGCCATGCTGCATGGca                                                                                                                                               >  2:1592251/3‑149 (MQ=255)
        actCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc                                                                                                                                               <  2:1530037/148‑1 (MQ=255)
        actCGGCGCCCTGGCCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGCCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATCATGCATGGcc                                                                                                                                               >  1:2243301/3‑150 (MQ=255)
             gcggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcac                                                                                                                                          <  1:1428222/150‑3 (MQ=255)
              cggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcacg                                                                                                                                         <  2:30554/150‑4 (MQ=255)
                 ccAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGcc                                                                                                                                      <  2:293066/150‑1 (MQ=255)
                 ccAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGcc                                                                                                                                      <  1:108263/150‑1 (MQ=255)
                   aGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCaa                                                                                                                                    <  2:794670/150‑1 (MQ=255)
                   aGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCaa                                                                                                                                    <  2:690801/150‑1 (MQ=255)
                        gcCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGa                                                                                                                               >  1:2271747/1‑150 (MQ=255)
                          cGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                             >  1:321847/1‑150 (MQ=255)
                               ggAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGtata                                                                                                                        <  1:932239/150‑1 (MQ=255)
                                gAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCttt                                                                                                                                         >  2:50590/1‑132 (MQ=255)
                                 aaGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGa                                                                                                                      >  1:1328620/1‑150 (MQ=255)
                                   ggACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATg                                                                                                                    >  2:1352090/1‑150 (MQ=255)
                                   ggACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATg                                                                                                                    <  1:1352090/150‑1 (MQ=255)
                                   ggACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATg                                                                                                                    <  1:1351065/150‑1 (MQ=255)
                                   ggACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATg                                                                                                                    >  2:2377359/1‑150 (MQ=255)
                                    gACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGt                                                                                                                                  <  2:2191542/135‑1 (MQ=255)
                                    gACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGt                                                                                                                                  >  1:2191542/1‑135 (MQ=255)
                                      ccGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGCTGAAg                                                                                                                 >  1:1855867/1‑150 (MQ=255)
                                         aGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGa                                                                                                              <  2:2108519/150‑1 (MQ=255)
                                          ggCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCCAGTCGAGCGTATAGATGAAGCGAg                                                                                                             >  2:318466/1‑150 (MQ=255)
                                             cAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGc                                                                                                          >  1:486446/1‑150 (MQ=255)
                                             cAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGc                                                                                                          >  1:527429/1‑150 (MQ=255)
                                                gcagcaCATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCGGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTc                                                                                                       >  2:1485713/1‑150 (MQ=255)
                                                gcagcaCATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTc                                                                                                       <  2:17230/150‑1 (MQ=255)
                                                    cacaTTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCt                                                                                                   <  1:2162408/150‑1 (MQ=255)
                                                    cacaTTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCt                                                                                                   <  1:1555744/150‑1 (MQ=255)
                                                              gCGTCGTGGCATAAAGGGTGGACATGAGTCACCACACGGCCAAGATGAAGCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGACATAGCGAGTGCGGCGCCTGTTGGTTATc                                                                                         <  1:563694/150‑1 (MQ=255)
                                                                    tGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGgctgc                                                                                   <  1:1531131/150‑1 (MQ=255)
                                                                       gATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCAGGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGATATCGGCTGCTtg                                                                                >  2:1299554/1‑150 (MQ=255)
                                                                        aTAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGa                                                                                                                               >  2:1542483/1‑102 (MQ=255)
                                                                        aTAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGa                                                                                                                               <  1:1542483/102‑1 (MQ=255)
                                                                         tAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTtgtg                                                                              >  2:164336/1‑150 (MQ=255)
                                                                              tgtgGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGc                                                                         >  2:1899407/1‑150 (MQ=255)
                                                                              tgtgGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGc                                                                         <  2:1615958/150‑1 (MQ=255)
                                                                                  gACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCtt                                                                     >  2:908579/1‑150 (MQ=255)
                                                                                  aaCATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGGCTAGATGTCCTTGCCAATTCATGGGGTTAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCtt                                                                     <  1:1485713/149‑1 (MQ=255)
                                                                                   aCATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGgct                                                                                     >  2:2023167/1‑133 (MQ=255)
                                                                                      tGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGAtt                                                                 <  2:486446/150‑1 (MQ=255)
                                                                                            aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTTGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTgg                                                           >  2:535630/1‑150 (MQ=255)
                                                                                            aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTgg                                                           <  1:1299554/150‑1 (MQ=255)
                                                                                                 aaGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCaa                                                      >  1:217389/1‑150 (MQ=255)
                                                                                                   ggCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTgg                                                                                              >  1:2357017/1‑108 (MQ=255)
                                                                                                     ccANAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGccc                                                  >  2:1506831/1‑150 (MQ=255)
                                                                                                     ccAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTACCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCg                                                                                        >  2:839659/1‑112 (MQ=255)
                                                                                                     ccAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTACCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCg                                                                                        <  1:839659/112‑1 (MQ=255)
                                                                                                      cAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCg                                                 >  1:2182944/1‑150 (MQ=255)
                                                                                                       aGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTg                                                                                                                                        >  1:2035080/1‑62 (MQ=255)
                                                                                                       aGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTg                                                                                                                                        <  2:2035080/62‑1 (MQ=255)
                                                                                                          aTGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGgc                                             >  1:932038/1‑150 (MQ=255)
                                                                                                              aaCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTg                                         >  1:2206770/1‑150 (MQ=255)
                                                                                                               aCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGc                                        >  1:2239430/1‑150 (MQ=255)
                                                                                                                ccGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGTGAGTGCGTCGCCTGTTGGTAATCGACTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGcc                                       >  1:1984412/1‑150 (MQ=255)
                                                                                                                 cGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGccc                                      >  2:282787/1‑150 (MQ=255)
                                                                                                                   gCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCctc                                    >  2:851039/1‑150 (MQ=255)
                                                                                                                     aGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGc                                                         >  2:2078556/1‑127 (MQ=255)
                                                                                                                     aGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGc                                                         <  1:2078556/127‑1 (MQ=255)
                                                                                                                      gCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGACATGTGCAATTCGAGCGTATAGATGAAGCACTTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgt                                 <  1:2377359/150‑1 (MQ=255)
                                                                                                                       ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTTCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgtg                                >  2:410293/1‑150 (MQ=255)
                                                                                                                       ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgtg                                <  1:1058849/150‑1 (MQ=255)
                                                                                                                       ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgtg                                >  1:1256270/1‑150 (MQ=255)
                                                                                                                             aGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTg                          <  1:282787/150‑1 (MQ=255)
                                                                                                                               tATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGa                        <  2:932038/150‑1 (MQ=255)
                                                                                                                               tATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGa                        <  2:2206770/150‑1 (MQ=255)
                                                                                                                                 tttGCCCGTGACCGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGCTATCGGCAGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                      <  2:1328620/150‑1 (MQ=255)
                                                                                                                                   tGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCt                    <  2:527429/150‑1 (MQ=255)
                                                                                                                                      ccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc                 >  1:814609/1‑150 (MQ=255)
                                                                                                                                      ccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc                 >  1:494300/1‑150 (MQ=255)
                                                                                                                                      ccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc                 <  2:844559/150‑1 (MQ=255)
                                                                                                                                      acGTGACAGACATGCTACATGGCCGACTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGGTCAAATCCAGAACCCAGGCCAACCCCGAGGCGCGTCCCTCTGTGGTGCGGGAGCCTGAc                 >  2:1106080/2‑150 (MQ=255)
                                                                                                                                        gTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGCAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACca               <  1:1106080/150‑1 (MQ=255)
                                                                                                                                         agACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcag              >  1:1939899/2‑150 (MQ=255)
                                                                                                                                         agACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcag              >  2:146554/2‑150 (MQ=255)
                                                                                                                                          gACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGGGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTTATTGACTGGGCCAATCCCGAGTCGCTGCCCTCTGTGGTGCTGTAGCCTGACcagc             >  1:638394/1‑150 (MQ=255)
                                                                                                                                           aCGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTCTTGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagca            <  1:1592251/150‑1 (MQ=255)
                                                                                                                                            cGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcag           >  2:1650009/1‑150 (MQ=255)
                                                                                                                                             ggCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcagc          >  1:137869/1‑150 (MQ=255)
                                                                                                                                              gCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGa                                                                   <  1:303827/92‑1 (MQ=255)
                                                                                                                                              gCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGa                                                                   >  2:303827/1‑92 (MQ=255)
                                                                                                                                              gCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcagca         <  2:1855867/150‑1 (MQ=255)
                                                                                                                                              gCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcagca         >  1:347454/1‑150 (MQ=255)
                                                                                                                                                cATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGt       <  1:164336/150‑1 (MQ=255)
                                                                                                                                                cATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGt       <  1:56870/150‑1 (MQ=255)
                                                                                                                                                  tgctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTgg     <  2:494300/150‑1 (MQ=255)
                                                                                                                                                  tgctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTgg     <  2:1561509/150‑1 (MQ=255)
                                                                                                                                                   gctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGc    >  1:52581/1‑150 (MQ=255)
                                                                                                                                                    ctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCt   <  1:1506831/150‑1 (MQ=255)
                                                                                                                                                     tgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCTg  <  1:1899407/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCTG  >  NC_002947/495950‑496248

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