Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 830,669 A→T 9.0% intergenic (‑26/‑203) 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,6690AT 9.0% 462.7 / 28.0 156intergenic (‑26/‑203)cecR/rhlEDNA‑binding transcriptional dual regulator CecR/ATP‑dependent RNA helicase RhlE
Reads supporting (aligned to +/- strand):  ref base A (85/57);  new base T (6/8);  total (91/65)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 163 total aligned reads displayed.

GCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCTTTT  >  NC_000913/830522‑830818
                                                                                                                                                   |                                                                                                                                                     
 ccTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgattg                                                                                                                                                    <  1:1212795/150‑1 (MQ=255)
  cTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgattga                                                                                                                                                   <  2:2049266/150‑1 (MQ=255)
   tGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgattgat                                                                                                                                                  >  2:1032198/1‑150 (MQ=255)
     ggcggcTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTa                                                                                                                                                >  1:677653/1‑150 (MQ=255)
      gcggcTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTaa                                                                                                                                               <  2:2105320/150‑1 (MQ=255)
       cggcTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAg                                                                                                                                              <  1:57726/150‑1 (MQ=255)
          cTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaat                                                                                                                                           <  1:2332303/150‑3 (MQ=255)
            aTCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTTTTTTTTCGCCTGTTCACCCTTGATTGTCATGGAAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTAt                                                                                                                                         <  1:1734343/150‑1 (MQ=255)
              ctcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGa                                                                                                                                       >  2:659444/1‑150 (MQ=255)
              ctcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGa                                                                                                                                       >  1:1860387/1‑150 (MQ=255)
               tcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGAc                                                                                                                                      <  1:2419794/150‑1 (MQ=255)
               tcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGAc                                                                                                                                      <  1:2162577/150‑1 (MQ=255)
                       ggtggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAgg                                                                                                                              >  2:256473/1‑150 (MQ=255)
                       ggtggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAgg                                                                                                                              <  1:949961/150‑1 (MQ=255)
                       ggtagCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAgg                                                                                                                              >  2:850757/1‑150 (MQ=255)
                          ggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGa                                                                                                                           >  1:1939701/1‑150 (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)
                                  ttCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCa                                                                                                                   <  2:1486558/150‑1 (MQ=255)
                                     cATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGtt                                                                                                                >  2:1548920/1‑149 (MQ=255)
                                     cATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTa                                                                                                                >  1:616138/1‑150 (MQ=255)
                                     cATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTa                                                                                                                <  2:1860387/150‑1 (MQ=255)
                                       tatTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATg                                                                                                              >  2:1693964/1‑150 (MQ=255)
                                       tatTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATg                                                                                                              <  1:1420251/150‑1 (MQ=255)
                                        atTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATgg                                                                                                             <  1:2108472/150‑1 (MQ=255)
                                          tCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCt                                                                                                           >  2:1244427/1‑150 (MQ=255)
                                           cACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTa                                                                                                          >  1:316708/1‑150 (MQ=255)
                                           cACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTa                                                                                                          >  1:387078/1‑150 (MQ=255)
                                           cACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTa                                                                                                          >  1:1869296/1‑150 (MQ=255)
                                            aCCAAACTGCGCCAGTGCGGCAGCAATCAGCTCTTTTTTCGCCTGTTCAAACTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTaa                                                                                                         >  2:773123/1‑150 (MQ=255)
                                             ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAGGATTATGAGTCCATAGGGGAGTTGTCCAGTTTGGGTAAg                                                                                                        >  1:2135023/1‑150 (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)
                                                     cgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAg                                                                                                >  1:1843709/1‑150 (MQ=255)
                                                      gcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGc                                                                                               <  1:2044172/150‑1 (MQ=255)
                                                       ccAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCa                                                                                              <  1:1561275/150‑1 (MQ=255)
                                                         aGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAAc                                                                                            >  1:924632/1‑150 (MQ=255)
                                                         aGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAAc                                                                                            >  1:1498094/1‑150 (MQ=255)
                                                            gCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGcc                                                                                         <  2:1635664/150‑1 (MQ=255)
                                                              ggCATCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGCCCAGTATGGCTAAGAATTTTAGCAACGCCAg                                                                                       <  1:869325/150‑1 (MQ=255)
                                                              ggCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAg                                                                                       >  2:926490/1‑150 (MQ=255)
                                                                cagcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTc                                                                                     <  2:961786/150‑1 (MQ=255)
                                                                  gcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTcac                                                                                   <  2:1413458/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)
                                                                               ttttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTAt                                                                      <  2:616138/150‑1 (MQ=255)
                                                                                  tCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATgcg                                                                   <  1:1874104/150‑1 (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)
                                                                                            aCCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTcaca                                                         <  1:1693964/150‑1 (MQ=255)
                                                                                             cccTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACaa                                                        <  2:931036/150‑1 (MQ=255)
                                                                                              ccTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACaaa                                                       >  1:688076/1‑150 (MQ=255)
                                                                                                  gATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACt                                                   <  2:677653/150‑1 (MQ=255)
                                                                                                     tGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCt                                                >  2:37970/1‑150 (MQ=255)
                                                                                                      gTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTa                                               >  2:626959/1‑150 (MQ=255)
                                                                                                      gTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTa                                               <  2:922633/150‑1 (MQ=255)
                                                                                                      gTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTa                                               >  2:1757855/1‑150 (MQ=255)
                                                                                                          tGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACt                                           >  1:42058/1‑150 (MQ=255)
                                                                                                           ggCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTc                                          >  2:460448/1‑150 (MQ=255)
                                                                                                           ggCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTc                                          >  2:1940466/1‑150 (MQ=255)
                                                                                                            gCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                         >  1:50264/1‑150 (MQ=255)
                                                                                                             cAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCg                                        <  1:773123/150‑1 (MQ=255)
                                                                                                             cAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCg                                        >  2:600882/1‑150 (MQ=255)
                                                                                                               ggATTATTCATCGCTCAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATATGGGAGTTGTCCAGTATCGCTAAGAATTTAAGCAACGCCAGTCACAGGGATAATTTATGCGATGCCTCACAAAAACTGCAACACTCAGCt                                      >  1:2018726/1‑150 (MQ=255)
                                                                                                               ggATTATTCATCGCAGAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCt                                      <  2:2427638/150‑1 (MQ=255)
                                                                                                                   tattCATCGCACAGTCAATTCTTAATCAAATGAATGATTAAGATTATGACTCCAAAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAAATGCTACACTCCGCTccca                                  >  2:1067741/1‑149 (MQ=255)
                                                                                                                       cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATTCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATg                              >  1:1540217/1‑150 (MQ=255)
                                                                                                                       cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCCCAAAAACTGCTACACTCCGCTCCCTCATg                              >  2:1334314/1‑150 (MQ=255)
                                                                                                                       cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATg                              <  2:2261624/150‑1 (MQ=255)
                                                                                                                         tCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGAc                            <  2:1843709/150‑1 (MQ=255)
                                                                                                                           gCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACAt                          >  1:431635/1‑150 (MQ=255)
                                                                                                                            tctaGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGGCCAGAATGGCTAAGAAATTTAGAAACGCAAGTCAAAGGGATAATTTATGAGCTGAGTCAAAAAAACTGCTACACTCCGCTCCCTCATGACAtt                         <  2:1079661/147‑1 (MQ=255)
                                                                                                                            cacaGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACAtt                         >  2:420087/1‑150 (MQ=255)
                                                                                                                             acaGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATtg                        <  2:2287805/150‑1 (MQ=255)
                                                                                                                               aGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATtgtg                      <  1:1452299/150‑1 (MQ=255)
                                                                                                                                    aatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtttt                 >  2:1674878/3‑150 (MQ=255)
                                                                                                                                    aTTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtttt                 >  2:879816/1‑150 (MQ=255)
                                                                                                                                    aTTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtttt                 <  1:673893/150‑1 (MQ=255)
                                                                                                                                     ttCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGttttt                >  1:126685/1‑150 (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)
                                                                                                                                          aaTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAt           >  1:1456077/1‑150 (MQ=255)
                                                                                                                                             cAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAtttt        <  1:812963/150‑1 (MQ=255)
                                                                                                                                                aTGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCt     <  2:489843/150‑1 (MQ=255)
                                                                                                                                                 tgattgatTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCtt    >  2:11687/1‑150 (MQ=255)
                                                                                                                                                   attgattAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCtttt  >  2:1309543/1‑150 (MQ=255)
                                                                                                                                                   attgattAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCtttt  >  2:987566/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCTTTT  >  NC_000913/830522‑830818

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