Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,719,803 A→T 5.5% intergenic (+35/‑246) yiaG → / → cspA HTH_CROC1 family putative transcriptional regulator/RNA chaperone and antiterminator, cold‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,719,8030AT 5.5% 465.9 / 9.0 156intergenic (+35/‑246)yiaG/cspAHTH_CROC1 family putative transcriptional regulator/RNA chaperone and antiterminator, cold‑inducible
Reads supporting (aligned to +/- strand):  ref base A (67/80);  new base T (5/4);  total (72/84)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

AGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGT  >  NC_000913/3719654‑3719951
                                                                                                                                                     |                                                                                                                                                    
aGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCtga                                                                                                                                                      >  2:1530894/1‑150 (MQ=255)
aGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCtga                                                                                                                                                      >  1:501423/1‑150 (MQ=255)
  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCtgatg                                                                                                                                                    <  1:260261/150‑1 (MQ=255)
  tCGCCACGGGAAAGGATTGGTAATCCAGACTCGTGAAGCCTTCCAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCtgatg                                                                                                                                                    <  2:356263/150‑1 (MQ=255)
  gcGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTTAGTAAGGAGTTGATGGAATAGACTTTTATCCACTTTATTGATGTTTACGGTCCtgatg                                                                                                                                                    <  1:55424/149‑1 (MQ=255)
   cGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCtgatga                                                                                                                                                   >  1:1439601/1‑150 (MQ=255)
         tgTAAAGGAATGTGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAATATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                                                                             <  1:653657/149‑1 (MQ=255)
         ggTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                                                                             >  1:1759011/1‑150 (MQ=255)
          gTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAcc                                                                                                                                            >  2:981006/1‑150 (MQ=255)
           tAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCg                                                                                                                                           >  2:1478178/1‑150 (MQ=255)
            aaaGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGTCATCAGGACCGt                                                                                                                                          <  1:974307/150‑1 (MQ=255)
               ggAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGtttt                                                                                                                                       >  1:1036821/1‑150 (MQ=255)
               ggAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGtttt                                                                                                                                       >  2:519101/1‑150 (MQ=255)
                gAATGGGAATCCAGACGCGTTAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGATTAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTc                                                                                                                                      <  1:1113049/150‑1 (MQ=255)
                 aaTGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGTCATCAGGACCGtaaaca                                                                                                                                     <  2:436017/150‑6 (MQ=255)
                 aaTGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTcc                                                                                                                                     <  2:175119/150‑1 (MQ=255)
                 aaTGGGAATCCAGACGCGGGAAGCCTTGAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCAGTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTcc                                                                                                                                     >  2:1061704/1‑150 (MQ=255)
                  aTGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATGAAGTAAGCAGGTGATGGAATAGACTTTTATCGACTTTATTGGTGTTTACGGTCGTGATGACAGGACCGTTTTcgc                                                                                                                                    >  2:1935421/1‑148 (MQ=255)
                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAATATTGCTGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAc                                                                                                                                  <  2:530799/150‑1 (MQ=255)
                     ggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAcc                                                                                                                                 >  1:1360484/1‑150 (MQ=255)
                      gAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCg                                                                                                                                >  1:356751/1‑150 (MQ=255)
                           cAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTaa                                                                                                                           <  1:798235/150‑1 (MQ=255)
                             gACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATc                                                                                                                         >  2:919217/1‑150 (MQ=255)
                             gACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATc                                                                                                                         >  1:626888/1‑150 (MQ=255)
                             gACGCGTGAACCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATc                                                                                                                         <  1:1259145/150‑1 (MQ=255)
                             gACGAGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGACTAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTACTc                                                                                                                         <  2:330337/150‑1 (MQ=255)
                                   tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAat                                                                                                                   <  2:1708123/150‑1 (MQ=255)
                                    gAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAata                                                                                                                  >  2:811800/1‑150 (MQ=255)
                                    gAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCAATAAGTAAGCAGTTCATGGAATAGACTTTTATCCAGTATATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAata                                                                                                                  >  2:399441/1‑150 (MQ=255)
                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAaat                                                                                                          <  1:1311753/150‑1 (MQ=255)
                                              aGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAaataa                                                                                                        >  2:772739/1‑150 (MQ=255)
                                               gTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAaataat                                                                                                       >  1:132212/1‑150 (MQ=255)
                                                  ccGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttgt                                                                                                    <  2:1512442/150‑1 (MQ=255)
                                                    gAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttgttg                                                                                                  <  1:1188411/150‑1 (MQ=255)
                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGc                                                                                                 <  2:129276/150‑1 (MQ=255)
                                                      aCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCa                                                                                                <  1:1768353/150‑1 (MQ=255)
                                                      aCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCa                                                                                                <  1:607373/150‑1 (MQ=255)
                                                       cTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATTGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCAt                                                                                               <  1:772726/150‑1 (MQ=255)
                                                          aaaTTGATGCGGTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAc                                                                                            <  2:1381430/150‑1 (MQ=255)
                                                             ttGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAcctg                                                                                         <  2:1642146/150‑3 (MQ=255)
                                                               gATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGTACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGcc                                                                                       <  1:993554/150‑1 (MQ=255)
                                                                  gCGTTTGCTTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAAt                                                                                    >  1:380690/1‑150 (MQ=255)
                                                                    gTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGCCTTTTTTCCACTTTATTGCTGTTTACTGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATTAAAAATAATTTTTTCATCACCCGCCAATGc                                                                                  >  1:71704/1‑150 (MQ=255)
                                                                       tGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTg                                                                               >  2:1102197/1‑150 (MQ=255)
                                                                        gATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTgg                                                                              <  1:2046264/150‑1 (MQ=255)
                                                                          ttCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAATAATAATTGTTGCATCACCCGCCAATGCGTGGCt                                                                            >  1:290624/1‑150 (MQ=255)
                                                                          ttCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCt                                                                            >  2:1131433/1‑150 (MQ=255)
                                                                           tCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCtt                                                                           <  2:1364590/150‑1 (MQ=255)
                                                                            cAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGCCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTa                                                                          >  2:52727/1‑150 (MQ=255)
                                                                              aGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAAt                                                                        >  1:1381695/1‑150 (MQ=255)
                                                                                ccAACCCGGCATTAAGTAAGCAGTTGATTGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGc                                                                      <  1:1987609/150‑1 (MQ=255)
                                                                                ccAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGc                                                                      <  1:285628/150‑1 (MQ=255)
                                                                                ccAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGc                                                                      <  2:1400875/150‑1 (MQ=255)
                                                                                  aaCCCGGCATTAAGTAAGCAGTTGATTGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCAAATGCGTGGCTTAATGcac                                                                    <  2:1009283/150‑1 (MQ=255)
                                                                                   aCCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGcaca                                                                   <  1:1633726/150‑1 (MQ=255)
                                                                                    cccGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGcacaa                                                                  >  1:592491/1‑149 (MQ=255)
                                                                                    cccGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACAt                                                                  >  1:850224/1‑150 (MQ=255)
                                                                                      cGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCa                                                                <  1:1202988/150‑1 (MQ=255)
                                                                                        gCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAAc                                                              <  2:1796723/150‑1 (MQ=255)
                                                                                           ttAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGt                                                           >  1:1955274/1‑150 (MQ=255)
                                                                                               gtaagCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGa                                                       <  2:779328/150‑1 (MQ=255)
                                                                                                 aagCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGCCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCCTCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGCCg                                                     >  2:847485/1‑150 (MQ=255)
                                                                                                 aagCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACg                                                     >  1:1302369/1‑150 (MQ=255)
                                                                                                    cAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATTAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTAc                                                  >  1:170031/1‑150 (MQ=255)
                                                                                                    cAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCATTGCGTGGCTTAATGCACATCAACGGTTTGACGTAc                                                  >  1:1113108/1‑150 (MQ=255)
                                                                                                     aGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACa                                                 >  1:1402632/1‑150 (MQ=255)
                                                                                                     aGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACa                                                 <  2:187096/150‑1 (MQ=255)
                                                                                                         gATTGAATAGACTTTTATCCACTTTATTGCTGTTTCAGGTCCTGATGACAGGACCGTTTTACAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCATCGGTTTGACGTACAGAcc                                             <  2:169439/150‑1 (MQ=255)
                                                                                                         gATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGAcc                                             <  1:102143/150‑1 (MQ=255)
                                                                                                          aTGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCa                                            <  2:874890/150‑1 (MQ=255)
                                                                                                          aTGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCa                                            <  2:1212825/150‑1 (MQ=255)
                                                                                                           tGGAATAGACTTTTATCAACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCAt                                           <  2:1225546/150‑1 (MQ=255)
                                                                                                            ggAATAGACTTTTATCCACTTTATTGCTGTTTACTGTCCTGATGACAGGACCGTTTTCCAACCGATTAATAATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGTCTTAATTCAAATCAAAGGTTTGACGTAAAGACCAtt                                          >  2:405867/1‑150 (MQ=255)
                                                                                                            ggAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCAGCAAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCAtt                                          >  1:1528974/1‑150 (MQ=255)
                                                                                                             gAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCTATTAATCATAATTATGAAAAATTATTGTTGCATCACCCGACAATGCGTGGCTTTATTCACATAAACGGTTTGACGTACATACCAtga                                         >  1:1776597/1‑148 (MQ=255)
                                                                                                              aaTAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATAACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTaa                                        <  2:338274/150‑1 (MQ=255)
                                                                                                               aTAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTaaa                                       >  2:1680066/1‑150 (MQ=255)
                                                                                                                tAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAg                                      >  1:461577/1‑150 (MQ=255)
                                                                                                                      tttATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAgtgt                                <  1:873007/150‑1 (MQ=255)
                                                                                                                       ttATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGATCGTTTTCCAACCGATTTATCTTAAATATGAAAAATATTTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTgta                               <  1:1478179/150‑1 (MQ=255)
                                                                                                                        tATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTgtag                              <  1:187396/150‑1 (MQ=255)
                                                                                                                            cACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAATTATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAg                          >  1:2010020/1‑150 (MQ=255)
                                                                                                                            cACTGTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCACCCGAGTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGCCGGACAGACAATTAACGCAGTggagtacg                          >  2:1729185/1‑148 (MQ=255)
                                                                                                                               tttATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCa                       <  2:1360483/150‑1 (MQ=255)
                                                                                                                               tttATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCa                       <  2:838835/150‑1 (MQ=255)
                                                                                                                                    tGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTcc                  <  2:597165/150‑1 (MQ=255)
                                                                                                                                     gCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTccc                 >  2:1150483/1‑150 (MQ=255)
                                                                                                                                     gCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTccc                 >  2:1162689/1‑150 (MQ=255)
                                                                                                                                       ggaaaaCGGTCCTGTCATCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCtt               >  1:770950/6‑150 (MQ=255)
                                                                                                                                        gaaaaCGGTCCTGTCATCAGGCCCGTTTTCCAACCGATTAATCATAAATATGATAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTAGACGTACAGACCATTAAAGCAGTGTAGTCAGGCAAGGCCCttg              >  1:307468/5‑149 (MQ=255)
                                                                                                                                         tttACGGTCCTGATGACAGGACCGTTTTGCAACCGATTAATCATAACTATGAAAAATAATTGTTGCATCACCAGCCAATGCGTGGCTGAATGCACATCAACGGTTTTACGTACAGACCATTAACGCAGTGTAGTAAGGCAAGTCCCTTCa             <  2:1568670/150‑1 (MQ=255)
                                                                                                                                         tttACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCa             <  1:519101/150‑1 (MQ=255)
                                                                                                                                          aaaCGGTCCTGTCATCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCAGCCAATGCGTGGCTTAATGCACATAAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCaa            >  2:897980/3‑150 (MQ=255)
                                                                                                                                           aaCGGTCCTGTCATCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAag           >  1:231447/2‑150 (MQ=255)
                                                                                                                                              ggTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGt        <  2:170031/150‑1 (MQ=255)
                                                                                                                                                tCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTa      <  1:984186/150‑1 (MQ=255)
                                                                                                                                                 ccTGATTACAGGAACGTTTTCGAACCGATTAATCATAAATATTAATATTAATTGTTGCATCACACGCCTCTGCGTGGCTTAATGCTCATCAACGGTTTGACGTAAAGACCATTAAAGCAGTGTAGTAATGCAAGTCCCTTAAAGAGTTAt     <  2:2044639/150‑1 (MQ=255)
                                                                                                                                                  cTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATc    <  1:1815427/150‑1 (MQ=255)
                                                                                                                                                  cTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATc    >  1:1830055/1‑150 (MQ=255)
                                                                                                                                                    gatgaCAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGt  <  2:660258/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACGTACAGACCATTAAAGCAGTGTAGTAAGGCAAGTCCCTTCAAGAGTTATCGT  >  NC_000913/3719654‑3719951

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