Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,217,119 C→A 5.6% intergenic (+123/‑209) ymgC → / → ycgG Blue light, low temperature and stress induced protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,217,1190CA 5.6% 289.0 / 5.5 106intergenic (+123/‑209)ymgC/ycgGBlue light, low temperature and stress induced protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (61/37);  new base A (4/2);  total (65/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 112 total aligned reads displayed.

TGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATG  >  NC_000913/1216972‑1217268
                                                                                                                                                   |                                                                                                                                                     
tGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGc                                                                                                                                                     >  2:1722044/1‑150 (MQ=255)
tGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGc                                                                                                                                                     >  1:1557830/1‑150 (MQ=255)
 gACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGct                                                                                                                                                    >  1:6481/1‑150 (MQ=255)
  aCAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATAAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGATTTCTTATTATGTCCGCAATATCATGTGTCAAGAATGGAGAGTTCTCGctc                                                                                                                                                   >  1:886089/1‑150 (MQ=255)
   cAGGGAATCCGTGCTCTCTTAGAGCTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGctct                                                                                                                                                  >  2:725300/1‑150 (MQ=255)
   cAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGctct                                                                                                                                                  <  2:527292/150‑1 (MQ=255)
    aGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCtctc                                                                                                                                                 >  1:1689056/1‑150 (MQ=255)
       gAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCAt                                                                                                                                              <  2:973596/150‑1 (MQ=255)
       gAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCAt                                                                                                                                              >  2:777013/1‑150 (MQ=255)
         aTCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGGCAAGAATGGAGAGTTCTCGCTCTCCATTc                                                                                                                                            <  1:1222212/150‑1 (MQ=255)
           ccGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGCGAGAACTCTCCATTCtt                                                                                                                                          <  2:1588797/150‑1 (MQ=39)
           ccGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACAATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCtt                                                                                                                                          >  2:1227587/1‑150 (MQ=255)
            cGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTg                                                                                                                                         <  2:1628117/150‑1 (MQ=255)
                ctctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTAATATGTCAGCAATATCAGGTGTCAATAATGGAGAGTTCTCGATCTCAATTCTTGAAGc                                                                                                                                     >  1:284634/1‑150 (MQ=255)
                 tctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGcc                                                                                                                                    <  1:1713444/150‑1 (MQ=255)
                  ctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt                                                                                                                                   >  2:951223/1‑150 (MQ=255)
                  ctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt                                                                                                                                   >  2:1202391/1‑150 (MQ=255)
                   tctTAGAACTTGAGTAAGGTATTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGACACCTg                                                                                                                                  <  1:167858/150‑1 (MQ=37)
                    ctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGa                                                                                                                                 >  2:632002/1‑150 (MQ=255)
                      tAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGata                                                                                                                               >  2:334345/1‑150 (MQ=255)
                        gAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATc                                                                                                                             >  1:1402313/1‑150 (MQ=255)
                        gAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATc                                                                                                                             >  1:1709551/1‑150 (MQ=255)
                          aCTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATccc                                                                                                                           >  2:6082/1‑150 (MQ=255)
                              gAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCt                                                                                                                       >  2:1545975/1‑150 (MQ=255)
                                gTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCAATTATTGACGCCTGATATCCAGCCTaa                                                                                                                     >  1:588085/1‑150 (MQ=255)
                                  atgtttttTGAATAAAAATCTGGTTTATCACCATGCTGAGGAATTTTGAGATTTTCAGGCTTTATTGTTTGTTTTTATGTCCGAAATATCTGGTGTCAAGTTTGGATTTTTCTGGCTCTCATTTATTGACGCCTGATATCCCGCCTAACt                                                                                                                   <  2:906699/146‑1 (MQ=255)
                                   aGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACtt                                                                                                                  >  1:793858/1‑150 (MQ=255)
                                        ttAGACTAAACATCAGCTGTATCACCATGCTCATGCAAAGTGAGCATTTCAGGCGTTAGGCTTTCTTATTATGTCCGCAATATCCGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTAtttt                                                                                                             <  1:337683/150‑1 (MQ=255)
                                        ttAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTAtttt                                                                                                             >  2:1642775/1‑150 (MQ=255)
                                         tAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTg                                                                                                            >  2:208143/1‑150 (MQ=255)
                                          ataaTAAACATAAGCTGTATCACAATGCTGATGAAAAGTGAGAATTTCTGGCTTTATGCTTTTTTATTATGTCCGAAATATCTGGTGTCAAGAATGGTGAGTTCTCGCTCTCCATTTTTGACGCCTGATATCCCGCCTAACTTATTTTGt                                                                                                           <  2:632100/148‑1 (MQ=255)
                                          aGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGt                                                                                                           >  2:1263998/1‑150 (MQ=255)
                                           gAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTa                                                                                                          >  1:1361258/1‑150 (MQ=255)
                                             aTAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACt                                                                                                        >  2:1514440/1‑150 (MQ=255)
                                               aaaCATCAGCTTTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGc                                                                                                      <  2:6893/150‑1 (MQ=255)
                                               aaaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGc                                                                                                      <  2:645412/150‑1 (MQ=255)
                                                aaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGcc                                                                                                     <  1:1168000/150‑1 (MQ=255)
                                                  cATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTa                                                                                                   >  1:1472111/1‑150 (MQ=255)
                                                   aTCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTAc                                                                                                  <  2:16908/150‑1 (MQ=255)
                                                       tcTGTATCACCATGCTTATGCAAATTGAGCATTTCAGGAGTTATGATTTATTATTTTGTCGGAAATATAAGGTGTCAAGAATGGACAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTa                                                                                              <  1:1094707/149‑1 (MQ=255)
                                                       gCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGAGATCCCGCCGAACTTAGTTTGTACTGCCTACAttt                                                                                              >  2:722692/1‑149 (MQ=255)
                                                           tATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTTTGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGTAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCt                                                                                          <  1:1007475/150‑1 (MQ=255)
                                                             tCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCACAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCa                                                                                        >  1:651556/1‑150 (MQ=255)
                                                              cACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTATCGCTCTACATTCTTGACTCCTGATATCCCGCATTACTTATTTTGTACTGCATACATTATGCTCaa                                                                                       >  1:75040/1‑150 (MQ=255)
                                                               aCCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATACCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaat                                                                                      >  2:872253/1‑150 (MQ=255)
                                                                ccATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaata                                                                                     >  2:1671083/1‑150 (MQ=255)
                                                                   tGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                                  >  1:664103/1‑150 (MQ=255)
                                                                     cTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATAACCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaattaa                                                                                >  2:329405/1‑150 (MQ=255)
                                                                       gATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATg                                                                              >  1:95402/1‑150 (MQ=255)
                                                                         tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                            >  1:1813952/1‑150 (MQ=255)
                                                                         tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                            >  2:1285600/1‑150 (MQ=255)
                                                                                 gAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCa                                                                    <  1:194321/150‑1 (MQ=255)
                                                                                 gAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCa                                                                    >  2:1350351/1‑150 (MQ=255)
                                                                                  aGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCaa                                                                   >  1:726892/1‑150 (MQ=255)
                                                                                   gCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAg                                                                  <  2:60959/150‑1 (MQ=255)
                                                                                   gCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGATCAATAATTAATGGGTTAGATCAAg                                                                  >  2:473235/1‑150 (MQ=255)
                                                                                    cATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTAAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGa                                                                 <  1:891707/150‑1 (MQ=255)
                                                                                     aTTTCAGGCGTTATGCTTTCTTTTTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTAGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCATACATTATGCTCAAAAATTAATGGGTTAGATCAAGAt                                                                >  1:1201036/1‑150 (MQ=255)
                                                                                      tttCAGGCTTTATGATTTCTAATTATGTCCTAAATATCAGGTGTCAATATTGGATATTTCTCGATCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATg                                                               <  1:631793/150‑1 (MQ=255)
                                                                                         cagggtttATGATTTATTATTATGTCCACAATATCAGGTTTAAATAATGGAGTTTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGgag                                                            <  1:762593/144‑1 (MQ=255)
                                                                                            gCGTTATGCTTTCTTATTATGTCCGAAATATCTGGTGTCAATATTGTTGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGaaa                                                         <  2:610153/150‑1 (MQ=255)
                                                                                               ttATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAAtt                                                      >  1:1252230/1‑150 (MQ=255)
                                                                                                tATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATTGGTTAGATCAAGATGGAGAAAAttt                                                     >  1:1620054/1‑150 (MQ=255)
                                                                                                tATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAAttt                                                     <  2:6481/150‑1 (MQ=255)
                                                                                                        cTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCg                                             <  1:1201271/150‑1 (MQ=255)
                                                                                                            ttatGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGt                                         >  1:528647/1‑150 (MQ=255)
                                                                                                             tatGTCCGAAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTg                                        <  2:456886/150‑1 (MQ=255)
                                                                                                              atGTCCGCAATATCTGGGGGCATTAATGTTGTGTTCTCGCTATCAATTCTTGTCGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACCTTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTgg                                       <  1:1242582/150‑1 (MQ=255)
                                                                                                                  ccGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGCCGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATAAAGATGGAGAAAATTTTATCTTCGGCGTGGatat                                   >  2:1674113/1‑150 (MQ=255)
                                                                                                                    gCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATAttt                                 >  1:981186/1‑150 (MQ=255)
                                                                                                                    gCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACAGCCTACATTATGCTCAATAATTAATCGGTTAGATAAATAAGGAGAAAAATTTATCTTCGGCGTGGATAttt                                 >  2:1630201/1‑150 (MQ=255)
                                                                                                                        tatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTa                             >  2:1012326/1‑150 (MQ=255)
                                                                                                                        tatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTa                             <  2:50186/150‑1 (MQ=255)
                                                                                                                         ttctggagtCAAGAATGGTGAGTTCTTGATCTCGATTCTTGACGCCTGATATCGAGCCTTACTTATTTTGTACTGCCTACATTATGCTCATTAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTAt                            <  1:54881/143‑1 (MQ=255)
                                                                                                                           cAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGc                          <  1:951223/150‑1 (MQ=255)
                                                                                                                             ggcgtCAAGAATGGAGAGCGAGAACTCTCCATTCTTGACACCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGctc                        >  2:1185292/4‑150 (MQ=37)
                                                                                                                             ggTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCGTGATTTCGGGCCTAACTTTGTTTGTACTGCCTAAATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGctc                        >  1:108304/1‑150 (MQ=255)
                                                                                                                                gtaaaGAATGGAGAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGctctct                     <  2:1769636/147‑1 (MQ=255)
                                                                                                                                gtCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGctctct                     <  1:632002/150‑1 (MQ=255)
                                                                                                                                 tCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTa                    >  2:840575/1‑150 (MQ=255)
                                                                                                                                  cAAGAATGGAGAGCGAGAACTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTTTCTTCGGCGTGGATATTTCTTATGCTCTCTAc                   >  1:1120647/1‑150 (MQ=38)
                                                                                                                                  cAAGAATGGAGAGCGAGAACTCTCCATTCTTGACACCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATTATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTAc                   >  2:805197/1‑150 (MQ=38)
                                                                                                                                   aaGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGATCTCTAcc                  >  2:354177/1‑150 (MQ=255)
                                                                                                                                     gAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCAt                >  1:1788945/1‑150 (MQ=255)
                                                                                                                                     gAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTAGGCGTGGATATTTCTTATGCTCTCTACCAt                >  1:1315660/1‑150 (MQ=255)
                                                                                                                                       aTGGATAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc              <  2:566123/150‑1 (MQ=255)
                                                                                                                                       aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc              >  1:1122920/1‑150 (MQ=255)
                                                                                                                                        tGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGct             >  1:1314819/1‑150 (MQ=255)
                                                                                                                                         ggAGAGTTCTCGTGCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctc            <  2:1135655/150‑1 (MQ=255)
                                                                                                                                          gagaGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctct           >  2:790495/1‑150 (MQ=255)
                                                                                                                                           ctttttCTCGCTCTCCTTTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctc          <  2:588085/146‑1 (MQ=255)
                                                                                                                                           agagTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctc          <  1:426511/150‑1 (MQ=255)
                                                                                                                                           agagTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGACTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctc          <  1:288285/150‑1 (MQ=255)
                                                                                                                                            gagTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctct         >  2:692869/1‑150 (MQ=255)
                                                                                                                                               ttCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATTGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTa      <  2:206745/150‑1 (MQ=255)
                                                                                                                                                tctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAg     <  1:821060/150‑1 (MQ=255)
                                                                                                                                                 ctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGa    >  1:1363989/1‑150 (MQ=255)
                                                                                                                                                   cGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATg  <  1:1263998/150‑1 (MQ=255)
                                                                                                                                                   cGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATg  >  2:1405186/1‑150 (MQ=255)
                                                                                                                                                   cGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATg  <  1:725554/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATG  >  NC_000913/1216972‑1217268

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: