Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,217,126 C→A 9.0% intergenic (+130/‑202) 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,1260CA 9.0% 296.7 / 3.6 133intergenic (+130/‑202)ymgC/ycgGBlue light, low temperature and stress induced protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (73/45);  new base A (2/13);  total (75/58)
Fisher's exact test for biased strand distribution p-value = 4.87e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCG  >  NC_000913/1216978‑1217274
                                                                                                                                                    |                                                                                                                                                    
ggAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCa                                                                                                                                                     <  2:1476846/150‑1 (MQ=255)
ggAATACGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATTCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTATTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCa                                                                                                                                                     >  1:351183/1‑150 (MQ=255)
 gAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCAt                                                                                                                                                    >  1:942852/1‑150 (MQ=255)
 gAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCAt                                                                                                                                                    <  1:1479885/150‑1 (MQ=255)
    tCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                                                                                 <  1:475307/150‑1 (MQ=255)
    tCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                                                                                 <  1:155468/150‑1 (MQ=255)
     ccGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGTTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCtt                                                                                                                                                >  1:259142/1‑150 (MQ=255)
      cgggCTCACTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTg                                                                                                                                               <  1:1063841/147‑1 (MQ=255)
         gCTCTCTTAGCCCTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCCGGTGCCAAGCATGGAGAGTTCTCGCTCTCCATTCTTGACg                                                                                                                                            <  1:2073342/150‑1 (MQ=255)
          ctctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGc                                                                                                                                           >  1:983316/1‑150 (MQ=255)
             tctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTg                                                                                                                                        >  1:923582/1‑150 (MQ=255)
              ctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGa                                                                                                                                       >  1:543601/1‑150 (MQ=255)
              ctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTACATTCTTGACGACTGa                                                                                                                                       >  2:1303245/1‑150 (MQ=255)
               ttAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGat                                                                                                                                      >  1:982437/1‑150 (MQ=255)
                tAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGata                                                                                                                                     <  2:1406455/150‑1 (MQ=255)
                tAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGata                                                                                                                                     >  1:2028882/1‑150 (MQ=255)
                 aGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGatat                                                                                                                                    >  2:1714528/1‑150 (MQ=255)
                 aGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGataa                                                                                                                                    >  2:737232/1‑149 (MQ=255)
                       ttattAAGTTTTTACAATAAAAATCAGCTTAATCCCCATGCTTATGAAAAGTTAGCATTTCAGGCGTTTTGCTTTCTTTTTATGTCCCAAATATCAGGGGTCAATAATGGACAGTTCTCGCTCTCCTTTCTTGACGCCTGATATCCCGcc                                                                                                                              <  2:960661/146‑1 (MQ=255)
                          gTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTaa                                                                                                                           >  1:106938/1‑150 (MQ=255)
                           tAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAAc                                                                                                                          <  2:337027/150‑1 (MQ=255)
                            atggTTTTATAATAAACATCAGCTGTATCACCATGCTGATGAAAAGTGTGCATTTCTGGCTTTATGGTTTCTTATTATGTCCGAAATATCTGGTGTCAAGATTGGAGAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACt                                                                                                                         <  2:1930978/148‑1 (MQ=255)
                              ggTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATCTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTa                                                                                                                       >  1:1627457/1‑150 (MQ=255)
                                 tttAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGAAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTAttt                                                                                                                    >  2:75918/1‑150 (MQ=255)
                                  ttAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTAtttt                                                                                                                   <  1:1093472/150‑1 (MQ=255)
                                   tAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCACCTAACTTATTTTg                                                                                                                  <  2:843342/150‑1 (MQ=255)
                                   tAGAAAAAACATCAGCTTTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGGCAAGAATGGAGCGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTg                                                                                                                  <  1:778300/150‑1 (MQ=255)
                                     taaTAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTa                                                                                                                <  1:1134107/149‑1 (MQ=255)
                                     gAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTa                                                                                                                >  2:619995/1‑150 (MQ=255)
                                     gAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTa                                                                                                                >  1:933191/1‑150 (MQ=255)
                                     gAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAACATTTCAGGCGTTATGCTTTCTTATTATGTCCTCAATATCAGGTGTCAAGAATTGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTa                                                                                                                >  1:1694985/1‑150 (MQ=255)
                                      aaTAAAAATATGCTTTATCACCTTGGTTTTGAAATGTTAGAATTTCTGGCTTTATGTTTTATTTTTATGTCCTAAATATCTGGTGTCAATAATGGTTAGTTCTCGCTCTCATTTTTTGACGCCTGATATCCCGCCTAACTTATTTTGTAc                                                                                                               <  2:942852/150‑1 (MQ=255)
                                        tAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCACGCCTAACTTATTTTGTACTg                                                                                                             >  2:1508570/1‑150 (MQ=255)
                                        tAAAAATCAGCTTTATCACCATGCTGTTGCAAAGTTAGCTTTTCAGGCGTTATCGTTTCTTTTTATGTCCGCAATATCAGGTGTCAAGAATGGAGCGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTg                                                                                                             <  2:633190/150‑1 (MQ=255)
                                         aaaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGc                                                                                                            <  1:1341526/150‑1 (MQ=255)
                                                aGCTGTATCACCATGCTGATGCATAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCTCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTTTACTGCCTACAtt                                                                                                     >  2:1727363/1‑150 (MQ=255)
                                                   tGTATCACCATGCTGATGCAAAGTGAGCATTTCCGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTGCACGCCAGATATACCGGCTAAAATATTTTGTACTTCCTACATTATg                                                                                                  >  1:1862203/1‑150 (MQ=255)
                                                   tGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAATAATGGAGAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATg                                                                                                  <  2:511487/150‑1 (MQ=255)
                                                    gTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGc                                                                                                 >  1:633203/1‑150 (MQ=255)
                                                      aTCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCta                                                                                               >  2:2031441/1‑149 (MQ=255)
                                                         aCCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCATATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTAGCTCTCCATTCTTGACGACTGATATCCCGCCTAACTTATTTTGTACTGCATACATTATGCTCaat                                                                                            >  1:228925/1‑150 (MQ=255)
                                                           cATGCTTTTGCAATGTGAGATTTTCAGGGGTTATGTTTTATTATTATGTCCGAAATATCAGGTGTCATCAATGGAGAGTTCTGGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaataa                                                                                          <  2:1778836/150‑1 (MQ=255)
                                                           cATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGGGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaataa                                                                                          <  2:621710/150‑1 (MQ=255)
                                                           cATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaataa                                                                                          <  2:1376951/150‑1 (MQ=255)
                                                           cATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGAAATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaataa                                                                                          >  2:196571/1‑150 (MQ=255)
                                                            aTGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaataat                                                                                         <  2:1131922/150‑1 (MQ=255)
                                                             tGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATACCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                                        >  1:1207784/1‑150 (MQ=255)
                                                              gCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGTAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatta                                                                                       <  1:1438565/150‑1 (MQ=255)
                                                                   tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                                  >  1:1925121/1‑150 (MQ=255)
                                                                   tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGATCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                                  >  1:36604/1‑150 (MQ=255)
                                                                    gCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGt                                                                                 >  2:180256/1‑150 (MQ=255)
                                                                          tGACCATTTCAGGCGTTATGCTTTCTTTTTATGTCCGCAATATCCGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATc                                                                           <  2:165765/150‑1 (MQ=255)
                                                                           gAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTGACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCa                                                                          >  2:1104149/1‑150 (MQ=255)
                                                                           gAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCa                                                                          >  2:1875497/1‑150 (MQ=255)
                                                                           gAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTATAGCTATCCATTCTTGACGCCTGATTTACCGCCTAAGTTTTTTTGTAATGCCTACATTTTGATCATTAATTAATGGGTTAGATCa                                                                          >  1:364510/1‑150 (MQ=255)
                                                                               aTTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGAt                                                                      <  1:114496/150‑1 (MQ=255)
                                                                                tttCTGGCGTTTTGCTTTTTTATTATGTCCGAATTATCCGTTGTCAACAATGTAGAGTTCTCGCTCTCAATTATTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATg                                                                     <  1:1773241/150‑1 (MQ=255)
                                                                                tttCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATg                                                                     <  2:1490076/150‑1 (MQ=255)
                                                                                tttCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATg                                                                     >  1:1547777/1‑150 (MQ=255)
                                                                                 ttCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGATTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATgg                                                                    <  2:403292/150‑1 (MQ=255)
                                                                                  tCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGga                                                                   <  2:2028882/150‑1 (MQ=255)
                                                                                            tgatttCTTTTTATGTCCGAAATATCTGGTTTAAAGAATGTAGTGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTa                                                         <  1:903058/147‑1 (MQ=255)
                                                                                              cTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGATAGTTCTCGCTCTCCATTCTTGCCGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATc                                                       >  2:1596849/1‑150 (MQ=255)
                                                                                              cTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATc                                                       >  2:724103/1‑150 (MQ=255)
                                                                                              cTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTACATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGATCAATAATTAATGGGTTAGATCAAGATTGAGAAAATTTTaaa                                                       >  2:1587834/1‑148 (MQ=255)
                                                                                                 tCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTc                                                    >  1:1772587/1‑150 (MQ=255)
                                                                                                       tatGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTTTTTTGTACTGCCTACATGTTGCTCTATTATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCggaggg                                              >  1:1142921/1‑146 (MQ=255)
                                                                                                        atGTCCGCAATATCAGGTGGCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTgg                                             <  2:1080380/150‑1 (MQ=255)
                                                                                                          gTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTTCTGCCTACATTATGCTCAATAATTAATGGGTTAGTTCAAGATGGAGAAAATTTTATCTTCGGCGTGGat                                           >  1:1756064/1‑150 (MQ=255)
                                                                                                          gTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGat                                           >  2:1198819/1‑150 (MQ=255)
                                                                                                             cGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGCCGCCTGATATCGGGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGGTCAAGATGGAGAAAATTTTATCTTCGGCGTGGATAtt                                        >  1:1940862/1‑150 (MQ=255)
                                                                                                                 atatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCtt                                    >  1:771733/1‑150 (MQ=255)
                                                                                                                  tatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTa                                   <  1:1349303/150‑1 (MQ=255)
                                                                                                                   atCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTTTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTAt                                  >  2:2007638/1‑150 (MQ=255)
                                                                                                                            cAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTAc                         >  2:472242/1‑150 (MQ=255)
                                                                                                                            cAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTAc                         >  2:1986630/1‑150 (MQ=255)
                                                                                                                            cAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTAc                         >  2:1730877/1‑150 (MQ=255)
                                                                                                                            cAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGACAAAATTTTATCTTCGGCGTGGATATTTCTTATGCGCTCTAc                         >  2:382354/1‑150 (MQ=255)
                                                                                                                             aaGAATGTATAGTTCTCTCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTAcc                        <  2:914342/150‑1 (MQ=255)
                                                                                                                               gAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCAt                      >  1:1783929/1‑150 (MQ=255)
                                                                                                                               gAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTAGGCGTGGATATTTCTTATGCTCTCTACCAt                      >  2:975876/1‑150 (MQ=255)
                                                                                                                                aaTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTAGATTATGCTCAATAATTAATGGGTTAGATCTAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATg                     <  1:1385339/150‑1 (MQ=255)
                                                                                                                                 ttGGCGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc                    <  2:1682102/149‑1 (MQ=255)
                                                                                                                                 aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc                    >  1:1589206/1‑150 (MQ=255)
                                                                                                                                    gagaGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctct                 >  1:2001495/1‑150 (MQ=255)
                                                                                                                                     ttatttgtcGCTATCATTTATTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGACTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctc                <  1:737232/143‑1 (MQ=255)
                                                                                                                                      gagTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATATTCGGCGTGGATATTTCTTATGCTCTCTACCATGctctct               >  2:1386098/1‑150 (MQ=255)
                                                                                                                                       tgTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTg              <  2:913974/149‑1 (MQ=255)
                                                                                                                                       agTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTg              <  1:1279725/150‑1 (MQ=255)
                                                                                                                                        gTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGt             <  1:740068/150‑1 (MQ=255)
                                                                                                                                          tctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGTTATTTCTTATGCTCTCTACCATGCTCTCTGTAg           >  2:1141637/1‑150 (MQ=255)
                                                                                                                                          tctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTat           >  2:1937070/1‑149 (MQ=255)
                                                                                                                                          tctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAg           >  1:1201268/1‑150 (MQ=255)
                                                                                                                                            tcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGAt         >  1:1922952/1‑150 (MQ=255)
                                                                                                                                                tctcCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTATCTGTAGATGCCa     >  2:1078569/1‑150 (MQ=255)
                                                                                                                                                 ctcCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAg    <  1:777305/150‑1 (MQ=255)
                                                                                                                                                  tcCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAgc   <  2:1685511/150‑1 (MQ=255)
                                                                                                                                                   ccATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCTgcg  >  2:816415/1‑150 (MQ=255)
                                                                                                                                                   ccATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGACAgcg  >  1:826471/1‑150 (MQ=255)
                                                                                                                                                   ccATTCTTGAAGCCTGATATCCCGCCTAACTTATATTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGAGTGTATATTTCTTATGATCTCTACCATGCTCTCTGTAGATGCCAgcg  >  2:1040265/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCG  >  NC_000913/1216978‑1217274

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: