Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 282,513 T→C 7.8% L79P (CTG→CCG)  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,5130TC 7.8% 349.0 / 17.1 129L79P (CTG→CCG) yagE2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (59/59);  new base C (5/5);  total (64/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-02

Only 100 of 130 total aligned reads displayed.

CGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCC  >  NC_000913/282364‑282661
                                                                                                                                                     |                                                                                                                                                    
cgccgcGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCt                                                                                                                                                      <  1:190357/150‑1 (MQ=255)
 gccgcGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTg                                                                                                                                                     >  1:1162120/1‑150 (MQ=255)
  ccgcGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGa                                                                                                                                                    <  2:1259419/150‑1 (MQ=255)
    gcgcTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATc                                                                                                                                                  >  1:1133460/1‑150 (MQ=255)
    gcgcTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATc                                                                                                                                                  <  1:543082/150‑1 (MQ=255)
     cgcTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGCTCg                                                                                                                                                 >  2:96643/1‑150 (MQ=255)
     cgcTGATCGACGATCTGATCAAAGAAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCCGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTAGATCGTCGCGTGCCGGTGCTGATCg                                                                                                                                                 >  1:1616784/1‑150 (MQ=255)
       cTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               <  2:1195390/150‑1 (MQ=255)
       cTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               >  1:1210358/1‑150 (MQ=255)
       cTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               >  1:1306528/1‑150 (MQ=255)
         gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCAc                                                                                                                                             >  2:429901/1‑150 (MQ=255)
          agcgacgaTCTGATCACCCCAGGCGGTGACGGCCTGTTCTTCCTTGGCAGCGGGGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAGCAcc                                                                                                                                            <  2:1290559/148‑1 (MQ=255)
            cgacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACcgg                                                                                                                                          >  2:823886/1‑150 (MQ=255)
             gacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGGGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAGCACcggc                                                                                                                                         <  1:63728/150‑1 (MQ=255)
             gacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAACGCCATTGCCCGCTTTGCTATCGATCATGTCGAACGTCGCGTGCCGGTGCTGATCGGCACcggc                                                                                                                                         >  2:106415/1‑150 (MQ=255)
             gacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAGCACcggc                                                                                                                                         <  2:488797/150‑1 (MQ=255)
                   ttgatgaATGCAGTCGTTGACGGCCTGTTATTCCTGGGGAGCGGGGGCGAGTTCTCCCAGCTCGGGGCCGAAGAGGGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCAcc                                                                                                                                   <  2:125969/149‑1 (MQ=255)
                      aTCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAAc                                                                                                                                <  2:623063/150‑1 (MQ=255)
                           aGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                               >  1:124348/1‑130 (MQ=255)
                           aGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCg                                                                                                                           <  1:1088719/150‑1 (MQ=255)
                           aGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCg                                                                                                                           >  1:1612546/1‑150 (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)
                                               tcttcCTTGGCAGCGGTGTCGAGTTCTCCCAGCTCGGCGCCGATGAGCGTTAAGCCATTGCCTGCTTTGCTATAGTTCTTGTCGATCGTCGCGGGCCGGTGCTGATCTGCAGCGGCGGCACCAACGCCCGGGAAACCATCGAACTcttcc                                                                                                       <  2:1514469/150‑5 (MQ=255)
                                               tcttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagcc                                                                                                       <  1:1701134/150‑1 (MQ=255)
                                                cttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagcca                                                                                                      >  1:653843/1‑150 (MQ=255)
                                                   ccTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCa                                                                                                   <  1:1671096/150‑1 (MQ=255)
                                                    cTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAc                                                                                                  <  1:813872/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)
                                                             ggTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATTGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGcagcag                                                                                         >  1:1667613/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                                                                                                            <  1:465473/128‑1 (MQ=255)
                                                                ggCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACt                                                                                                            >  2:465473/1‑128 (MQ=255)
                                                                  cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCCCCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                    >  1:1715024/1‑150 (MQ=255)
                                                                  cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                    >  1:919311/1‑150 (MQ=255)
                                                                  cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                    <  1:1200223/150‑1 (MQ=255)
                                                                  cGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                    <  1:570445/150‑1 (MQ=255)
                                                                   gAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCCGCCAGCACGCGCAGCAggcgggc                                                                                   >  1:229243/1‑150 (MQ=255)
                                                                      ttCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGcgcg                                                                                <  1:404313/150‑1 (MQ=255)
                                                                          cccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACg                                                                            >  1:376040/1‑150 (MQ=255)
                                                                           ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGCACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACgg                                                                           >  1:854041/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:1629214/150‑1 (MQ=255)
                                                                                         aaGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCa                                                             <  2:1181899/150‑1 (MQ=255)
                                                                                          agagCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATAGTGGTGATCaa                                                            <  2:918424/150‑1 (MQ=255)
                                                                                             gCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAccc                                                         <  1:1107451/150‑1 (MQ=255)
                                                                                             gCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAccc                                                         <  2:1210358/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)
                                                                                                     ccATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTactact                                                 <  1:511313/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)
                                                                                                        ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCATGGTGATCAACCCCTACTACTGGa                                              <  2:1052637/150‑1 (MQ=255)
                                                                                                         tGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaa                                             <  1:1344689/150‑1 (MQ=255)
                                                                                                          gCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaaa                                            >  2:1513206/1‑150 (MQ=255)
                                                                                                          gCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCAGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaaa                                            >  1:844451/1‑150 (MQ=255)
                                                                                                           cccGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtg                                                                    <  1:59879/125‑1 (MQ=255)
                                                                                                           cccGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtg                                                                    >  2:59879/1‑125 (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)
                                                                                                                                          cgTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGa            <  2:1280086/150‑1 (MQ=255)
                                                                                                                                           gTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAg           <  2:159664/150‑1 (MQ=255)
                                                                                                                                           ccgCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCCGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTCCTACTGGAAAGTGTCGGAAGCGCACCTGATCCGCTATTTCGAg           >  1:307774/3‑150 (MQ=255)
                                                                                                                                           ccgCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAg           >  1:1164661/3‑150 (MQ=255)
                                                                                                                                            tGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGc          >  1:1205982/1‑150 (MQ=255)
                                                                                                                                               cGGTGCGGAGCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGCACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAgg       <  2:1120726/150‑1 (MQ=255)
                                                                                                                                                ggTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGt      >  1:498404/1‑150 (MQ=255)
                                                                                                                                                 gTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTg     >  2:1584255/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCC  >  NC_000913/282364‑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: ‑