Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,398,566 T→G 5.1% R385R (CGT→CGG mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,398,5660TG 5.1% 499.7 / 13.2 138R385R (CGT→CGGmutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (59/71);  new base G (2/5);  total (62/76)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 140 total aligned reads displayed.

CTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGT  >  NC_000913/4398423‑4398713
                                                                                                                                               |                                                                                                                                                   
cTACCGCTGGACGATGAACCCCATCCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGc                                                                                                                                                <  2:1150574/149‑1 (MQ=255)
  aCCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTg                                                                                                                                              <  2:1483175/149‑1 (MQ=255)
    cGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGCCGGACGACTGcc                                                                                                                                            <  2:218830/149‑1 (MQ=255)
    cGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGCCGGACGACTGcc                                                                                                                                            <  1:841567/149‑1 (MQ=255)
      cTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGCCGGACGACTGcctg                                                                                                                                          <  4:539777/149‑3 (MQ=255)
         gACGATGAGCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTg                                                                                                                                       >  3:1833568/1‑149 (MQ=255)
         gACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTg                                                                                                                                       >  1:1315281/1‑149 (MQ=255)
         gACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTg                                                                                                                                       >  1:1561926/1‑149 (MQ=255)
            gATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGcc                                                                                                                                    >  2:1636033/1‑149 (MQ=255)
            gATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGcc                                                                                                                                    >  1:1640290/1‑149 (MQ=255)
            gATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGcc                                                                                                                                    >  1:398137/1‑149 (MQ=255)
                aaCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAAt                                                                                                                                <  2:1315281/149‑1 (MQ=255)
                aaCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAAt                                                                                                                                <  2:398137/149‑1 (MQ=255)
                 aCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATg                                                                                                                               <  3:494137/149‑1 (MQ=255)
                  ccccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgc                                                                                                                              >  3:779808/1‑149 (MQ=255)
                   cccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGaa                                                                                                                                 <  3:579002/145‑1 (MQ=255)
                   cccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGaa                                                                                                                                 >  4:579002/1‑145 (MQ=255)
                    ccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgcgc                                                                                                                            <  4:1833568/149‑1 (MQ=255)
                    ccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgcgc                                                                                                                            >  2:105904/1‑149 (MQ=255)
                    ccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgcgc                                                                                                                            >  1:718344/1‑149 (MQ=255)
                      aaCCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAg                                                                                                                          >  3:763132/1‑149 (MQ=255)
                      aaCCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAg                                                                                                                          >  2:969426/1‑149 (MQ=255)
                        ccTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGcc                                                                                                                        >  4:716052/1‑149 (MQ=255)
                          tGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                      <  2:1561926/149‑1 (MQ=255)
                            cACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGc                                                                                                                    <  1:1636033/149‑1 (MQ=255)
                            cACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGc                                                                                                                    <  4:717755/149‑1 (MQ=255)
                            cACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGc                                                                                                                    <  2:718344/149‑1 (MQ=255)
                             aCCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAgggg                                                                                                                   <  1:464339/149‑3 (MQ=255)
                                  gTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGCCGGACGACTGcc                                                                                                                                            >  3:200091/1‑119 (MQ=255)
                                  gTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGCCGGACGACTGcc                                                                                                                                            <  4:200091/119‑1 (MQ=255)
                                     ccATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGaaa                                                                                                           <  1:436263/149‑1 (MQ=255)
                                       aTTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGCCGGACGACTGcctga                                                                                                                                         <  2:1136959/117‑4 (MQ=255)
                                        ttCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTAc                                                                                                                 <  3:1708528/140‑1 (MQ=255)
                                        ttCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTAc                                                                                                                 >  4:1708528/1‑140 (MQ=255)
                                         tCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGc                                                                                                       <  1:105904/149‑1 (MQ=255)
                                           cGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCaa                                                                                                     <  1:969426/149‑1 (MQ=255)
                                           cGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCaa                                                                                                     <  3:665543/149‑1 (MQ=255)
                                              aaaaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTAcc                                                                                                                >  2:872618/1‑135 (MQ=255)
                                              aaaaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTAcc                                                                                                                <  1:872618/135‑1 (MQ=255)
                                              aaaaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAg                                                                                                    >  4:360513/1‑147 (MQ=255)
                                              aaaaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAg                                                                                                    <  3:360513/147‑1 (MQ=255)
                                                aaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                      >  1:1006706/1‑127 (MQ=255)
                                                aaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                      <  2:1006706/127‑1 (MQ=255)
                                                  ccGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAg                                                                                              >  1:1722118/1‑149 (MQ=255)
                                                    gcgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                      <  1:239868/123‑1 (MQ=255)
                                                      gTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTa                                                                                          <  3:716052/149‑1 (MQ=255)
                                                        ggcggcGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATc                                                                                        <  4:1053449/149‑1 (MQ=255)
                                                             cggGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAg                                                                                   <  4:779808/149‑1 (MQ=255)
                                                                 gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAg                                                                                                        >  4:1307399/1‑124 (MQ=255)
                                                                 gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAg                                                                                                        <  3:1307399/124‑1 (MQ=255)
                                                                 gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAg                                                                                                    >  3:987801/1‑128 (MQ=255)
                                                                 gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGaa                                                                                               <  3:532189/133‑1 (MQ=255)
                                                                 gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCttt                                                                               >  2:202846/1‑149 (MQ=255)
                                                                  cgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                              <  2:271601/149‑1 (MQ=255)
                                                                     aaTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                      <  4:1603940/106‑1 (MQ=255)
                                                                     aaTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                      >  3:1603940/1‑106 (MQ=255)
                                                                       tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                              <  4:1766913/144‑1 (MQ=255)
                                                                       tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                              >  3:1766913/1‑144 (MQ=255)
                                                                       tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCaaa                                                                         >  2:776000/1‑149 (MQ=255)
                                                                       tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCaaa                                                                         <  3:1656938/149‑1 (MQ=255)
                                                                         aCTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATc                                                                                        <  2:978605/132‑1 (MQ=255)
                                                                          cTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACa                                                                                                         <  4:1142649/114‑1 (MQ=255)
                                                                          cTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACa                                                                                                         >  3:1142649/1‑114 (MQ=255)
                                                                           tttGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAAt                                                                                                                                >  1:1432615/1‑90 (MQ=255)
                                                                                  aaCCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCg                                                              >  1:433532/1‑149 (MQ=255)
                                                                                     cGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGaaa                                                                                                           >  3:224192/1‑101 (MQ=255)
                                                                                          gCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaa                                                      <  2:433532/149‑1 (MQ=255)
                                                                                           cTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaaa                                                     <  1:202846/149‑1 (MQ=255)
                                                                                            tCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAAt                                                    >  2:939008/1‑149 (MQ=255)
                                                                                              gTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTa                                                  >  4:231451/1‑149 (MQ=255)
                                                                                               tGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGAt                                                            <  1:352624/138‑1 (MQ=255)
                                                                                               tGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaa                                                 >  2:283047/1‑149 (MQ=255)
                                                                                                gAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaaa                                                <  1:776000/149‑1 (MQ=255)
                                                                                                    cGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgcg                                            <  1:876113/149‑1 (MQ=255)
                                                                                                        aGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCgg                                        <  1:283047/149‑1 (MQ=255)
                                                                                                          cTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAgg                                                                                                   >  2:1618359/1‑88 (MQ=255)
                                                                                                           tCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAAc                                     >  2:471897/1‑149 (MQ=255)
                                                                                                                  tACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGc                                                                            >  4:101974/1‑103 (MQ=255)
                                                                                                                  tACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGc                                                                            <  3:101974/103‑1 (MQ=255)
                                                                                                                      cTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCt                          >  4:1667598/1‑149 (MQ=255)
                                                                                                                       tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgc                                             >  4:373823/1‑129 (MQ=255)
                                                                                                                       tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgc                                             <  3:373823/129‑1 (MQ=255)
                                                                                                                         cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAgtgt                                                                                           <  2:128413/81‑1 (MQ=255)
                                                                                                                         cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAgtgt                                                                                           >  1:128413/1‑81 (MQ=255)
                                                                                                                         cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCaaa                                                                         >  2:369622/1‑99 (MQ=255)
                                                                                                                         cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCa                       <  4:8631/149‑1 (MQ=255)
                                                                                                                            cgcCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCg                                                              >  4:1077999/1‑107 (MQ=255)
                                                                                                                            cgcCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCg                                                              <  3:1077999/107‑1 (MQ=255)
                                                                                                                            cgcCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCg                                                              >  2:1044744/1‑107 (MQ=255)
                                                                                                                            cgcCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCg                                                              <  1:1044744/107‑1 (MQ=255)
                                                                                                                                  cATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaa                                                 <  2:949953/114‑1 (MQ=255)
                                                                                                                                  cATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaa                                                 >  1:949953/1‑114 (MQ=255)
                                                                                                                                   aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                              <  1:964441/84‑1 (MQ=255)
                                                                                                                                        gCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGt        >  2:634119/1‑149 (MQ=255)
                                                                                                                                        gCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGt        <  1:471897/149‑1 (MQ=255)
                                                                                                                                        gCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCACGAACCTGCACTTGCGGCGAGCAGt        <  1:1597575/149‑1 (MQ=255)
                                                                                                                                          aGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                              >  3:719008/1‑77 (MQ=255)
                                                                                                                                          aGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                              <  4:719008/77‑1 (MQ=255)
                                                                                                                                            tcgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCaga    <  2:1722118/149‑1 (MQ=255)
                                                                                                                                              gtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGt  >  2:157589/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGT  >  NC_000913/4398423‑4398713

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