Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,217,116 T→G 6.3% 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 6.3% 356.6 / 13.7 125intergenic (+120/‑212)ymgC/pdeGprotein YmgC/putative c‑di‑GMP phosphodiesterase PdeG
Reads supporting (aligned to +/- strand):  ref base T (52/65);  new base G (2/6);  total (54/71)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 131 total aligned reads displayed.

AATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAG  >  NC_000913/1216970‑1217265
                                                                                                                                                  |                                                                                                                                                     
aaTGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtctc                                                                                                                                                    <  2:1060989/150‑1 (MQ=255)
aaTGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtctc                                                                                                                                                    <  1:815741/150‑1 (MQ=255)
 aTGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCg                                                                                                                                                   >  1:1051689/1‑150 (MQ=255)
  agACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTCTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGc                                                                                                                                                  <  2:652946/149‑1 (MQ=255)
   gACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGct                                                                                                                                                 <  1:659694/150‑1 (MQ=255)
     cAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGctct                                                                                                                                               <  2:796896/150‑1 (MQ=255)
      aGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCtctc                                                                                                                                              <  1:835673/150‑1 (MQ=255)
       gggAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTcc                                                                                                                                             <  1:404055/150‑1 (MQ=255)
       gggAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTcc                                                                                                                                             <  1:1152933/150‑1 (MQ=255)
           aTCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTc                                                                                                                                         >  2:59682/1‑150 (MQ=255)
           aTCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTc                                                                                                                                         >  2:53597/1‑150 (MQ=255)
            tCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                                                                        <  2:1504235/150‑1 (MQ=255)
            tCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                                                                        <  1:1372537/150‑1 (MQ=255)
            tCCGTGCTCTCTTAGAACTTGAGTAAGGTTGTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                                                                        <  2:985748/150‑1 (MQ=255)
            tCCGTGCGCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                                                                        <  2:1537512/150‑1 (MQ=255)
             ccGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCtt                                                                                                                                       <  1:917238/150‑1 (MQ=255)
              cGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTg                                                                                                                                      <  1:1037334/150‑1 (MQ=38)
               gTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGa                                                                                                                                     <  1:212308/150‑1 (MQ=38)
                 gCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACg                                                                                                                                   >  1:1662712/1‑150 (MQ=255)
                   tctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGcc                                                                                                                                 <  1:1411474/150‑1 (MQ=255)
                     tctTAGACCTTGACTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAACGTGAGCTTTTCCGGCGTTATGCTTTCTTATTATGGCCGCAATATCAGGTGTCAAGAATGGAGAGTGCTCGCGCGCAATTCTTGACGCCTg                                                                                                                               <  2:203036/150‑1 (MQ=255)
                     tctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTg                                                                                                                               <  2:1229121/150‑1 (MQ=255)
                      ctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGa                                                                                                                              >  1:1644972/1‑150 (MQ=255)
                       ttAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGACACCTGat                                                                                                                             <  1:256427/150‑1 (MQ=37)
                               tGAGTAAGGTTTTAGAATAAACATCAGCTTTATCACCATGCTTATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAATAATGTAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGcc                                                                                                                     <  2:890995/150‑1 (MQ=255)
                                 aGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTa                                                                                                                   >  2:1241043/1‑150 (MQ=255)
                                   tAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAAc                                                                                                                 <  2:220837/150‑1 (MQ=255)
                                   tAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGAAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAAc                                                                                                                 <  2:642941/150‑1 (MQ=255)
                                     aGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACtt                                                                                                               <  1:755823/150‑1 (MQ=255)
                                       gTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCCCCATTCTTGACGCCTGATATCCCGCCTAACTTAt                                                                                                             <  1:820173/150‑1 (MQ=255)
                                         tttAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTAttt                                                                                                           <  1:1241043/150‑1 (MQ=255)
                                           tAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTg                                                                                                         <  1:1162076/150‑1 (MQ=255)
                                            aGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGt                                                                                                        >  2:1659394/1‑150 (MQ=255)
                                              aaTAAACATCAGCTGTATCACCATGCTGATGAAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGAAATATCAGGTGTCAATAATGGATAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTAc                                                                                                      <  2:1626738/150‑1 (MQ=255)
                                                  aaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGcc                                                                                                  <  2:1427830/150‑1 (MQ=255)
                                                    cATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTa                                                                                                <  2:669421/150‑1 (MQ=255)
                                                      tCAGCTGTATCACCATGCTGATGCAAATTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACa                                                                                              <  2:561799/150‑1 (MQ=255)
                                                      tCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACa                                                                                              <  1:1644329/150‑1 (MQ=255)
                                                         gCTGTATCAACATGCTGATGAAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTACATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTAATGCCTACATTa                                                                                           >  2:86369/1‑150 (MQ=255)
                                                              aTCACCATGCTGATGCAAAGTGCGCATTTCAGGCGTTATGCTTTCTTATTATGCCCGCAATATCCGGTGTCAACAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTc                                                                                      <  1:172863/150‑1 (MQ=255)
                                                                cACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTATCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaa                                                                                    >  2:205230/1‑150 (MQ=255)
                                                                cACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTGCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaa                                                                                    >  2:1052252/1‑150 (MQ=255)
                                                                 aCCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaat                                                                                   <  1:44311/150‑1 (MQ=255)
                                                                 aCCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaat                                                                                   >  1:1103156/1‑150 (MQ=255)
                                                                  ccATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCaata                                                                                  >  1:1472257/1‑150 (MQ=255)
                                                                   cATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCAAACTTATTTTGTACTGCCTACATTATGCTCaataa                                                                                 >  2:1489380/1‑150 (MQ=255)
                                                                     tGCTGATGCAAAGTTAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                               <  1:1136240/150‑1 (MQ=255)
                                                                     tGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                               <  2:1051689/150‑1 (MQ=255)
                                                                     tGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGCAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                               <  2:1274216/150‑1 (MQ=255)
                                                                       cTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaattaa                                                                             <  1:689817/150‑1 (MQ=255)
                                                                          aTGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATgg                                                                          <  1:53597/150‑1 (MQ=255)
                                                                          aTGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATgg                                                                          <  2:1297455/150‑1 (MQ=255)
                                                                          aTGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATgg                                                                          <  1:59682/150‑1 (MQ=255)
                                                                           tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                         >  1:326851/1‑150 (MQ=255)
                                                                           tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                         <  1:1437404/150‑1 (MQ=255)
                                                                           tGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATggg                                                                         >  1:662311/1‑150 (MQ=255)
                                                                              aaaGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                               <  2:908359/141‑1 (MQ=255)
                                                                              aaaGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaatt                                                                               >  1:908359/1‑141 (MQ=255)
                                                                               aaGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAg                                                                     >  2:1037393/1‑150 (MQ=255)
                                                                                aGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGa                                                                    >  2:1526783/1‑150 (MQ=255)
                                                                                      cATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGa                                                              >  1:960034/1‑150 (MQ=255)
                                                                                      cATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGa                                                              >  1:540973/1‑150 (MQ=255)
                                                                                       aTTGCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGCAGAGTTCTCGCTCTCCTTTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGAt                                                             <  2:1411015/150‑1 (MQ=255)
                                                                                        tttCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATg                                                            >  2:1189023/1‑150 (MQ=255)
                                                                                          tCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGga                                                          >  1:967080/1‑150 (MQ=255)
                                                                                                 ttATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAAtt                                                   >  1:1616557/1‑150 (MQ=255)
                                                                                                  tATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAATATTGTAGAGTTCTCGCTATCCATTATTGACGCCTTATATCCCGCCTAACTTATTTTGTACTTCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAAttt                                                  <  1:86369/150‑1 (MQ=255)
                                                                                                     gCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACAGTATGCTCAATAATTAATGGGTTAGATGAAGATGGAGAATATTTTAt                                               >  2:1428392/1‑150 (MQ=255)
                                                                                                      cTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATc                                              >  1:1301343/1‑150 (MQ=255)
                                                                                                          cTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCg                                          >  2:253056/1‑150 (MQ=255)
                                                                                                          cTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCg                                          <  1:1141797/150‑1 (MQ=255)
                                                                                                               tatGGCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGGCTGATATCCCGCCTAACTTATGTTGTACTGCCCACATTATGCTCAATACTTAATGGGTTAGATCAAGATGGAGAAAATTCTATCTTCGGCGTg                                     <  2:1238307/150‑1 (MQ=255)
                                                                                                                 tGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATgg                                                           <  1:455550/126‑1 (MQ=255)
                                                                                                                   tCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGata                                 >  1:792288/1‑150 (MQ=255)
                                                                                                                    ccGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGatat                                >  2:665497/1‑150 (MQ=255)
                                                                                                                    ccGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGatat                                >  2:872355/1‑150 (MQ=255)
                                                                                                                    ccGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGatat                                >  1:1447654/1‑150 (MQ=255)
                                                                                                                      gCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGGCGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATAttt                              >  2:387424/1‑150 (MQ=255)
                                                                                                                        gatatCAGGCGTCAAGAATGGAGAGCGAGAACTCTCCATTCTTGACACCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCt                            >  2:859731/2‑150 (MQ=37)
                                                                                                                        aaTATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCt                            <  2:624915/150‑1 (MQ=255)
                                                                                                                         atatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCtt                           <  2:1103156/150‑1 (MQ=255)
                                                                                                                          tatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTa                          >  2:446948/1‑150 (MQ=255)
                                                                                                                             cAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGc                       >  1:749124/1‑150 (MQ=255)
                                                                                                                                   tCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTa                 <  2:46551/150‑1 (MQ=255)
                                                                                                                                   tCAAGAAATGAGAGTTCTAGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATATTATCTTCGGAGTGGATATTTCTTATTCACTCTa                 >  2:1663932/1‑150 (MQ=255)
                                                                                                                                    cAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATAAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTAc                >  2:128600/1‑150 (MQ=255)
                                                                                                                                     aaGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTAcc               <  2:1644972/150‑1 (MQ=255)
                                                                                                                                      aGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTActa              >  2:1142305/1‑148 (MQ=255)
                                                                                                                                      aGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCa              <  2:529267/150‑1 (MQ=255)
                                                                                                                                        aaTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATg            >  2:465747/1‑150 (MQ=255)
                                                                                                                                         aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc           <  1:168489/150‑1 (MQ=255)
                                                                                                                                         aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc           <  1:1189023/150‑1 (MQ=255)
                                                                                                                                         aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc           >  1:381463/1‑150 (MQ=255)
                                                                                                                                         aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc           >  1:1228228/1‑150 (MQ=255)
                                                                                                                                         aTGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGc           <  2:1472257/150‑1 (MQ=255)
                                                                                                                                            gagaGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGctct        <  2:540973/150‑1 (MQ=255)
                                                                                                                                                gTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGt    >  1:961628/1‑150 (MQ=255)
                                                                                                                                                gTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGt    <  1:842309/150‑1 (MQ=255)
                                                                                                                                                  tctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAg  <  2:1662712/150‑1 (MQ=255)
                                                                                                                                                  tctcGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAg  <  1:100120/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAG  >  NC_000913/1216970‑1217265

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