Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 830,666 A→T 5.4% 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 5.4% 424.8 / 8.9 130intergenic (‑23/‑206)cecR/rhlEDNA‑binding transcriptional dual regulator CecR/ATP‑dependent RNA helicase RhlE
Reads supporting (aligned to +/- strand):  ref base A (64/59);  new base T (4/3);  total (68/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 131 total aligned reads displayed.

CCCGGCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCC  >  NC_000913/830518‑830814
                                                                                                                                                    |                                                                                                                                                    
cccGGCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAt                                                                                                                                                     >  1:965861/1‑150 (MQ=255)
    gCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                 <  1:1038509/150‑1 (MQ=255)
    gCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                 <  1:37531/150‑1 (MQ=255)
    gCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                 <  2:1458932/150‑1 (MQ=255)
     ccTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgattg                                                                                                                                                <  1:862686/150‑1 (MQ=255)
     ccTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgattg                                                                                                                                                <  2:1310018/150‑1 (MQ=255)
        gggCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATgattgatt                                                                                                                                             >  1:1602960/1‑150 (MQ=255)
         ggcggcTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTa                                                                                                                                            <  2:305586/150‑1 (MQ=255)
              cTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaat                                                                                                                                       <  2:1494149/150‑3 (MQ=255)
               tATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaatg                                                                                                                                      <  2:716799/150‑4 (MQ=255)
                  ctcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTAATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGa                                                                                                                                   >  2:685772/1‑150 (MQ=255)
                  ctcGCGACTGCTGCAGATCATTCCATATTCAACAAACTGCGCCAGAGCGGCAGCAATCAGCTGTTTTTTCGCATGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGa                                                                                                                                   >  2:230000/1‑150 (MQ=255)
                   tcGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGAc                                                                                                                                  >  1:507037/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)
                            gtggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAggg                                                                                                                         <  1:1518137/150‑1 (MQ=255)
                                cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                     <  2:1338815/150‑1 (MQ=255)
                                cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                     >  1:1055955/1‑150 (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:1008798/150‑1 (MQ=255)
                                       tCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAg                                                                                                              <  1:808715/150‑1 (MQ=255)
                                        ccATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGt                                                                                                             <  1:181532/150‑1 (MQ=255)
                                        ccATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGt                                                                                                             >  2:331279/1‑150 (MQ=255)
                                             ttCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGc                                                                                                        >  1:327259/1‑150 (MQ=255)
                                              tCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCt                                                                                                       >  2:86410/1‑150 (MQ=255)
                                               cACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTa                                                                                                      >  1:2343/1‑150 (MQ=255)
                                               cACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTa                                                                                                      >  2:1140787/1‑150 (MQ=255)
                                                aCCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTAt                                                                                                           >  2:1083934/1‑144 (MQ=255)
                                                aCCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTaa                                                                                                     >  1:541358/1‑150 (MQ=255)
                                                 ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAgggg                                                                                                                        <  1:120653/130‑1 (MQ=255)
                                                 ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAg                                                                                                    <  1:849129/150‑1 (MQ=255)
                                                 ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAg                                                                                                    >  2:553902/1‑150 (MQ=255)
                                                 ccAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAg                                                                                                    <  2:1147716/150‑1 (MQ=255)
                                                       tGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAAtttt                                                                                              <  1:1433824/150‑1 (MQ=255)
                                                        gcgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTa                                                                                             <  1:1140787/150‑1 (MQ=255)
                                                        gcgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTa                                                                                             <  2:1465104/150‑1 (MQ=255)
                                                         cgcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAg                                                                                            >  2:96324/1‑150 (MQ=255)
                                                          gcCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGc                                                                                           <  2:1199817/150‑1 (MQ=255)
                                                                  ggCAGCAATCAGATGTATTTTCGACTGTTCACCCTAGATTGTCATGGCAGGATTATTCATCGCACAGTCCATACTTAATCAAATCATAGATTAAGATTATGACTCCAAAGGGGAGTTGTCAAGTAAGCCTAAGAATTTAAGCAACGCCAg                                                                                   >  2:1596913/1‑150 (MQ=255)
                                                                    cagcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTc                                                                                 <  1:553902/150‑1 (MQ=255)
                                                                     agcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTca                                                                                >  1:1475611/1‑150 (MQ=255)
                                                                      gcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTcac                                                                               <  2:916204/150‑1 (MQ=255)
                                                                        aaTCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAg                                                                             <  2:387568/150‑1 (MQ=255)
                                                                        aaTCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAg                                                                             <  1:1493532/150‑1 (MQ=255)
                                                                        aaTCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAg                                                                             <  2:563567/150‑1 (MQ=255)
                                                                           cAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGa                                                                          >  2:567616/1‑150 (MQ=255)
                                                                            aGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGAt                                                                         >  1:123221/1‑150 (MQ=255)
                                                                              cTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCCGTCACAGGGATaa                                                                       >  2:1234833/1‑150 (MQ=255)
                                                                              cTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATaa                                                                       >  2:642862/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)
                                                                                    tttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATg                                                                 <  1:256634/150‑1 (MQ=255)
                                                                                     ttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATgc                                                                <  1:843781/150‑1 (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)
                                                                                                        ttGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGc                                             >  1:515080/1‑150 (MQ=255)
                                                                                                            cATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTaca                                         >  2:165758/1‑150 (MQ=255)
                                                                                                               ggCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGAGCTGCGTCACAAAAACTGCTACACTc                                      >  2:1404865/1‑150 (MQ=255)
                                                                                                                gCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                     >  2:664222/1‑150 (MQ=255)
                                                                                                                   ggATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCt                                  >  2:843793/1‑150 (MQ=255)
                                                                                                                    gATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTc                                 <  1:642862/150‑1 (MQ=255)
                                                                                                                    gATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTc                                 <  1:270533/150‑1 (MQ=255)
                                                                                                                    gATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAAGTTATGCGCTGCGTCACAAACACTGCTACACTCCGCTc                                 >  2:1591727/1‑150 (MQ=255)
                                                                                                                     attattCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                     >  2:1309049/1‑145 (MQ=255)
                                                                                                                      ttattCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTccc                               >  2:1001038/1‑150 (MQ=255)
                                                                                                                        attCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTc                             <  2:123221/150‑1 (MQ=255)
                                                                                                                        attCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTc                             <  1:12170/150‑1 (MQ=255)
                                                                                                                         ttCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCa                            >  1:278279/1‑150 (MQ=255)
                                                                                                                           cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCACCGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAACAACTGCTACACTCCGCTCCCTCATg                          >  2:1355621/1‑150 (MQ=255)
                                                                                                                           cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATg                          >  2:673764/1‑150 (MQ=255)
                                                                                                                               gCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACAt                      >  2:242109/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:1363734/150‑1 (MQ=255)
                                                                                                                                    gTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTgg                 <  1:1643921/150‑1 (MQ=255)
                                                                                                                                    gTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTgg                 <  1:567616/150‑1 (MQ=255)
                                                                                                                                     cataatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGt                >  2:1191805/6‑150 (MQ=255)
                                                                                                                                       cATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGATGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGttt              >  2:1628651/1‑150 (MQ=255)
                                                                                                                                        aatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtttt             >  2:1289170/3‑150 (MQ=255)
                                                                                                                                         ttCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGAGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGttttt            >  2:777729/1‑150 (MQ=255)
                                                                                                                                          tCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTg           >  2:1620182/1‑150 (MQ=255)
                                                                                                                                          tCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTg           <  1:682658/150‑1 (MQ=255)
                                                                                                                                            ttAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTc         <  2:1475611/150‑1 (MQ=255)
                                                                                                                                            ttAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTc         <  2:1055955/150‑1 (MQ=255)
                                                                                                                                            ttAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGCATAATTTATGCGCTGCGTCACAAAACCTGCTACACTCCGCTCCATCATGACATTGTGGTTTTTGTc         >  2:1045342/1‑150 (MQ=255)
                                                                                                                                              aaTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAt       <  1:340022/150‑1 (MQ=255)
                                                                                                                                              aaTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAt       <  1:333960/150‑1 (MQ=255)
                                                                                                                                               aTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAtt      <  1:1642045/150‑1 (MQ=255)
                                                                                                                                                   aaTGATTGATTAAGATTATGACTCCATATGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTcc  <  1:1185579/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCGGCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCC  >  NC_000913/830518‑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: ‑