Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 830,666 A→T 8.9% intergenic (‑23/‑206) cecR ← / → rhlE DNA‑binding transcriptional dual regulator CecR/ATP‑dependent RNA helicase RhlE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913830,6660AT 8.9% 466.1 / 27.7 157intergenic (‑23/‑206)cecR/rhlEDNA‑binding transcriptional dual regulator CecR/ATP‑dependent RNA helicase RhlE
Reads supporting (aligned to +/- strand):  ref base A (84/59);  new base T (6/8);  total (90/67)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

Only 100 of 163 total aligned reads displayed.

GCCCGGCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCT  >  NC_000913/830517‑830815
                                                                                                                                                     |                                                                                                                                                     
gCCCGGCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCaaa                                                                                                                                                       >  1:675082/1‑150 (MQ=255)
gCCCGGCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCaaa                                                                                                                                                       >  1:1635664/1‑150 (MQ=255)
  ccGGCCTGGGCGGCTATCTCGCGAGTGGGGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCATCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtg                                                                                                                                                     <  1:1971915/150‑1 (MQ=255)
     gCCTGGACGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                  <  2:913882/150‑1 (MQ=255)
      ccTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgattg                                                                                                                                                 <  1:1212795/150‑1 (MQ=255)
             ggcTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGa                                                                                                                                          <  2:2117796/150‑1 (MQ=255)
              gcTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGAt                                                                                                                                         <  2:1411932/150‑1 (MQ=255)
               cTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaat                                                                                                                                        <  1:2332303/150‑3 (MQ=255)
               cTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGAtt                                                                                                                                        <  1:1947498/150‑1 (MQ=255)
                tATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGTTTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaatg                                                                                                                                       <  1:853856/150‑4 (MQ=255)
                tATCTCGCCAGGGCTGGCGTTCATTCCATAATCACCAAACAGCTCCATTGAGTCAGCAATCACCTGTTTTTTCTCCTGTTCACCCATGATTGTCATGGCAGGATTATTCTTCTCACAGTCCATTCTTAATCAATCATTTGATTAAGaatg                                                                                                                                       <  2:1715802/150‑4 (MQ=255)
                 aTCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaatgg                                                                                                                                      <  1:950811/150‑5 (MQ=255)
                   ctcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGa                                                                                                                                    >  1:1860387/1‑150 (MQ=255)
                    tcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGAc                                                                                                                                   <  1:2162577/150‑1 (MQ=255)
                      gcgAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTc                                                                                                                                 <  2:28661/150‑1 (MQ=255)
                               ggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGa                                                                                                                        >  2:2077055/1‑150 (MQ=255)
                               ggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGa                                                                                                                        >  2:673893/1‑150 (MQ=255)
                                gCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAg                                                                                                                       <  1:2373138/150‑1 (MQ=255)
                                 cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                      >  2:2241512/1‑150 (MQ=255)
                                 cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                      >  1:997446/1‑150 (MQ=255)
                                    tCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGt                                                                                                                   <  1:938365/150‑1 (MQ=255)
                                    tCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGt                                                                                                                   >  2:835334/1‑150 (MQ=255)
                                     cATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTc                                                                                                                  <  1:1250183/150‑1 (MQ=255)
                                          cATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTa                                                                                                             >  1:616138/1‑150 (MQ=255)
                                            tatTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATg                                                                                                           <  1:1420251/150‑1 (MQ=255)
                                             atTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATgg                                                                                                          <  1:2108472/150‑1 (MQ=255)
                                               tCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCt                                                                                                        >  2:1053785/1‑150 (MQ=255)
                                                 aCCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTCTCATGGCAGGATTAATCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTACAGTATGGCTaa                                                                                                      >  2:928721/1‑150 (MQ=255)
                                                  ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAg                                                                                                     <  2:629359/150‑1 (MQ=255)
                                                   cAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGa                                                                                                    >  1:1939110/1‑150 (MQ=255)
                                                   aaaaCTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGa                                                                                                    >  2:1839239/2‑150 (MQ=255)
                                                       cTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAAttt                                                                                                >  2:941462/1‑150 (MQ=255)
                                                         gcgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTa                                                                                              <  2:1210541/150‑1 (MQ=255)
                                                         gcgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCATGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTa                                                                                              <  2:1633944/150‑1 (MQ=255)
                                                          cgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAg                                                                                             >  1:2468087/1‑150 (MQ=255)
                                                           gcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGc                                                                                            <  1:2044172/150‑1 (MQ=255)
                                                            ccAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCa                                                                                           <  1:1561275/150‑1 (MQ=255)
                                                            ccAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCa                                                                                           <  2:792344/150‑1 (MQ=255)
                                                              aGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAAc                                                                                         >  1:1498094/1‑150 (MQ=255)
                                                              aGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAAc                                                                                         >  1:924632/1‑150 (MQ=255)
                                                                 gCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGcc                                                                                      >  1:337189/1‑150 (MQ=255)
                                                                 gCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGcc                                                                                      <  2:87755/150‑1 (MQ=255)
                                                                   ggCATCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGCCCAGTATGGCTAAGAATTTTAGCAACGCCAg                                                                                    <  1:869325/150‑1 (MQ=255)
                                                                              gCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATa                                                                         >  1:506709/1‑150 (MQ=255)
                                                                                tGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAAt                                                                       >  2:1281811/1‑150 (MQ=255)
                                                                                 gTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTAATCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAAtt                                                                      >  2:1873024/1‑150 (MQ=255)
                                                                                   tttttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAAACAAATGATTGATTAAGATTATGACTCCATAGGGGAGATGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTa                                                                    >  1:684588/1‑150 (MQ=255)
                                                                                    ttttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTAt                                                                   <  2:616138/150‑1 (MQ=255)
                                                                                      ttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACGCCATAGGGGAGTTGTCAAGTATGGCTAAGAATTTTAGAAACGCCAGTAACAGGGATAATTTATgg                                                                 >  2:452086/1‑149 (MQ=255)
                                                                                        cGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATgcgc                                                               >  1:2427492/1‑150 (MQ=255)
                                                                                        cGCCTGTTCACCCTTGATTGTCATGGAAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAACATTATGACTCCATAGGGGACTTGTCCAGTAAGGCTAAGAATATTAGAAAAGCAAGTCACAGGGATAAATTATgagc                                                               >  1:511233/1‑150 (MQ=255)
                                                                                         gCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCt                                                              >  1:922633/1‑150 (MQ=255)
                                                                                         gCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGAAACGCCAGTCACAGGGATAATTTATGCGCt                                                              >  2:1475573/1‑150 (MQ=255)
                                                                                          ccTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTg                                                             >  1:328907/1‑150 (MQ=255)
                                                                                            tGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCg                                                           >  1:1926202/1‑150 (MQ=255)
                                                                                                 aCCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTcaca                                                      <  1:1693964/150‑1 (MQ=255)
                                                                                                  cccTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACaa                                                     >  2:2262947/1‑150 (MQ=255)
                                                                                                   ccTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACaaa                                                    >  1:688076/1‑150 (MQ=255)
                                                                                                     ttGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACaaaaa                                                  <  2:1308069/150‑1 (MQ=255)
                                                                                                       gATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACt                                                <  2:677653/150‑1 (MQ=255)
                                                                                                          tGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCt                                             >  2:37970/1‑150 (MQ=255)
                                                                                                           gTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTa                                            >  2:523771/1‑150 (MQ=255)
                                                                                                           gTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTa                                            <  2:922633/150‑1 (MQ=255)
                                                                                                           gTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTa                                            >  2:985243/1‑150 (MQ=255)
                                                                                                            tCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTac                                           <  1:659444/150‑1 (MQ=255)
                                                                                                             cATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGAGTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCAAGTCACAGGTATAATTTATGCGCTGCGTAACAAAAACTGCTaaa                                          >  1:1848762/1‑148 (MQ=255)
                                                                                                               tGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACt                                        >  1:42058/1‑150 (MQ=255)
                                                                                                                ggCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTc                                       >  2:460448/1‑150 (MQ=255)
                                                                                                                 gCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                      >  2:1673292/1‑150 (MQ=255)
                                                                                                                 gCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                      >  2:1815030/1‑150 (MQ=255)
                                                                                                                 gCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                      >  2:1993911/1‑150 (MQ=255)
                                                                                                                  cAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCg                                     <  2:14654/150‑1 (MQ=255)
                                                                                                                  cAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCg                                     <  1:773123/150‑1 (MQ=255)
                                                                                                                   aGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGc                                    >  2:234410/1‑150 (MQ=255)
                                                                                                                     gATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTc                                  <  2:1869296/150‑1 (MQ=255)
                                                                                                                     gATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTc                                  >  2:931500/1‑150 (MQ=255)
                                                                                                                        tattCATCGCACAGTCAATTCTTAATCAAATGAATGATTAAGATTATGACTCCAAAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAAATGCTACACTCCGCTccca                               >  2:1067741/1‑149 (MQ=255)
                                                                                                                          ttCATCGCACAGTCCATTTTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCa                             <  2:1232312/150‑1 (MQ=255)
                                                                                                                            cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATTCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATg                           >  1:1540217/1‑150 (MQ=255)
                                                                                                                               cGCACAGTCAATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCAGCGTAACAAAAACTGCTACACTCCGCTCCATCATGACa                        >  2:1856221/1‑150 (MQ=255)
                                                                                                                                gCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACAt                       >  2:1544435/1‑150 (MQ=255)
                                                                                                                                gCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACAt                       >  1:431635/1‑150 (MQ=255)
                                                                                                                                 tctaGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGGCCAGAATGGCTAAGAAATTTAGAAACGCAAGTCAAAGGGATAATTTATGAGCTGAGTCAAAAAAACTGCTACACTCCGCTCCCTCATGACAtt                      <  2:1079661/147‑1 (MQ=255)
                                                                                                                                  acaGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATtg                     <  1:523771/150‑1 (MQ=255)
                                                                                                                                    aGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATtgtg                   <  1:1452299/150‑1 (MQ=255)
                                                                                                                                       ataatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtt                >  1:2127760/5‑150 (MQ=255)
                                                                                                                                         aatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtttt              >  2:1674878/3‑150 (MQ=255)
                                                                                                                                         aatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtttt              >  1:485887/3‑150 (MQ=255)
                                                                                                                                          ttCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGttttt             >  1:126685/1‑150 (MQ=255)
                                                                                                                                          ttCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGttttt             <  1:1873024/150‑1 (MQ=255)
                                                                                                                                          atCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGttttt             >  1:1028760/2‑150 (MQ=255)
                                                                                                                                          atCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGttttt             >  1:959172/2‑150 (MQ=255)
                                                                                                                                           tCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTg            <  2:2356065/150‑1 (MQ=255)
                                                                                                                                            cTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGt           >  2:2045409/1‑150 (MQ=255)
                                                                                                                                              tAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCa         >  1:28498/1‑150 (MQ=255)
                                                                                                                                               aaTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAt        >  1:1456077/1‑150 (MQ=255)
                                                                                                                                                 tCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAttt      >  2:61358/1‑150 (MQ=255)
                                                                                                                                                  cAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTAAtttt     >  2:1572347/1‑150 (MQ=255)
                                                                                                                                                   aaaTGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTc    >  1:1743495/1‑150 (MQ=255)
                                                                                                                                                    aaTGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTcc   <  2:2468087/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCGGCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCT  >  NC_000913/830517‑830815

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