Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,217,117 C→A 5.2% intergenic (+121/‑211) 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,1170CA 5.2% 274.1 / 4.6 106intergenic (+121/‑211)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/39);  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 = 6.17e-01

Only 100 of 110 total aligned reads displayed.

TAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGA  >  NC_000913/1216969‑1217266
                                                                                                                                                    |                                                                                                                                                     
taATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtct                                                                                                                                                      >  1:1393003/1‑150 (MQ=255)
   tGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGc                                                                                                                                                   >  1:1557830/1‑150 (MQ=255)
   tGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGc                                                                                                                                                   >  2:1722044/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:1202391/1‑150 (MQ=255)
                     ctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt                                                                                                                                 >  2:951223/1‑150 (MQ=255)
                      tctTAGAACTTGAGTAAGGTATTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGACACCTg                                                                                                                                <  1:167858/150‑1 (MQ=37)
                       ctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGa                                                                                                                               >  2:632002/1‑150 (MQ=255)
                       ctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGa                                                                                                                               <  2:1476078/150‑1 (MQ=255)
                         tAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGata                                                                                                                             >  2:334345/1‑150 (MQ=255)
                           gAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATc                                                                                                                           >  1:1709551/1‑150 (MQ=255)
                           gAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATc                                                                                                                           >  1:1402313/1‑150 (MQ=255)
                             aCTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATccc                                                                                                                         >  2:6082/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)
                                           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)
                                                       tCAGTTGTTTCACCTTTCTTTTTCAAATTTAGAATTTCTGGGGTTATGCTTTCTTATTATGTCCGAAATATCAGGTGTAAATAATGTAGATTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACa                                                                                               <  2:1675665/150‑1 (MQ=255)
                                                          tcTGTATCACCATGCTTATGCAAATTGAGCATTTCAGGAGTTATGATTTATTATTTTGTCGGAAATATAAGGTGTCAAGAATGGACAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTa                                                                                            <  1:1094707/149‑1 (MQ=255)
                                                             gTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGc                                                                                         >  2:434917/1‑150 (MQ=255)
                                                              tATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTTTGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGTAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCt                                                                                        <  1:1007475/150‑1 (MQ=255)
                                                                tCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCACAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCa                                                                                      >  1:651556/1‑150 (MQ=255)
                                                                 cACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaa                                                                                     >  2:618121/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)
                                                                         tGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCAATTCTTGACGACTTATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatttat                                                                             >  1:626148/1‑150 (MQ=255)
                                                                          gATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATg                                                                            >  1:95402/1‑150 (MQ=255)
                                                                            tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                          >  2:1285600/1‑150 (MQ=255)
                                                                            tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                          >  1:1813952/1‑150 (MQ=255)
                                                                                    gAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCa                                                                  >  2:1350351/1‑150 (MQ=255)
                                                                                    gAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCa                                                                  <  1:194321/150‑1 (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)
                                                                                                 gTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAAt                                                     >  1:1236578/1‑150 (MQ=255)
                                                                                                  ttATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAAtt                                                    >  1:1252230/1‑150 (MQ=255)
                                                                                                   tATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAAttt                                                   <  2:6481/150‑1 (MQ=255)
                                                                                                     tgatttCTTATTATGTCGCCAATATCAGGTGTCAAGAATGGAGAGTTCTCTCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTa                                                 <  1:180022/147‑1 (MQ=255)
                                                                                                           cTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCg                                           <  1:1201271/150‑1 (MQ=255)
                                                                                                           attattaTGTCCGCAATATCAGGTGTCAATAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCg                                           <  2:741588/149‑1 (MQ=255)
                                                                                                               ttatGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGt                                       >  1:528647/1‑150 (MQ=255)
                                                                                                                tatGTCCGAAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTg                                      <  2:456886/150‑1 (MQ=255)
                                                                                                                     ccGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGCCGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATAAAGATGGAGAAAATTTTATCTTCGGCGTGGatat                                 >  2:1674113/1‑150 (MQ=255)
                                                                                                                     ccGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGatat                                 >  1:750723/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)
                                                                                                                                              agagTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctc        <  1:426511/150‑1 (MQ=255)
                                                                                                                                              agagTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGACTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctc        <  1:288285/150‑1 (MQ=255)
                                                                                                                                               gagTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctct       >  2:692869/1‑150 (MQ=255)
                                                                                                                                                   tctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAg   <  1:821060/150‑1 (MQ=255)
                                                                                                                                                    ctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGa  >  1:1363989/1‑150 (MQ=255)
                                                                                                                                                    ctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGa  >  1:1435019/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGA  >  NC_000913/1216969‑1217266

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: