Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 282,518 G→A 8.3% G81S (GGC→AGC)  yagE → 2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913282,5180GA 8.3% 365.1 / 14.3 122G81S (GGC→AGC) yagE2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (53/57);  new base A (5/5);  total (59/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

Only 100 of 124 total aligned reads displayed.

CGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCC  >  NC_000913/282369‑282661
                                                                                                                                                     |                                                                                                                                               
cgcTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGCTCg                                                                                                                                                 >  2:96643/1‑150 (MQ=255)
cgcTGATCGACGATCTGATCAAAGAAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCCGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTAGATCGTCGCGTGCCGGTGCTGATCg                                                                                                                                                 >  1:1616784/1‑150 (MQ=255)
  cTGATCGACGATCTGATCACAGCAGGCGTGGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               <  2:1575387/150‑1 (MQ=255)
  cTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               <  2:1195390/150‑1 (MQ=255)
  cTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               >  1:1306528/1‑150 (MQ=255)
  cTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               >  1:1210358/1‑150 (MQ=255)
    gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCAc                                                                                                                                             >  2:429901/1‑150 (MQ=255)
     aTCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCAcc                                                                                                                                            >  1:159664/1‑150 (MQ=255)
        gacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGGGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAGCACcggc                                                                                                                                         <  1:63728/150‑1 (MQ=255)
        gacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAACGCCATTGCCCGCTTTGCTATCGATCATGTCGAACGTCGCGTGCCGGTGCTGATCGGCACcggc                                                                                                                                         >  2:106415/1‑150 (MQ=255)
        gacgaTCTGATCAAAGCAGGCGGTGACCGCCTGTTCTTCCTGGGCAGCGGGGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTTGATCGTCGCGTGCCGGTGCCGATCAGCACcggc                                                                                                                                         <  1:270165/150‑1 (MQ=255)
           gaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAGCACcggcacg                                                                                                                                      <  1:1175814/150‑4 (MQ=255)
              ttgatgaATGCAGTCGTTGACGGCCTGTTATTCCTGGGGAGCGGGGGCGAGTTCTCCCAGCTCGGGGCCGAAGAGGGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCAcc                                                                                                                                   <  2:125969/149‑1 (MQ=255)
                gATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCaa                                                                                                                                 >  2:687644/1‑150 (MQ=255)
                 aTCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAAc                                                                                                                                >  1:791963/1‑150 (MQ=255)
                  tCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACg                                                                                                                               <  2:1282743/150‑1 (MQ=255)
                      aGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               <  2:124348/130‑1 (MQ=255)
                      aGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               >  1:124348/1‑130 (MQ=255)
                      aGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCg                                                                                                                           <  1:1088719/150‑1 (MQ=255)
                        cAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCggg                                                                                                                         <  2:1306528/150‑1 (MQ=255)
                            cGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAAc                                                                                                                     <  1:774224/150‑1 (MQ=255)
                              ttGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCa                                                                                                                   >  2:615380/1‑150 (MQ=255)
                              ttGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTATCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCACGGGAAACCa                                                                                                                   >  1:428735/1‑150 (MQ=255)
                               tGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCAt                                                                                                                  <  2:1638990/150‑1 (MQ=255)
                               tGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCAt                                                                                                                  <  2:1632377/150‑1 (MQ=255)
                                 aaggCCTTTTCTTCCTGGTCGGAGGTGGCGAGTTCTCCCAGCGCGGCGGCGAAGAGCGTAAAGCCATTGCCCGCTATGCTATAGATCATCTCGATCGTCGCGTGCCGGTGCTGCTCGGCACAGGCGGCACCTACGACCGGGAAACCTTCg                                                                                                                <  1:860696/148‑1 (MQ=255)
                                   ggCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGaa                                                                                                              <  1:467737/150‑1 (MQ=255)
                                    gCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTCTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCCCCGGCGGCACCAACGCCCGGGAAACCATCGAAc                                                                                                             >  2:821511/1‑150 (MQ=255)
                                       tGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGCGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGCTCGTCGCGTCCCGGTGCTGATCGGCACCGGCGGCACCACCGCCCGCGCAACCATCGAACTca                                                                                                          >  1:67847/1‑150 (MQ=255)
                                       tGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcc                                                                                                          >  1:1428525/1‑149 (MQ=255)
                                        gTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcag                                                                                                         >  2:1716167/1‑150 (MQ=255)
                                          tcttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagcc                                                                                                       <  1:1701134/150‑1 (MQ=255)
                                           cttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagcca                                                                                                      >  1:653843/1‑150 (MQ=255)
                                             tcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGCTCGGCACCGGCGGCACCGACGCCCGGGCTCCCAGCGACCTccgccagc                                                                                                    >  1:444128/1‑150 (MQ=255)
                                              ccTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCa                                                                                                   <  1:1671096/150‑1 (MQ=255)
                                               cTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAc                                                                                                  <  1:1614222/150‑1 (MQ=255)
                                                  ggCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAcgcg                                                                                               <  2:500998/150‑1 (MQ=255)
                                                  ggCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAcgcg                                                                                               <  2:314808/150‑1 (MQ=255)
                                                  ggCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGGGCCGGTGCTGAGCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAcgcg                                                                                               <  1:908125/150‑1 (MQ=255)
                                                    cAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACgcgcc                                                                                             >  2:311786/1‑149 (MQ=255)
                                                      gCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCgcagc                                                                                           >  2:409453/1‑150 (MQ=255)
                                                        ggTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATAGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCGCCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGcagcag                                                                                         >  2:751130/1‑150 (MQ=255)
                                                         gTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAgg                                                                                        <  1:1161468/150‑1 (MQ=255)
                                                         gTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGAGCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAgg                                                                                        <  1:1155407/150‑1 (MQ=255)
                                                           ggCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAATGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGGGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcg                                                                                      <  2:1541626/150‑1 (MQ=255)
                                                           ggCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACt                                                                                                            >  2:465473/1‑128 (MQ=255)
                                                             cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCCCCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                    >  1:1715024/1‑150 (MQ=255)
                                                             cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                    <  1:570445/150‑1 (MQ=255)
                                                             cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                    <  1:1200223/150‑1 (MQ=255)
                                                             cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                    <  1:956044/150‑1 (MQ=255)
                                                              gAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcgggc                                                                                   >  1:1629214/1‑150 (MQ=255)
                                                                     cccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACg                                                                            >  1:376040/1‑150 (MQ=255)
                                                                      ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGCACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACgg                                                                           >  1:854041/1‑150 (MQ=255)
                                                                         gCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCAt                                                                        >  1:538299/1‑150 (MQ=255)
                                                                            cGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgt                                                                     <  2:1133460/150‑1 (MQ=255)
                                                                               cgcCGATGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtggt                                                                  >  1:777971/1‑150 (MQ=255)
                                                                                  cGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGAt                                                               >  1:165807/1‑150 (MQ=255)
                                                                                  cGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGAt                                                               <  2:529689/150‑1 (MQ=255)
                                                                                  cGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGAt                                                               <  2:1629214/150‑1 (MQ=255)
                                                                                    aaGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCa                                                             <  2:1181899/150‑1 (MQ=255)
                                                                                    aaGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCa                                                             <  2:1626096/150‑1 (MQ=255)
                                                                                     agagCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATAGTGGTGATCaa                                                            <  2:918424/150‑1 (MQ=255)
                                                                                        gCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAccc                                                         <  1:1107451/150‑1 (MQ=255)
                                                                                        ccGTAAAGCCATTGACCGCTTTGCTATCGATCATGTCGATCGACGCGTGCCGGTGCTGATCTGCAACGGCGGAACCAACGCCCGCGTACCCATCGAACTCAGCCAGCCCGCGCAGCAGGCGGTCGAGGACGGCATCGAGGTGATCAAcac                                                         >  2:250184/2‑148 (MQ=255)
                                                                                           tAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGCCGGCATCGTGGTGATCAACCCcta                                                      >  1:434997/1‑150 (MQ=255)
                                                                                             aaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctact                                                    >  1:925334/1‑150 (MQ=255)
                                                                                                ccATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGTAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTactact                                                 <  2:1616784/150‑1 (MQ=255)
                                                                                                 cATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTg                                                >  1:566421/1‑150 (MQ=255)
                                                                                                 cATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGAACGGCACCGGCGGGCCCAACGCCCGGGAAACCAACGAAATAAGCCAGCACGCGCCGCCGGCGGGCGCGGACGGCCACGTGGTGATCCACCCcgacgtcag                                                >  2:1539027/1‑145 (MQ=255)
                                                                                                  aTTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTgg                                               <  2:1363979/150‑1 (MQ=255)
                                                                                                   ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGa                                              <  1:1159847/150‑1 (MQ=255)
                                                                                                   ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCATGGTGATCAACCCCTACTACTGGa                                              <  2:1052637/150‑1 (MQ=255)
                                                                                                    tGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaa                                             <  1:1344689/150‑1 (MQ=255)
                                                                                                     gCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAACCCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACAGGaaa                                            >  1:402919/1‑150 (MQ=255)
                                                                                                     gCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaaa                                            >  2:1513206/1‑150 (MQ=255)
                                                                                                     gCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCAGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaaa                                            >  1:844451/1‑150 (MQ=255)
                                                                                                      cccGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtg                                                                    >  2:59879/1‑125 (MQ=255)
                                                                                                        cGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAgtg                                         <  2:1125342/150‑1 (MQ=255)
                                                                                                             tGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGa                                    >  1:1383540/1‑150 (MQ=255)
                                                                                                             tGCTATCGAGCATGTCGATCGTCGCGGGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGTAAACCATCGTACTCAGCCAGCACGCGGAGCAGGCGGGCGCGGTCGGCATCGTGGTGATCAAGCCCTACTACTGGAAAGTGTCggt                                    <  2:969445/150‑2 (MQ=255)
                                                                                                                 aTCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGAGCAACCCCTACTACTGGAAAGTGTCGGAAGCg                                >  1:349394/1‑150 (MQ=255)
                                                                                                                   cGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCGCCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGaa                              <  1:751130/150‑1 (MQ=255)
                                                                                                                        aTGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAGCCTgc                         >  1:1269080/1‑149 (MQ=255)
                                                                                                                         tGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGAt                        >  2:675917/1‑150 (MQ=255)
                                                                                                                            cGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCCGCAGGCGGGCGCGGACGGCATCGTGGTGATCACCCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCg                     >  1:1547397/1‑150 (MQ=255)
                                                                                                                               tcgtcgCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTa                  >  1:533692/1‑150 (MQ=255)
                                                                                                                                cgtcgCGTGCCGGGGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTAt                 <  2:41180/150‑1 (MQ=255)
                                                                                                                                 gtcgCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTAtt                >  1:472487/1‑150 (MQ=255)
                                                                                                                                    gcgTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCg             <  1:429901/150‑1 (MQ=255)
                                                                                                                                     gccgCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGa            >  1:489005/4‑150 (MQ=255)
                                                                                                                                     gccgCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGa            >  2:833588/4‑150 (MQ=255)
                                                                                                                                     cgTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGa            <  2:1280086/150‑1 (MQ=255)
                                                                                                                                      ccgCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCCGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTCCTACTGGAAAGTGTCGGAAGCGCACCTGATCCGCTATTTCGAg           >  1:307774/3‑150 (MQ=255)
                                                                                                                                      ccgCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAg           >  2:1510665/3‑150 (MQ=255)
                                                                                                                                       tGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGc          >  1:1205982/1‑150 (MQ=255)
                                                                                                                                        gCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCa         <  1:619888/150‑1 (MQ=255)
                                                                                                                                           ggTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGt      >  1:498404/1‑150 (MQ=255)
                                                                                                                                            gTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTg     >  2:1584255/1‑150 (MQ=255)
                                                                                                                                               cTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGCACCTGATCCGCTATTTCGAGCAGGTGGcc  >  1:89888/1‑150 (MQ=255)
                                                                                                                                               cTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGcc  <  1:1703930/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCC  >  NC_000913/282369‑282661

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