Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 282,513 T→C 5.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 5.8% 344.5 / 8.3 103L79P (CTG→CCG) yagE2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (50/47);  new base C (3/3);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 103 total aligned reads displayed.

CGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGG  >  NC_000913/282369‑282659
                                                                                                                                                |                                                                                                                                                  
cgcTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATc                                                                                                                                                <  3:709093/149‑1 (MQ=255)
   tGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGc                                                                                                                                             <  3:149571/149‑1 (MQ=255)
     aTCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCAc                                                                                                                                           >  2:460309/1‑149 (MQ=255)
        gacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACcgg                                                                                                                                        >  2:961203/1‑149 (MQ=255)
        gacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACcgg                                                                                                                                        >  1:269047/1‑149 (MQ=255)
            aTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCggcggc                                                                                                                                    >  2:1138315/1‑149 (MQ=255)
                 aTCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCaa                                                                                                                               >  2:739471/1‑149 (MQ=255)
                  tCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAAc                                                                                                                              <  2:1128247/149‑1 (MQ=255)
                           gCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGaa                                                                                                                     >  3:764039/1‑149 (MQ=255)
                             gTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAAc                                                                                                                   >  4:279096/1‑149 (MQ=255)
                              ttGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAAcc                                                                                                                  <  1:713325/149‑1 (MQ=255)
                                gACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCAt                                                                                                                <  1:961203/149‑1 (MQ=255)
                                   ggCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGa                                                                                                             >  1:243292/1‑149 (MQ=255)
                                    gCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCg                                                                                                                         <  2:564803/136‑1 (MQ=255)
                                    gCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCg                                                                                                                         >  1:564803/1‑136 (MQ=255)
                                    gCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGaa                                                                                                            <  3:354096/149‑1 (MQ=255)
                                    gCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGaa                                                                                                            <  4:525919/149‑1 (MQ=255)
                                      cTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACt                                                                                                          >  2:956677/1‑149 (MQ=255)
                                         ttcttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcag                                                                                                       >  2:65256/1‑149 (MQ=255)
                                         ttcttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcag                                                                                                       >  2:186961/1‑149 (MQ=255)
                                            ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGc                                                                                                                            >  1:1148140/1‑125 (MQ=255)
                                            ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGc                                                                                                                            <  2:1148140/125‑1 (MQ=255)
                                            ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagcca                                                                                                    <  1:739471/149‑1 (MQ=255)
                                            ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagcca                                                                                                    >  4:1052978/1‑149 (MQ=255)
                                             tcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagccag                                                                                                   >  2:412226/1‑149 (MQ=255)
                                               cTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCa                                                                                                 <  4:725565/149‑1 (MQ=255)
                                                 gggCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCAcg                                                                                               >  4:479663/1‑149 (MQ=255)
                                                              gAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcggg                                                                                  >  3:500325/1‑149 (MQ=255)
                                                                 ttCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGATTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGgcgc                                                                               >  4:238976/1‑149 (MQ=255)
                                                                    tcCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGa                                                                            >  2:1073007/1‑149 (MQ=255)
                                                                      ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCg                                                                                                              >  4:1139887/1‑113 (MQ=255)
                                                                      ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCg                                                                                                              <  3:1139887/113‑1 (MQ=255)
                                                                      ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagc                                                                                                      <  4:1123141/121‑1 (MQ=255)
                                                                      ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagc                                                                                                      >  3:1123141/1‑121 (MQ=255)
                                                                      ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGCGCCGGTGCCGATCAGCACcggca                                                                                                                                      <  4:473364/89‑2 (MQ=39)
                                                                      ccAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGCGCCGGTGCCGATCAGCACcggca                                                                                                                                      >  3:473364/1‑88 (MQ=39)
                                                                        aGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggc                                                                                     >  4:380948/1‑136 (MQ=255)
                                                                        aGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggc                                                                                     <  3:380948/136‑1 (MQ=255)
                                                                        aGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGc                                                                        >  2:451871/1‑149 (MQ=255)
                                                                           tCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGaa                                                                                                            >  3:114693/1‑110 (MQ=255)
                                                                           tCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGaa                                                                                                            <  4:114693/110‑1 (MQ=255)
                                                                             ggCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGTACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgt                                                                   >  4:669281/1‑149 (MQ=255)
                                                                              gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtg                                                                  >  3:579545/1‑149 (MQ=255)
                                                                              gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtg                                                                  <  4:281535/149‑1 (MQ=255)
                                                                               cgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                                                                 <  1:451871/149‑1 (MQ=255)
                                                                                gcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtggt                                                                <  2:532924/149‑1 (MQ=255)
                                                                                  cGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGa                                                              <  4:79206/149‑1 (MQ=255)
                                                                                       agCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAc                                                         <  4:500325/149‑1 (MQ=255)
                                                                                           tAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCct                                                     >  2:375671/1‑149 (MQ=255)
                                                                                           tAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCct                                                     >  1:97837/1‑149 (MQ=255)
                                                                                             aaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctac                                                   <  3:1113200/149‑1 (MQ=255)
                                                                                              aGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCa                                                                       >  2:871510/1‑128 (MQ=255)
                                                                                              aGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCa                                                                       <  1:871510/128‑1 (MQ=255)
                                                                                               gCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCaa                                                                                                                               <  3:557222/71‑1 (MQ=255)
                                                                                               gCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCaa                                                                                                                               >  4:557222/1‑71 (MQ=255)
                                                                                                ccATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCtactac                                                <  4:764039/149‑1 (MQ=255)
                                                                                                ccATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCtactac                                                <  2:606588/149‑1 (MQ=255)
                                                                                                   ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTgg                                             >  4:344161/1‑149 (MQ=255)
                                                                                                   ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTgg                                             >  1:1126898/1‑149 (MQ=255)
                                                                                                   ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTgg                                             >  4:246029/1‑149 (MQ=255)
                                                                                                        cGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAgt                                        <  1:375671/149‑1 (MQ=255)
                                                                                                           tttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTc                                     <  4:569281/149‑1 (MQ=255)
                                                                                                           tttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTc                                     <  1:186961/149‑1 (MQ=255)
                                                                                                           tttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTc                                     >  2:577378/1‑149 (MQ=255)
                                                                                                           tttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTc                                     <  1:412226/149‑1 (MQ=255)
                                                                                                           tttGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTc                                     <  2:97837/149‑1 (MQ=255)
                                                                                                               cTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGaa                                 <  1:460309/149‑1 (MQ=255)
                                                                                                               cTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGaa                                 <  3:723885/149‑1 (MQ=255)
                                                                                                                tATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAg                                >  3:467322/1‑149 (MQ=255)
                                                                                                                tATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAg                                >  2:817413/1‑149 (MQ=255)
                                                                                                                 aTCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGc                               >  4:844180/1‑149 (MQ=255)
                                                                                                                  tCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCgg                                                                                                                        <  2:1092350/59‑1 (MQ=255)
                                                                                                                  tCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCgg                                                                                                                        >  1:1092350/1‑59 (MQ=255)
                                                                                                                  tCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCg                              <  1:65256/149‑1 (MQ=255)
                                                                                                                    gATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGaa                            >  4:379270/1‑149 (MQ=255)
                                                                                                                    gATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGaa                            >  4:699661/1‑149 (MQ=255)
                                                                                                                    gATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGaa                            >  2:573021/1‑149 (MQ=255)
                                                                                                                      tCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAAcc                          <  1:817413/149‑1 (MQ=255)
                                                                                                                       cATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCt                         >  2:945830/1‑149 (MQ=255)
                                                                                                                         tGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGa                       >  4:56837/1‑149 (MQ=255)
                                                                                                                         tGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGa                       <  3:1052978/149‑1 (MQ=255)
                                                                                                                          gTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGAt                      <  1:1058482/149‑1 (MQ=255)
                                                                                                                          gTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGAt                      >  2:1058482/1‑149 (MQ=255)
                                                                                                                          gTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGAt                      >  1:274445/1‑149 (MQ=255)
                                                                                                                            cGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATcc                    <  2:1010430/149‑1 (MQ=255)
                                                                                                                            cGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATcc                    <  4:637378/149‑1 (MQ=255)
                                                                                                                                  tcgCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTAtt              <  1:1073007/149‑1 (MQ=255)
                                                                                                                                    gcgTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTc            <  3:479663/149‑1 (MQ=255)
                                                                                                                                      gTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAAc                           <  4:810978/132‑1 (MQ=255)
                                                                                                                                      gTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAAc                           >  3:810978/1‑132 (MQ=255)
                                                                                                                                      gTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGa          <  3:844180/149‑1 (MQ=255)
                                                                                                                                       tGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAg         <  3:56837/149‑1 (MQ=255)
                                                                                                                                        gCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                                                                 >  3:917824/1‑92 (MQ=255)
                                                                                                                                        gCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTc                                     <  1:1119457/120‑1 (MQ=255)
                                                                                                                                        gCCGGTGCCGATCAGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTc                                     >  2:1119457/1‑120 (MQ=255)
                                                                                                                                            gTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGt    >  1:1051363/1‑149 (MQ=255)
                                                                                                                                             tGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAAcc                          >  1:481946/1‑126 (MQ=255)
                                                                                                                                             tGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAAcc                          <  2:481946/126‑1 (MQ=255)
                                                                                                                                             tGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTg   <  3:238976/149‑1 (MQ=255)
                                                                                                                                              gCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTgg  >  1:668932/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGG  >  NC_000913/282369‑282659

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