Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 830,669 A→T 5.6% 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 5.6% 396.9 / 7.9 122intergenic (‑26/‑203)cecR/rhlEDNA‑binding transcriptional dual regulator CecR/ATP‑dependent RNA helicase RhlE
Reads supporting (aligned to +/- strand):  ref base A (62/53);  new base T (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.75e-01

Only 100 of 130 total aligned reads displayed.

GCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCTTTT  >  NC_000913/830522‑830818
                                                                                                                                                   |                                                                                                                                                     
gCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                     <  1:1038509/150‑1 (MQ=255)
gCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                     <  1:37531/150‑1 (MQ=255)
gCCTGGGCGGCTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAAtgatt                                                                                                                                                     <  2:1458932/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)
         gcTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGAt                                                                                                                                            >  2:1185579/1‑150 (MQ=255)
          cTATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaat                                                                                                                                           <  2:1494149/150‑3 (MQ=255)
           tATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAATCATTTGATTAAGaatg                                                                                                                                          <  2:716799/150‑4 (MQ=255)
           tATCTCGCGAGTGGTGGCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTa                                                                                                                                          >  1:563567/1‑150 (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)
                        gtggCGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAggg                                                                                                                             <  1:1518137/150‑1 (MQ=255)
                            cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                         >  1:1055955/1‑150 (MQ=255)
                            cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                         >  2:843781/1‑150 (MQ=255)
                            cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                         <  2:1442910/150‑1 (MQ=255)
                            cGTTCATTCCATATTCACCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGt                                                                                                                         <  2:1338815/150‑1 (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                                                                                                          >  2:1140787/1‑150 (MQ=255)
                                            aCCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTAt                                                                                                               >  2:1083934/1‑144 (MQ=255)
                                            aCCAAACTGCGCCAGTGCGGCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTaa                                                                                                         >  1:541358/1‑150 (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                                                                                                 <  1:1140787/150‑1 (MQ=255)
                                                              ggCAGCAATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAg                                                                                       >  2:860858/1‑150 (MQ=255)
                                                              ggCAGCAATCAGATGTATTTTCGACTGTTCACCCTAGATTGTCATGGCAGGATTATTCATCGCACAGTCCATACTTAATCAAATCATAGATTAAGATTATGACTCCAAAGGGGAGTTGTCAAGTAAGCCTAAGAATTTAAGCAACGCCAg                                                                                       >  2:1596913/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:916204/150‑1 (MQ=255)
                                                                  gcaATCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTcac                                                                                   <  2:201001/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)
                                                                      tCAGCTGTTTTTTCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAggg                                                                               <  2:1100583/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)
                                                                               ttttCGCCTGTTCACCCTTGATGGTCATGGCAGGATTATTAATCGCACAGTCCATTCTTAATCAAATGATTGAATAACATTATGCCTCCATAGGGGAATTTTCCAGTATGCCTAATAATTTTAGCAAAGGCACTCACAGCGATAATTTAt                                                                      >  2:1463274/1‑150 (MQ=255)
                                                                                tttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATg                                                                     <  1:256634/150‑1 (MQ=255)
                                                                                 ttCGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATgc                                                                    <  1:843781/150‑1 (MQ=255)
                                                                                   cGCCTGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTTAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATgcgc                                                                  >  2:1302330/1‑150 (MQ=255)
                                                                                       tGTTCACCCTTGATTGTCATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCg                                                              >  2:1041660/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)
                                                                                             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                                             >  2:682658/1‑150 (MQ=255)
                                                                                                        cATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTaca                                             >  2:165758/1‑150 (MQ=255)
                                                                                                        cATGGCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTaca                                             >  1:277653/1‑150 (MQ=255)
                                                                                                           ggCAGGATTATTCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGAGCTGCGTCACAAAAACTGCTACACTc                                          >  2:1404865/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)
                                                                                                                 attattCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTcc                                         <  1:1309049/145‑1 (MQ=255)
                                                                                                                  ttattCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTccc                                   >  2:1001038/1‑150 (MQ=255)
                                                                                                                  ttattCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTccc                                   <  1:745503/150‑1 (MQ=255)
                                                                                                                    attCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTc                                 <  2:123221/150‑1 (MQ=255)
                                                                                                                    attCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTc                                 <  1:12170/150‑1 (MQ=255)
                                                                                                                    attCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTc                                 <  1:519044/150‑1 (MQ=255)
                                                                                                                    attCATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTc                                 <  2:327259/150‑1 (MQ=255)
                                                                                                                       cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCACCGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAACAACTGCTACACTCCGCTCCCTCATg                              >  2:1355621/1‑150 (MQ=255)
                                                                                                                       cATCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATg                              >  2:673764/1‑150 (MQ=255)
                                                                                                                        aTCGCACAGTCCATTCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGa                             >  2:1343964/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)
                                                                                                                                 cataatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGt                    >  2:1191805/6‑150 (MQ=255)
                                                                                                                                 cataatCTTAATCAATCATTTGATTAAGATTATGACTCCATAAGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGt                    >  2:1589197/6‑150 (MQ=255)
                                                                                                                                    aatCTTAATCAATCATTTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGtttt                 >  2:1289170/3‑150 (MQ=255)
                                                                                                                                      tCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTg               <  1:682658/150‑1 (MQ=255)
                                                                                                                                      tCTTAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTg               >  2:1620182/1‑150 (MQ=255)
                                                                                                                                        ttAATCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTc             <  2:1475611/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)
                                                                                                                                           aTCAAATGATTGATTAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCAtt          <  1:1642045/150‑1 (MQ=255)
                                                                                                                                               aaTGATTGATTAAGATTATGACTCCATATGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTcc      <  1:1185579/150‑1 (MQ=255)
                                                                                                                                                   attgattAAGATTATGACTCCATAGGGGAGTTGTCCAGTATGGCTAAGAATTTTAGCAACGCCAGTCACAGGGATAATTTATGCGCTGCGTCACAAAAACTGCTACACTCCGCTCCCTCATGACATTGTGGTTTTTGTCATTTTCCtttt  >  2:191094/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 < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑