Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 830,668 G→A 5.7% intergenic (‑25/‑204) 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,6680GA 5.7% 382.5 / 7.9 122intergenic (‑25/‑204)cecR/rhlEDNA‑binding transcriptional dual regulator CecR/ATP‑dependent RNA helicase RhlE
Reads supporting (aligned to +/- strand):  ref base G (60/53);  new base A (4/3);  total (66/56)
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.36e-01

Only 100 of 129 total aligned reads displayed.

GCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCC  >  NC_000913/830522‑830814
                                                                                                                                                  |                                                                                                                                                  
gCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                 <  1:1038509/150‑1 (MQ=255)
gCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                 <  1:37531/150‑1 (MQ=255)
 ccTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgattg                                                                                                                                                <  2:1310018/150‑1 (MQ=255)
 ccTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgattg                                                                                                                                                <  1:862686/150‑1 (MQ=255)
    gggCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATgattgatt                                                                                                                                             >  1:1602960/1‑150 (MQ=255)
    gggCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATgattgatt                                                                                                                                             <  2:600097/150‑1 (MQ=255)
     ggcggcTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTa                                                                                                                                            <  2:305586/150‑1 (MQ=255)
          cTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaat                                                                                                                                       <  2:1494149/150‑3 (MQ=255)
           tATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaatg                                                                                                                                      <  2:716799/150‑4 (MQ=255)
           tATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTa                                                                                                                                      >  1:563567/1‑150 (MQ=255)
           tATCTCGCGAGTGGTGGCGTTCATTACATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGATCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaatg                                                                                                                                      <  1:1489217/150‑4 (MQ=255)
              ctcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTAATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGa                                                                                                                                   >  2:685772/1‑150 (MQ=255)
              ctcGCGACTGCTGCAGATCATTCCATATTCAACAAACTGCGCCAGAGCGGCAGCAATCAGCTGTTTTTTCGCATGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGa                                                                                                                                   >  2:230000/1‑150 (MQ=255)
               tcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGAc                                                                                                                                  >  1:1199817/1‑150 (MQ=255)
                cgcgAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACt                                                                                                                                 >  1:1442910/1‑150 (MQ=255)
                cgcgAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACt                                                                                                                                 >  1:1100583/1‑150 (MQ=255)
                       ggtggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAgg                                                                                                                          >  2:1581925/1‑150 (MQ=255)
                       ggtggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAgg                                                                                                                          >  1:359104/1‑150 (MQ=255)
                        gtggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAggg                                                                                                                         <  1:1518137/150‑1 (MQ=255)
                            cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                     <  2:1442910/150‑1 (MQ=255)
                            cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                     >  1:1055955/1‑150 (MQ=255)
                            cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                     <  2:1338815/150‑1 (MQ=255)
                            cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                     >  2:843781/1‑150 (MQ=255)
                             gTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGACTGATTAAGATTATGACTCCATAGGGGAGtt                                                                                                                    >  2:1191951/1‑150 (MQ=255)
                                cATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTc                                                                                                                 >  2:56962/1‑150 (MQ=255)
                                   tCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAg                                                                                                              <  1:808715/150‑1 (MQ=255)
                                   tCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAg                                                                                                              <  1:1008798/150‑1 (MQ=255)
                                    ccATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGt                                                                                                             <  1:181532/150‑1 (MQ=255)
                                    ccATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGt                                                                                                             >  2:331279/1‑150 (MQ=255)
                                          tCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCt                                                                                                       >  2:86410/1‑150 (MQ=255)
                                           cACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTa                                                                                                      >  2:1140787/1‑150 (MQ=255)
                                           cACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTa                                                                                                      >  1:2343/1‑150 (MQ=255)
                                            aCCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTAt                                                                                                           >  2:1083934/1‑144 (MQ=255)
                                             ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAgggg                                                                                                                        >  2:120653/1‑130 (MQ=255)
                                             ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAgggg                                                                                                                        <  1:120653/130‑1 (MQ=255)
                                             ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAg                                                                                                    >  2:553902/1‑150 (MQ=255)
                                             ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAg                                                                                                    <  1:849129/150‑1 (MQ=255)
                                             ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAg                                                                                                    <  2:1147716/150‑1 (MQ=255)
                                                   tGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAAtttt                                                                                              <  1:1433824/150‑1 (MQ=255)
                                                    gcgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTa                                                                                             <  2:1465104/150‑1 (MQ=255)
                                                    gcgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTa                                                                                             <  1:1140787/150‑1 (MQ=255)
                                                      gcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGc                                                                                           <  2:1199817/150‑1 (MQ=255)
                                                              ggCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAg                                                                                   >  2:860858/1‑150 (MQ=255)
                                                               gcagcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGt                                                                                  >  2:1642045/1‑150 (MQ=255)
                                                                cagcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTc                                                                                 <  1:553902/150‑1 (MQ=255)
                                                                 agcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTca                                                                                >  1:1475611/1‑150 (MQ=255)
                                                                  gcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTcac                                                                               <  2:201001/150‑1 (MQ=255)
                                                                  gcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTcac                                                                               <  2:916204/150‑1 (MQ=255)
                                                                    aaTCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAg                                                                             <  2:563567/150‑1 (MQ=255)
                                                                        aGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGAt                                                                         >  1:123221/1‑150 (MQ=255)
                                                                          cTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCCGTCACAGGGATaa                                                                       >  2:1234833/1‑150 (MQ=255)
                                                                          cTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATaa                                                                       >  2:642862/1‑150 (MQ=255)
                                                                           tGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAAt                                                                      >  1:1614598/1‑150 (MQ=255)
                                                                            gTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAAtt                                                                     <  2:1669141/150‑1 (MQ=255)
                                                                             ttttttCGCCTGTTCACCCTTGATTGTGTTGGCAGGATTATTGATCTCACATTCCATTCTTAATCAAATGATTGATTAAGATTATGACTGCATAGGGGATTTGTCCAGTATGTGTAAGTATTTTAGGAACGCCAGTCACAGGGATAAttt                                                                    <  2:114631/150‑1 (MQ=255)
                                                                              tttttCGCCTATTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTa                                                                   <  1:750161/150‑1 (MQ=255)
                                                                               ttttCGCCTGTTCACCCTTGATGGTCATGGCAGGATTATTAATCGCACAGTCCATTCTTAATCAAATGATTGAATAACATTATGCCTCCATAGGGGAATTTTCCAGTATGCCTAATAATTTTAGCAAAGGCACTCACAGCGATAATTTAt                                                                  >  2:1463274/1‑150 (MQ=255)
                                                                                tttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATg                                                                 >  2:124625/1‑150 (MQ=255)
                                                                                 ttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATgc                                                                <  1:1208459/150‑1 (MQ=255)
                                                                                       tGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCg                                                          >  2:1041660/1‑150 (MQ=255)
                                                                                        gTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGt                                                         >  2:406312/1‑150 (MQ=255)
                                                                                         ttCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTc                                                        <  1:1261216/150‑1 (MQ=255)
                                                                                         ttCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCAGCGTc                                                        >  2:1666989/1‑150 (MQ=255)
                                                                                          tCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTca                                                       <  2:78892/150‑1 (MQ=255)
                                                                                           cACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTcac                                                      <  1:537487/150‑1 (MQ=255)
                                                                                             cccTTGATTGTGATGGCAGGATTATTAATCGCACAGTCCATTATTTATAAAATTATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACaa                                                    >  2:44647/1‑150 (MQ=255)
                                                                                             cccTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACaa                                                    >  2:788783/1‑150 (MQ=255)
                                                                                                 tGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAAc                                                <  1:230000/150‑1 (MQ=255)
                                                                                                    ttGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGc                                             >  1:515080/1‑150 (MQ=255)
                                                                                                      gTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTa                                           >  2:1470114/1‑150 (MQ=255)
                                                                                                        cATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTaca                                         >  1:277653/1‑150 (MQ=255)
                                                                                                        cATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTaca                                         >  2:682658/1‑150 (MQ=255)
                                                                                                           ggCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGAGCTGCGTCACAAAAACTGCTACACTc                                      >  2:1404865/1‑150 (MQ=255)
                                                                                                                gATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTc                                 <  1:270533/150‑1 (MQ=255)
                                                                                                                gATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAAGTTATGCGCTGCGTCACAAACACTGCTACACTCCGCTc                                 >  2:1591727/1‑150 (MQ=255)
                                                                                                                 attattCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                     >  2:1309049/1‑145 (MQ=255)
                                                                                                                 attattCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                     <  1:1309049/145‑1 (MQ=255)
                                                                                                                  ttattCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTccc                               <  1:745503/150‑1 (MQ=255)
                                                                                                                    attCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTc                             <  1:519044/150‑1 (MQ=255)
                                                                                                                    attCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTc                             <  2:327259/150‑1 (MQ=255)
                                                                                                                     ttCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCa                            >  1:278279/1‑150 (MQ=255)
                                                                                                                       cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATg                          >  2:673764/1‑150 (MQ=255)
                                                                                                                        aTCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGa                         >  2:1343964/1‑150 (MQ=255)
                                                                                                                           gCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACAt                      >  2:1405477/1‑150 (MQ=255)
                                                                                                                             acaGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATtg                    <  1:1041660/150‑1 (MQ=255)
                                                                                                                               gtcataatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATtgtg                  >  1:1125835/8‑150 (MQ=255)
                                                                                                                                gTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTgg                 <  1:1643921/150‑1 (MQ=255)
                                                                                                                                gTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTgg                 <  1:567616/150‑1 (MQ=255)
                                                                                                                                 cataatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGt                >  2:1191805/6‑150 (MQ=255)
                                                                                                                                 cataatCTTAATCAATCATTTGATTAAGATTATGACTCCATAAGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGt                >  2:1589197/6‑150 (MQ=255)
                                                                                                                                   cATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGATGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGttt              >  2:1628651/1‑150 (MQ=255)
                                                                                                                                    aatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtttt             >  2:1289170/3‑150 (MQ=255)
                                                                                                                                      tCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTg           >  2:1620182/1‑150 (MQ=255)
                                                                                                                                      tCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTg           <  1:682658/150‑1 (MQ=255)
                                                                                                                                        ttAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTc         <  2:1055955/150‑1 (MQ=255)
                                                                                                                                        ttAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGCATAATTTATGCGCTGCGTCACAAAACCTGCTACACTCCGCTCCATCATGACATTGTGGTTTTTGTc         >  2:1045342/1‑150 (MQ=255)
                                                                                                                                          aaTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAt       <  1:333960/150‑1 (MQ=255)
                                                                                                                                          aaTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAt       <  1:340022/150‑1 (MQ=255)
                                                                                                                                          aaTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAt       <  1:843793/150‑1 (MQ=255)
                                                                                                                                               aaTGATTGATTAAGATTATGACTCCATATGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTcc  <  1:1185579/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCC  >  NC_000913/830522‑830814

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: ‑