Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,403,906 T→G 16.1% G202G (GGT→GGG hflC → HflB protease modulator specific for phage lambda cII repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,403,9060TG16.1% 170.0 / 2.9 75G202G (GGT→GGGhflCHflB protease modulator specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base T (3/53);  new base G (15/2);  total (20/55)
Fisher's exact test for biased strand distribution p-value = 6.90e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCGTGGTGAAGGCGATGCCGAAGCAG  >  NC_000913/4403670‑4404138
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
tttgcccttCTGCGTTCTTAAAGTGGTCGCCTGGGCGTGAAAGATTACGTCCCCGATTCCCGTGGTCGGCTGAACCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTAt                                                                                                                                                                                                                                          <  1:1296178/231‑1 (MQ=255)
    aCCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGCGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCCGGGGATTGa                                                                                                                                                                                                                                       >  2:400543/1‑238 (MQ=255)
       gTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGGAtgtatata                                                                                                                                                                                                                                   >  2:672844/1‑232 (MQ=255)
         cTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGGGTTGactgtct                                                                                                                                                                                                                                 >  2:849975/1‑233 (MQ=255)
          tGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCGGGGGATTGAAGTTGtg                                                                                                                                                                                                                               >  2:776689/1‑238 (MQ=255)
              ttgtGAAGTGGGTCGCCTGGGCGTGAAAGATATCGTGCCCGATTCCCGTGGTCGTCGGACCCGAGAAGGAGGGGACCGGCGGAACGCCGGTTCTGCGGGTAGAGAAGATGAAGGTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGAt                                                                                                                                                                                                                            <  1:1234348/235‑1 (MQ=255)
              ttCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGCAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGGGCCGGGAAtttggatgttgcgag                                                                                                                                                                                                                            >  2:623273/1‑226 (MQ=255)
               tCTGATATTGGGTGCCTTGCAGTGAAAAGTATCGGCAGCGATGCCCGGGGTCGGCTGACCCTCGAAGGACGGTGCGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGAtgt                                                                                                                                                                                                                          <  1:1048979/239‑1 (MQ=255)
                cTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGGGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGAtgt                                                                                                                                                                                                                          <  2:27591/238‑1 (MQ=255)
                    aaTTTGTCCCCTTGCCGTGAAAGATATCGGCCCCGATTCCCTTGGGCGTGTGACCCTGGAAGTACGTGACGCGCTGAACGCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTa                                                                                                                                                                                                                     <  1:611318/239‑1 (MQ=255)
                        ggTCGCCTGGACGTGAAAGATTTCGTCACCGATTCCCGGGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCaa                                                                                                                                                                                                                 <  2:146547/239‑1 (MQ=255)
                        ggTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGGGCGCGTAACGGCTGAGACGAAGGGAAAAGTTCCGGCCCTCACCCCGAAAAGTATGGCGGGCGCGGGGGTTGAAGAAGTCGAGGTGCgggtaaa                                                                                                                                                                                                                 >  1:1273011/1‑233 (MQ=255)
                         gTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGCCATCAACCCGAACAGATTGACGGCCCCGGGGATTGGAggttacgcgtgggggctcaca                                                                                                                                                                                                                >  2:735703/1‑219 (MQ=255)
                           cGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGACGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCGGGGGATTGTAGGTGGCGGTGTGCGGATCaaacc                                                                                                                                                                                                              >  2:1266467/1‑236 (MQ=255)
                            gCCTGGACGGTAAAGTTATCGTCACCGATTCCCGGTGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCa                                                                                                                                                                                                              <  2:221624/238‑1 (MQ=255)
                                       aaaGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCACTGGGCGTTGAAAGTGTTGGTTTGCGGAGCCAGCAAATCAAACttc                                                                                                                                                                                                   >  2:156451/1‑236 (MQ=255)
                                            tgtgctccacGATTCCCGTGTTCTTCTGTCTCGCCACGTTTGTGCCGCGCTGAGCCTCGGTTCTGCGGGTACAGAAGATGAAGTTGCTCCCCCGGCGGCAGAGAGCGCCATGGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCTGCGCTGGGGATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGccgac                                                                                                                                                                                               <  1:280936/228‑1 (MQ=255)
                                             atcggccacGATTCCCGTGGTAGTCTGACCCTCGAAGTACGGTACGCGCTGAACTCCGGTTCTGCGGGTAAAGAAGATGAAATTACTCCCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGccgaccga                                                                                                                                                                                            <  2:702404/231‑1 (MQ=255)
                                             accggccacGATTCCCGTGGTCGTCCGACCCTCTAAGTAGGTGCCGCGATGAACTCCGGTTCTGCGGGTTCAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGccgaccga                                                                                                                                                                                            <  2:772300/231‑1 (MQ=255)
                                                     cGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGGTGGGGATTGAGAGTTTCGATTTGCGGCACAAGCAGATTAACCTGCCGGCCGAAggggttg                                                                                                                                                                                    >  2:279621/1‑233 (MQ=255)
                                                        ttCCCGTGGTTCGTTGACCCGCGGAGGAGCTTGCGCGCGAAACTCGGGTTCTGCGGGTACAGAGGAATGAATTACTGCCCAGGCGGCGAACAGCGCCATTGCCGAAGTGGCAGCGCGAGTGACGGCTGAGGCGAAGGGCAAAGTTCCGATCATCAACCCGAACAGTATGGCGACGCTGGGTATTGAAGTTGTCTATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGGAg                                                                                                                                                                                 <  1:938793/239‑1 (MQ=255)
                                                        ttCCCGTGGTCTTTTCACCCCCGATGTACGTGACGCCCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTTCTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGGGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  1:1012538/239‑1 (MQ=255)
                                                        ttCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:1036802/239‑1 (MQ=255)
                                                        ttCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGGTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:1316457/239‑1 (MQ=255)
                                                        ttCCCGTGGTCGTCTGACCCGCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:960830/239‑1 (MQ=255)
                                                        ttCCCGTGGGCGTCTTACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:1059443/239‑1 (MQ=255)
                                                        ttCCCGTGGGCCGCTGCCCCCCGAATCCGATGACGCGCGTAAATCCGGTTCCGCGGGTCCAAAAGATGAAGGTGCCACCCCGGCGGCAGATAACCCCAATGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  1:372730/239‑1 (MQ=255)
                                                        ttCCCGGTGTAGTCCGACCCTCGAAGTACGGTACGCGCTGAACGCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:1325597/239‑1 (MQ=255)
                                                        tccccGTTGCCGTTTGACCCTGGAAGTACTTGACGCGCTGAACTCCGGTTCGGGGGGTAAAGAAGGTGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  1:1152397/237‑1 (MQ=255)
                                                         tCCCGTTGTCGTCTGACCCTCGAAGTCCGTGACGCGCTGAACTCCGGGTCTGCGGGTAAAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:557241/238‑1 (MQ=255)
                                                         tCCCGTGGTCGTCTGACCCTGGAAGTACGTGACGCGCTGAACGCCGGTTCGGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:571509/238‑1 (MQ=255)
                                                         tCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGGTGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:217880/238‑1 (MQ=255)
                                                         tCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:945741/238‑1 (MQ=255)
                                                         tCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:1273011/238‑1 (MQ=255)
                                                         tCCCGGGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTCCCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:485089/238‑1 (MQ=255)
                                                         tCCCGGGGTCGTCATCCCCTCGGAGTACGTGACGCGCTGAACCGCGGTTCTGCGGGTAAAGAAGGTGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  1:800402/238‑1 (MQ=255)
                                                          gcggTGGTCGTCTGACCCTCGAAGGACGGGACGGGCTGAACTCCGGTTCGGCGGGGACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCg                                                                                                                                                                               <  1:941/236‑1 (MQ=255)
                                                          cccGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCg                                                                                                                                                                               <  2:1303440/239‑1 (MQ=255)
                                                                  tcgtcTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTAc                                                                                                                                                                         <  2:1000883/238‑1 (MQ=255)
                                                                       tGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGGATTGAAGGTGGCGGTGGGCGGATCAAGCAGATCAACCTGCCGCCCGCAGTGTCTGGAGGGATCTACAAcca                                                                                                                                                                    >  2:941/1‑237 (MQ=255)
                                                                            cTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCCCTGGGGATTGAAGTTTTCGATGTGCGGTTCAAGCAGATCAACCTGCCCACCCCAGGGTCTTAAGGGATCTACAACCCTATgc                                                                                                                                                               >  2:544783/1‑238 (MQ=255)
                                                                                 aaatacttgacgccctggctgccGGTTGTGCGGGTACAGAAGGTGAAGTTGCTACCCCGGCGTCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCACAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGcc                                                                                                                                                          <  1:802689/218‑1 (MQ=255)
                                                                                   tACGTGACGCGCTGAACTCCGGTTCGGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAg                                                                                                                                                                                 <  2:1260084/213‑1 (MQ=255)
                                                                                         ggccccgtggtcttcggtttcgggggTACAGGAGATGAAGTTACTCCCCCGGCGGCAGATAACCCCATGGCCGAAGCGGCAGGGCTGATAACGGCGGGGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGGATGGCGGCGCTGGGTATGGAAGTTGGCGAGGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTg                                                                                                                                                   <  1:685541/215‑1 (MQ=255)
                                                                                                          ttCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGcgtcgt                                                                                                                                  <  2:785457/238‑1 (MQ=255)
                                                                                                           ttctgcggggtcaagaagggAAGGTACGCCCCCGGCGGGAGATAACGGCATTGGCGAAGCGGCAGGGCGAGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCAGCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGcgtcgt                                                                                                                                  <  1:143044/219‑1 (MQ=255)
                                                                                                           tatGCGGGGAAAGAAGATGAAGTTACTGCCCCGGCGGCAGATAACGCAGTTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCAGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCa                                                                                                                                <  2:549424/237‑1 (MQ=255)
                                                                                                                    aCAGAAGATGAGGTTACTCCCCCGGCGGGAGATAGCGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGAGCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTcaca                                                                                                                       <  1:68182/239‑1 (MQ=255)
                                                                                                                    aCAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCCCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTcaca                                                                                                                       <  2:695937/239‑1 (MQ=255)
                                                                                                                     aagaagaTGCAGTTACCACCCTGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGGAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGACCAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTcaca                                                                                                                       <  2:703374/237‑1 (MQ=255)
                                                                                                                      gcgagggggggttttggtccccTGCGGGAGATAACGCAATTGCCGAAGCGGCAGGGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGGCATCAACCCGATCAGTATGGCGGCGCTGGGGATTGAAGTTGTCGATGTGCGTATCAAGCAGATCATCCTGGCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTcaca                                                                                                                       <  1:169016/220‑1 (MQ=255)
                                                                                                                       gaagaTGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGGAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGGGGCTGAAGGGACCTACAACCGAATTCGGCCCGACGGTGAAGCGGGAGACGGTCGTCACCCTTccccacg                                                                                                                    >  2:113963/1‑233 (MQ=255)
                                                                                                                          aaatgaattaccccccccGGCGGCAGATAACGCCTCTGCCGAAGCGGAAAGGCGCGGAACGGCTGAGACGTAGGGGAAAGTTCCGGTCAACAACCCGAACAGTCCGGCGGCGCTGGGTATTGAAGTTGTCGATGGGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGAGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGt                                                                                                                   <  1:1226664/224‑1 (MQ=255)
                                                                                                                              gAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGGATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGAGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGAGTGAAACGGGAGCGCGTCGTGACCGTTCACAAGGTCAGga                                                                                                              >  2:390853/1‑239 (MQ=255)
                                                                                                                                             gcggcAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGGATTGTAGGTGTCGGCTGGCGTCTCAAGCAGATCAACCCGCCGCCCGACGTGGCTGCTGCGGTCTCTAACCCCAAGCGCGCCACGCGAGAAGAGCGGGAGCGCCGTCACCATCCACAAgcgcggcgagaaaagggaaaaaa                                                                                               >  1:1088414/1‑217 (MQ=255)
                                                                                                                                                 cgcatacggCCCTTGCCCAAGCGCGAGGGCGCGTAACGGCGTGGACGAAGGGCAATGTCCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATGGAAGTTGTCGATGTGCGTAGCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTgcg                                                                                           <  1:978876/230‑1 (MQ=255)
                                                                                                                                                         gCCGTTGCCTAAGCGGCAGAGCGGGTAACGGGTGAGACGATGGGCAAGGTGGCGGCACGCACACCGAACAGTATGGCGGCGCTGGCTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGGCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGc                                                                                  <  1:735703/239‑1 (MQ=255)
                                                                                                                                                             ttGCCGAAGGGGCAAGGCGCTGAACGGCTGATACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGAc                                                                              <  2:233351/239‑1 (MQ=255)
                                                                                                                                                             ttGCCGAAGCGGCAGAGCGCGTACCGGCTGATACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGAc                                                                              <  2:184313/239‑1 (MQ=255)
                                                                                                                                                              tGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGGATTTAAGTTGTCGATGTGCGTATCAAACAGATCAACCTGCCGACCGAAGGGTCTGAAGCGATCTACAACCGTATGCGCCCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCTGGAAGAAACGGGAAAACTGCCCCCGACTGCCgg                                                                               >  2:991013/1‑236 (MQ=255)
                                                                                                                                                              tGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCCGTGTGGCTTGAAGTTGTCGACGTGCGTTTCAAGCAGAACAACCTCCCGCCCGTAGAGTCTGAGGCGATCTACAACACTCTGGTGCCCCAGCGTGAAGAGGGAGGGGCGCGGCACCCGTCACAAGGTTACGAGAAACGCGGaaaaagacgccgcgcccgcgccc                                                                              >  1:1242877/1‑221 (MQ=255)
                                                                                                                                                              tGCCGAAGCGGCAGAGCGCGTAACGGCTGAAACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACt                                                                             <  2:1242877/239‑1 (MQ=255)
                                                                                                                                                                          agagCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAACAGATCAAACTGCCGACCGAAGTGTCTCAAGCGATCTACAAACCGATGGGCGCCCGGCGTGAAAGGGTAGCGGGGCGGCACCCTTCACAAGGGTAGGGAGAGACGGGAAAAATGCGCGCGACCGCCGACCCAGAAGGGACCa                                                                 >  1:332909/1‑239 (MQ=255)
                                                                                                                                                                              cgcgTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTcac                                                                                                                        <  2:783770/181‑1 (MQ=255)
                                                                                                                                                                               ccttagcGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCGGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAGCCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACg                                                            <  1:1301808/233‑1 (MQ=255)
                                                                                                                                                                                     cGGCTGAGACGAAGGGGAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGGTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCa                                                      <  2:897820/239‑1 (MQ=255)
                                                                                                                                                                                              cGAAGGGCAAAGTCCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTAGTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTg                                               <  2:784450/237‑1 (MQ=255)
                                                                                                                                                                                                 ggggCACTGTTCCGGACATTAACCCAACCAGGTGGGGGGCGCTGGCTATGTATGTTGTCGATGTGCGTATGAAGCAGATCAACCTGCCGACCGAAGTGTCTGACGCGCTCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGt                                          <  1:828417/238‑1 (MQ=255)
                                                                                                                                                                                                   ggCAAAGTTCCGGTCATAAACCCAAACGGTAGGGCGGCGCTGGTTATGTAAGTTGTGGGTGTGCGTAACAAGCACACAAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGt                                          <  1:1135909/237‑1 (MQ=255)
                                                                                                                                                                                                       aagggccagTCATCAACCCGAACAGGTTGGGGGCGCTGGGTATTGAAGTTTTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGc                                    <  1:16083/231‑1 (MQ=255)
                                                                                                                                                                                                             cGGTCATCAACCCGAACAGTATTGCGGCGCGGGGTATTGAAGTTGTCGATGGGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGcatc                              <  2:1215676/239‑1 (MQ=255)
                                                                                                                                                                                                                           aaCAGTATGGCGGCGCTGGGGATTTAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGGCAGGAAGAAGCGGGAAAACTGCGCGCGGCTGCCGACTATTAAGGGACCAGAAAGCTGGCAGAAGCTGAGCGTCAGGGGCGCATCATgg                          >  2:16083/1‑228 (MQ=255)
                                                                                                                                                                                                                                tATGGCGGCGATGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGc                                 <  2:594140/217‑1 (MQ=255)
                                                                                                                                                                                                                                         cgcTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGGGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGCCTGCCGGCTATGAAGTGACCCAAACGCTGGCAGAAGCTAAGCGCCAGGGCCCCACAATGCGTGGGGAAGGCGATGCCGaaca    >  2:847007/1‑236 (MQ=255)
                                                                                                                                                                                                                                           cTGGGTATTGAAGTTGTCGATGTGCGTAACAAGCAGATCAACCTGCCGACCGAAGTGTCGGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCGTGGTGAAGGCGATGCCGAagcag  <  2:864272/238‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
TCTGACCGTCTGCGTTCTGAAATTGGTCGCCTGGACGTGAAAGATATCGTCACCGATTCCCGTGGTCGTCTGACCCTCGAAGTACGTGACGCGCTGAACTCCGGTTCTGCGGGTACAGAAGATGAAGTTACTACCCCGGCGGCAGATAACGCCATTGCCGAAGCGGCAGAGCGCGTAACGGCTGAGACGAAGGGCAAAGTTCCGGTCATCAACCCGAACAGTATGGCGGCGCTGGGTATTGAAGTTGTCGATGTGCGTATCAAGCAGATCAACCTGCCGACCGAAGTGTCTGAAGCGATCTACAACCGTATGCGCGCCGAGCGTGAAGCGGTAGCGCGTCGTCACCGTTCACAAGGTCAGGAAGAAGCGGAAAAACTGCGCGCGACTGCCGACTATGAAGTGACCAGAACGCTGGCAGAAGCTGAGCGTCAGGGCCGCATCATGCGTGGTGAAGGCGATGCCGAAGCAG  >  NC_000913/4403670‑4404138

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