Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,048,928 C→A 93.3% E903* (GAA→TAA)  torS ← hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,048,9280CA93.3% 852.4 / 46.6 270E903* (GAA→TAA) torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base C (9/9);  new base A (126/126);  total (135/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 270 total aligned reads displayed.

GGCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTT  >  CP009273/1048781‑1049075
                                                                                                                                                   |                                                                                                                                                   
ggCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACa                                                                                                                                                   <  1:3897415/150‑1 (MQ=255)
ggCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACa                                                                                                                                                   >  1:201004/1‑150 (MQ=255)
  ccgccgCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACaga                                                                                                                                                 >  1:970880/1‑150 (MQ=255)
   cgccgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACagag                                                                                                                                                >  2:1064015/1‑150 (MQ=255)
   cgccgcTGCCAGCGGGTTCTTTTAGTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACagag                                                                                                                                                <  2:1553711/150‑1 (MQ=255)
      cgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCgg                                                                                                                                             >  2:1463708/1‑150 (MQ=255)
      cgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCgg                                                                                                                                             >  1:2118420/1‑150 (MQ=255)
      cacTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTGTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCgg                                                                                                                                             <  2:4314542/148‑1 (MQ=255)
       gcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGc                                                                                                                                            >  1:1726766/1‑150 (MQ=255)
        cTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCa                                                                                                                                           <  2:3682487/150‑1 (MQ=255)
         tGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCaa                                                                                                                                          <  2:844475/150‑1 (MQ=255)
             aGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACt                                                                                                                                      <  1:2700414/150‑1 (MQ=255)
                   ttCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgtgt                                                                                                                                <  1:1792947/150‑1 (MQ=255)
                   ttCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgtgt                                                                                                                                >  2:2680312/1‑150 (MQ=255)
                   ttCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgtgt                                                                                                                                <  2:2194729/150‑1 (MQ=255)
                       tttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAAtt                                                                                                                            <  2:4487731/150‑1 (MQ=255)
                        ttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAAttt                                                                                                                           <  2:2529046/150‑1 (MQ=255)
                         tATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttc                                                                                                                          >  2:610225/1‑150 (MQ=255)
                            tctcTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttcttc                                                                                                                       >  1:4092405/1‑150 (MQ=255)
                             ctctGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCg                                                                                                                      >  2:4452665/1‑150 (MQ=255)
                                tGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCgtgt                                                                                                                   <  1:2354078/150‑1 (MQ=255)
                                tGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCATGTAATTTCTTCgtgt                                                                                                                   >  1:1651771/1‑150 (MQ=255)
                                  ccATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTgg                                                                                                                 >  1:2615661/1‑150 (MQ=255)
                                  ccATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTgg                                                                                                                 >  2:559446/1‑150 (MQ=255)
                                    aTGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                               <  1:2174990/150‑1 (MQ=255)
                                    aTGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                               <  1:2174843/150‑1 (MQ=255)
                                    aTGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                               <  1:3975481/150‑1 (MQ=255)
                                          ttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAgggggg                                                                                                         >  2:3756828/1‑150 (MQ=255)
                                          ttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAgggggg                                                                                                         >  2:333230/1‑150 (MQ=255)
                                            gCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCa                                                                                                       <  2:3159003/150‑1 (MQ=255)
                                            gCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCa                                                                                                       >  1:4474042/1‑150 (MQ=255)
                                             ccAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAg                                                                                                      >  1:1399559/1‑150 (MQ=255)
                                              cAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGa                                                                                                     >  1:847818/1‑150 (MQ=255)
                                              cAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGa                                                                                                     >  1:1130320/1‑150 (MQ=255)
                                              cAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGa                                                                                                     >  1:1914419/1‑150 (MQ=255)
                                                gTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATa                                                                                                   <  2:3819071/150‑1 (MQ=255)
                                                 tCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATaa                                                                                                  <  2:4488313/150‑1 (MQ=255)
                                                   cGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATg                                                                                                >  2:2001799/1‑150 (MQ=255)
                                                   cGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATg                                                                                                >  1:3443375/1‑150 (MQ=255)
                                                   cGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATg                                                                                                >  1:4375124/1‑150 (MQ=255)
                                                    gTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGGGGACGGGGGGCAGATACTgg                                                                                               >  1:1313774/1‑150 (MQ=255)
                                                      tAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgct                                                                                             >  2:2083197/1‑150 (MQ=255)
                                                      tAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgct                                                                                             <  1:431039/150‑1 (MQ=255)
                                                            gTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgctgct                                                                                       <  2:910488/150‑1 (MQ=255)
                                                             tataCCACATGCCACGATTCGAGTGGTTCTGAACTGGCTGATATTGATTCAACTCGCGTTCAGGTCTTTCGTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTc                                                                                      <  2:3439045/150‑1 (MQ=255)
                                                               taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCa                                                                                    <  1:1162418/150‑1 (MQ=255)
                                                                 ccAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGt                                                                                  >  2:4705480/1‑150 (MQ=255)
                                                                         cACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                                          >  2:1792911/1‑150 (MQ=255)
                                                                          aCTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                                         >  1:1528440/1‑150 (MQ=255)
                                                                                 tAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                                  <  2:3259753/150‑1 (MQ=255)
                                                                                    ttgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGctggctgg                                                               <  1:4321969/150‑1 (MQ=255)
                                                                                     tgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                              >  2:2513526/1‑150 (MQ=255)
                                                                                     tgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                              <  1:30710/150‑1 (MQ=255)
                                                                                     tgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                              >  1:4508508/1‑150 (MQ=255)
                                                                                         cTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGata                                                          <  1:269/150‑1 (MQ=255)
                                                                                          ttAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                                         <  1:1411105/150‑1 (MQ=255)
                                                                                          ttAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                                         >  1:251766/1‑150 (MQ=255)
                                                                                             aCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCa                                                      <  1:528589/150‑1 (MQ=255)
                                                                                                ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTc                                                   >  1:2891409/1‑150 (MQ=255)
                                                                                                 gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                                  <  1:1254111/150‑1 (MQ=255)
                                                                                                  cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTccc                                                 <  2:3982840/150‑1 (MQ=255)
                                                                                                    gATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                               >  2:1804829/1‑150 (MQ=255)
                                                                                                     atatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGa                                              >  1:1182999/1‑150 (MQ=255)
                                                                                                     atatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGa                                              >  2:2415009/1‑150 (MQ=255)
                                                                                                      tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGAc                                             >  1:2510435/1‑150 (MQ=255)
                                                                                                          gATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGa                                         >  1:3506256/1‑150 (MQ=255)
                                                                                                          gATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGa                                         >  1:3845326/1‑150 (MQ=255)
                                                                                                           aTTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                                        <  1:63439/150‑1 (MQ=255)
                                                                                                            ttCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                                       <  2:1091034/150‑1 (MQ=255)
                                                                                                            ttCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                                       >  2:2594453/1‑150 (MQ=255)
                                                                                                             tCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                                      <  1:2386857/150‑1 (MQ=255)
                                                                                                              cAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCaa                                     >  2:3026545/1‑150 (MQ=255)
                                                                                                              cAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCaa                                     <  1:323823/150‑1 (MQ=255)
                                                                                                                aaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                                   <  2:1634522/150‑1 (MQ=255)
                                                                                                                      gTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                             >  1:3274174/1‑150 (MQ=255)
                                                                                                                        tCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTa                           <  1:4610520/150‑1 (MQ=255)
                                                                                                                          aGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc                         >  2:2933921/1‑150 (MQ=255)
                                                                                                                           ggTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt                        <  1:1064015/150‑1 (MQ=255)
                                                                                                                           ggTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt                        <  1:4452665/150‑1 (MQ=255)
                                                                                                                            gTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg                       <  2:2910235/150‑1 (MQ=255)
                                                                                                                            gTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg                       >  1:4200530/1‑150 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                     >  2:3930307/1‑150 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                     >  2:2876481/1‑150 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                     >  2:4235817/1‑150 (MQ=255)
                                                                                                                                ttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgc                   >  2:1721026/1‑150 (MQ=255)
                                                                                                                                  cttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcgc                 <  2:4044403/150‑1 (MQ=255)
                                                                                                                                   ttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGcgcg                <  1:559446/150‑1 (MQ=255)
                                                                                                                                    tATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg               >  1:3813380/1‑150 (MQ=255)
                                                                                                                                      tGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCa             >  1:822885/1‑150 (MQ=255)
                                                                                                                                        cAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACg           >  2:4581296/1‑150 (MQ=255)
                                                                                                                                        cAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACg           >  2:1233550/1‑150 (MQ=255)
                                                                                                                                             aTGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTa      <  1:3056615/150‑1 (MQ=255)
                                                                                                                                              tGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAt     >  2:3621121/1‑150 (MQ=255)
                                                                                                                                               gCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtt    <  1:3170350/150‑1 (MQ=255)
                                                                                                                                               gCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtt    <  1:1423627/150‑1 (MQ=255)
                                                                                                                                               gCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtt    <  2:2372059/150‑1 (MQ=255)
                                                                                                                                               gCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtt    >  2:1030553/1‑150 (MQ=255)
                                                                                                                                                cTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAttt   >  2:2989026/1‑150 (MQ=255)
                                                                                                                                                 ttCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtttt  >  2:466900/1‑150 (MQ=255)
                                                                                                                                                 ttCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtttt  >  2:2810989/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTT  >  CP009273/1048781‑1049075

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: