Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,217,116 T→G 5.2% intergenic (+120/‑212) ymgC → / → pdeG protein YmgC/putative c‑di‑GMP phosphodiesterase PdeG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,217,1160TG 5.2% 433.2 / 11.3 149intergenic (+120/‑212)ymgC/pdeGprotein YmgC/putative c‑di‑GMP phosphodiesterase PdeG
Reads supporting (aligned to +/- strand):  ref base T (81/60);  new base G (3/5);  total (84/65)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

ATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCT  >  NC_000913/1216967‑1217259
                                                                                                                                                     |                                                                                                                                               
   aaTGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtctc                                                                                                                                              >  1:244573/1‑150 (MQ=255)
     tGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGc                                                                                                                                            >  2:1544007/1‑150 (MQ=255)
     tGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGc                                                                                                                                            >  1:1539986/1‑150 (MQ=255)
              aTCCTTTCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCTCTCTCCATTc                                                                                                                                   <  2:1910278/150‑1 (MQ=255)
               tCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                                                                  <  1:1094504/150‑1 (MQ=255)
                 cGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTg                                                                                                                                >  1:193075/1‑150 (MQ=255)
                 cGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTg                                                                                                                                >  1:1985383/1‑150 (MQ=255)
                 cGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTg                                                                                                                                <  2:742997/150‑1 (MQ=255)
                 cGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTg                                                                                                                                <  1:4279/150‑1 (MQ=39)
                   tGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGAc                                                                                                                              >  1:766833/1‑150 (MQ=255)
                     ctctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGAcac                                                                                                                            >  1:1493759/1‑148 (MQ=38)
                     ctctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGAcac                                                                                                                            <  1:1940727/150‑3 (MQ=38)
                      tctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGcc                                                                                                                           <  2:1278718/150‑1 (MQ=255)
                      tctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGcc                                                                                                                           >  2:1003357/1‑150 (MQ=255)
                      tctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGACAcc                                                                                                                           <  2:1493759/150‑1 (MQ=37)
                      tctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGACAcc                                                                                                                           <  1:1862509/150‑1 (MQ=37)
                       ctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt                                                                                                                          <  2:442542/150‑1 (MQ=255)
                       ctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt                                                                                                                          >  2:2016273/1‑150 (MQ=255)
                       ctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt                                                                                                                          >  2:2312036/1‑150 (MQ=255)
                          ttAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGat                                                                                                                       >  2:1438144/1‑150 (MQ=255)
                             gAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATc                                                                                                                    >  2:2439691/1‑150 (MQ=255)
                                cTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCg                                                                                                                 >  1:1557574/1‑150 (MQ=255)
                                  tGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAATAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGcc                                                                                                               <  2:558318/150‑1 (MQ=255)
                                  tGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGcc                                                                                                               >  2:1175331/1‑150 (MQ=255)
                                    aGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTa                                                                                                             <  1:2133010/150‑1 (MQ=255)
                                        aGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACtt                                                                                                         >  1:955258/1‑150 (MQ=255)
                                           ttttAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTAtt                                                                                                      <  1:1179253/150‑1 (MQ=255)
                                              tAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTAttttt                                                                                                   >  1:2344766/1‑149 (MQ=255)
                                              tAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCATAACTTATTTTg                                                                                                   >  1:1596578/1‑150 (MQ=255)
                                               aGAATAAACATCAGCTTTATCACCATGCTTATGCAAATTTATCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAATAATGGAGAGTTCTCGCTCTCAATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGt                                                                                                  <  1:1570493/150‑1 (MQ=255)
                                               aGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGt                                                                                                  <  1:664305/150‑1 (MQ=255)
                                               aGAATAAACATAAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGt                                                                                                  >  2:2057833/1‑150 (MQ=255)
                                                gAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTa                                                                                                 <  2:2476911/150‑1 (MQ=255)
                                                    aaaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACt                                                                                                          <  2:1626046/137‑1 (MQ=255)
                                                    aaaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGc                                                                                             <  1:396285/150‑1 (MQ=255)
                                                    aaaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGc                                                                                             >  1:1727803/1‑150 (MQ=255)
                                                    aaaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGc                                                                                             <  2:10705/150‑1 (MQ=255)
                                                     aaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGcc                                                                                            >  2:1383610/1‑150 (MQ=255)
                                                      aCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGATATTTCTAGCTCTCCATTCTTGAAGCCTGATATCCCGCCTAACTTATTTTGTACTGCCt                                                                                           <  2:868128/150‑1 (MQ=255)
                                                      aCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCt                                                                                           <  1:1285993/150‑1 (MQ=255)
                                                        aTCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTAc                                                                                         <  1:1438144/150‑1 (MQ=255)
                                                        aTCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTAc                                                                                         >  1:295891/1‑150 (MQ=255)
                                                         tCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTGTGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACa                                                                                        >  1:1366639/1‑150 (MQ=255)
                                                         tCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACa                                                                                        >  2:483484/1‑150 (MQ=255)
                                                          cAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACAt                                                                                       >  2:1207379/1‑150 (MQ=255)
                                                          cAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACAt                                                                                       >  2:1265275/1‑150 (MQ=255)
                                                             cggTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTAt                                                                                    <  1:87103/148‑1 (MQ=255)
                                                               gTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTCGTACTGCCTACATTATGc                                                                                  <  1:1805984/150‑1 (MQ=255)
                                                                 aTCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTc                                                                                <  2:1624087/150‑1 (MQ=255)
                                                                 aTCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTc                                                                                <  1:160463/150‑1 (MQ=255)
                                                                   cACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaa                                                                              <  2:799214/150‑1 (MQ=255)
                                                                     ccATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaata                                                                            >  2:1351666/1‑150 (MQ=255)
                                                                     ccATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaata                                                                            >  1:1894669/1‑150 (MQ=255)
                                                                       aTGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaataat                                                                          <  1:1012256/150‑1 (MQ=255)
                                                                        tGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                         <  1:1544007/150‑1 (MQ=255)
                                                                        tGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                         >  2:1005336/1‑150 (MQ=255)
                                                                           tGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaattaat                                                                      >  2:2342289/1‑150 (MQ=255)
                                                                           tGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaattaat                                                                      <  1:388142/150‑1 (MQ=255)
                                                                            gATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATg                                                                     >  2:2379716/1‑150 (MQ=255)
                                                                              tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                   <  2:937861/150‑1 (MQ=255)
                                                                              tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                   >  1:2195093/1‑150 (MQ=255)
                                                                              tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                   >  1:2398052/1‑150 (MQ=255)
                                                                               gCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGt                                                                  >  1:2255720/1‑150 (MQ=255)
                                                                                cAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTATTTTGTACTGCCGACATTATGCTGAATTATTAATGGGttt                                                                >  2:1782463/1‑149 (MQ=255)
                                                                                       aGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCTTACATTATGCTCAATAATTAATGGGTTAGATCaa                                                          <  2:147807/150‑1 (MQ=255)
                                                                                        gCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCGTACTTATTTTGAACTGCATACATTATGGTCAATAATTAATGGGTTAGATCAAg                                                         >  2:230528/1‑150 (MQ=255)
                                                                                        gCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTATCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAg                                                         >  2:1889227/1‑150 (MQ=255)
                                                                                         cATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGa                                                        <  1:1270600/150‑1 (MQ=255)
                                                                                          aTTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGAt                                                       >  2:1606683/1‑150 (MQ=255)
                                                                                             tCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGga                                                    <  1:1351666/150‑1 (MQ=255)
                                                                                               aGGCGTTATGCTTTCTTATTATGTCCGCATTATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGgaga                                                  <  2:7060/150‑1 (MQ=255)
                                                                                                ggCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGaa                                                 >  2:1556151/1‑150 (MQ=255)
                                                                                                    ttATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAAtt                                             <  1:2108858/150‑1 (MQ=255)
                                                                                                       tGCTTTCTTATTATGTCCGCAATATCAGTTGTCAAGAATGGAGATTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGTTTAGATCAAGATGGAGAAAATTTTa                                          <  2:696270/150‑1 (MQ=255)
                                                                                                         cTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATc                                        >  1:1773782/1‑150 (MQ=255)
                                                                                                              ttattaTGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCgg                                   <  2:222501/150‑1 (MQ=255)
                                                                                                                  tatGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTg                               <  2:1985383/150‑1 (MQ=255)
                                                                                                                   atGTCCGCAATATCCGGGGCCAAGAATGGAGAGTGCTCGCGCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTgg                              <  1:2422639/150‑1 (MQ=255)
                                                                                                                     gTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGCTATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATGTTCGGCGTGGat                            >  2:1732509/1‑150 (MQ=255)
                                                                                                                     gTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGat                            >  1:924735/1‑150 (MQ=255)
                                                                                                                     gTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGat                            >  2:2283421/1‑150 (MQ=255)
                                                                                                                      tCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGata                           >  2:1722222/1‑150 (MQ=255)
                                                                                                                      cccGCAATATCAGGTGTCAAGAATGTAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGata                           <  1:1968963/149‑1 (MQ=255)
                                                                                                                       ccGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGatat                          >  1:751808/1‑150 (MQ=255)
                                                                                                                        cGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATAtt                         <  2:1366639/150‑1 (MQ=255)
                                                                                                                         gCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATgg                                                                    <  1:817988/106‑1 (MQ=255)
                                                                                                                         gCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATAAAGATGGAGAAAATTTTATGTTCGGCGTGGATtttt                        >  2:2145958/1‑150 (MQ=255)
                                                                                                                             tatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTa                    <  1:2016273/150‑1 (MQ=255)
                                                                                                                              atCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTAt                   <  1:1197461/150‑1 (MQ=255)
                                                                                                                                cAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGc                 >  2:170946/1‑150 (MQ=255)
                                                                                                                                atgttttCAATAATGTAGTTTTCTCGCTCTCAATTTTTGACGACTTATATCCCGCCTAACTTATTTTGTACTTCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGc                 <  1:635681/144‑1 (MQ=255)
                                                                                                                                     gtCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGctctct            <  1:170946/150‑1 (MQ=255)
                                                                                                                                      tCAAGAATGGAGAGCGAGAACTCTCCATTCTTGACACCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTa           >  1:1091217/1‑150 (MQ=38)
                                                                                                                                       cAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTAc          >  2:2306038/1‑150 (MQ=255)
                                                                                                                                         aGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCa        >  2:1756571/1‑150 (MQ=255)
                                                                                                                                          gAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCAt       <  2:955258/150‑1 (MQ=255)
                                                                                                                                            aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc     >  1:2474299/1‑150 (MQ=255)
                                                                                                                                            aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc     <  1:1005336/150‑1 (MQ=255)
                                                                                                                                               gagaGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctct  >  2:2278175/1‑150 (MQ=255)
                                                                                                                                               gagaGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctct  >  2:2270894/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCT  >  NC_000913/1216967‑1217259

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