Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,398,563 T→C 5.1% S384S (AGT→AGC 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,5630TC 5.1% 482.1 / 10.3 138S384S (AGT→AGCmutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (60/71);  new base C (2/5);  total (62/76)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 147 total aligned reads displayed.

AACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGA  >  NC_000913/4398416‑4398711
                                                                                                                                                   |                                                                                                                                                    
aaCGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAgtc                                                                                                                                                     <  2:1864652/149‑1 (MQ=255)
       cTACCGCTGGACGATGAACCCCATCCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGc                                                                                                                                              <  2:1150574/149‑1 (MQ=255)
         aCCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTg                                                                                                                                            <  2:1483175/149‑1 (MQ=255)
           cGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGCCGGACGACTGcc                                                                                                                                          <  1:841567/149‑1 (MQ=255)
               ggACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAgtc                                                                                                                                                     >  2:786088/1‑134 (MQ=255)
               ggACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAgtc                                                                                                                                                     <  1:786088/134‑1 (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:398137/1‑149 (MQ=255)
                   gATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGcc                                                                                                                                  >  3:717755/1‑149 (MQ=255)
                   gATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGcc                                                                                                                                  >  1:1640290/1‑149 (MQ=255)
                        aCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATg                                                                                                                             <  3:494137/149‑1 (MQ=255)
                        aCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATg                                                                                                                             <  2:1200674/149‑1 (MQ=255)
                         ccccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgc                                                                                                                            >  3:779808/1‑149 (MQ=255)
                          cccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGaa                                                                                                                               >  4:579002/1‑145 (MQ=255)
                          cccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGaa                                                                                                                               <  4:1803131/145‑1 (MQ=255)
                           ccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgcgc                                                                                                                          >  2:105904/1‑149 (MQ=255)
                           ccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgcgc                                                                                                                          <  4:1833568/149‑1 (MQ=255)
                           ccAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgcgc                                                                                                                          >  1:718344/1‑149 (MQ=255)
                             aaCCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAg                                                                                                                        >  2:969426/1‑149 (MQ=255)
                                   cACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGc                                                                                                                  <  2:718344/149‑1 (MQ=255)
                                   cACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGc                                                                                                                  <  4:717755/149‑1 (MQ=255)
                                   cACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGc                                                                                                                  <  1:1636033/149‑1 (MQ=255)
                                    aCCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAgggg                                                                                                                 <  1:464339/149‑3 (MQ=255)
                                      cgcgTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAgt                                                                                                                                                      >  4:1532372/1‑110 (MQ=255)
                                         gTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGCCGGACGACTGcc                                                                                                                                          >  3:200091/1‑119 (MQ=255)
                                            ccATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGaaa                                                                                                         <  2:1047039/149‑1 (MQ=255)
                                               ttCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTAc                                                                                                               <  3:1708528/140‑1 (MQ=255)
                                               ttCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTAc                                                                                                               >  4:1708528/1‑140 (MQ=255)
                                                tCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGc                                                                                                     <  1:105904/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                                                                                                  <  3:360513/147‑1 (MQ=255)
                                                       aaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                    <  2:1006706/127‑1 (MQ=255)
                                                       aaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                    >  1:1006706/1‑127 (MQ=255)
                                                         ccGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAg                                                                                            >  1:1722118/1‑149 (MQ=255)
                                                           gcgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                    >  2:239868/1‑123 (MQ=255)
                                                                    cggGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAg                                                                                 <  4:779808/149‑1 (MQ=255)
                                                                        gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAg                                                                                                      >  4:1307399/1‑124 (MQ=255)
                                                                        gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAg                                                                                                      <  3:1307399/124‑1 (MQ=255)
                                                                        gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGaa                                                                                             <  3:532189/133‑1 (MQ=255)
                                                                        gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGaa                                                                                             >  4:532189/1‑133 (MQ=255)
                                                                        gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCttt                                                                             >  2:202846/1‑149 (MQ=255)
                                                                         cgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                            <  2:271601/149‑1 (MQ=255)
                                                                           cAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGc                                                                          <  2:1640290/149‑1 (MQ=255)
                                                                            aaTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                    <  4:1603940/106‑1 (MQ=255)
                                                                            aaTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAg                                                                                                                    >  3:1603940/1‑106 (MQ=255)
                                                                            aaTCACTTTGCAGAACCGGCAGCTCGGGAGCCGGTAGCTCCGCGCTACACTCCTGCACCAGCATCAGGCAGTCGACCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCa                                                                         >  3:1174653/1‑149 (MQ=255)
                                                                              tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgc                                                                                                                            >  3:1050521/1‑96 (MQ=255)
                                                                              tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgc                                                                                                                            <  4:1050521/96‑1 (MQ=255)
                                                                              tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                            >  3:1766913/1‑144 (MQ=255)
                                                                              tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCaaa                                                                       <  3:1656938/149‑1 (MQ=255)
                                                                                aCTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATc                                                                                      >  1:978605/1‑132 (MQ=255)
                                                                                aCTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATc                                                                                      <  2:978605/132‑1 (MQ=255)
                                                                                 cTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACa                                                                                                       <  4:1142649/114‑1 (MQ=255)
                                                                                  tttGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAAt                                                                                                                              >  1:1432615/1‑90 (MQ=255)
                                                                                  tttGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAAt                                                                                                                              <  2:1432615/90‑1 (MQ=255)
                                                                                        gAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAg                                                                                            <  2:711078/118‑1 (MQ=255)
                                                                                        gAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAg                                                                                            >  1:711078/1‑118 (MQ=255)
                                                                                        gAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGcc                                                             >  2:1597575/1‑149 (MQ=255)
                                                                                         aaCCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCg                                                            >  1:433532/1‑149 (MQ=255)
                                                                                            cGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGaaa                                                                                                         <  4:224192/101‑1 (MQ=255)
                                                                                            cGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGaaa                                                                                                         >  3:224192/1‑101 (MQ=255)
                                                                                                 gCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaa                                                    <  2:433532/149‑1 (MQ=255)
                                                                                                  cTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaaa                                                   <  1:202846/149‑1 (MQ=255)
                                                                                                   tCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAgtc                                                                                                                                                     <  1:1623679/50‑1 (MQ=255)
                                                                                                   tCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAAt                                                  >  3:8631/1‑149 (MQ=255)
                                                                                                   tCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAAt                                                  >  2:939008/1‑149 (MQ=255)
                                                                                                     gTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTa                                                >  4:231451/1‑149 (MQ=255)
                                                                                                      tGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGAt                                                          >  2:352624/1‑138 (MQ=255)
                                                                                                           cGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgcg                                          <  1:876113/149‑1 (MQ=255)
                                                                                                               aGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCgg                                      <  1:283047/149‑1 (MQ=255)
                                                                                                                 cTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAgg                                                                                                 >  2:1618359/1‑88 (MQ=255)
                                                                                                                 cTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAgg                                                                                                 <  1:1618359/88‑1 (MQ=255)
                                                                                                                      agcTACCCTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                               <  1:249433/148‑1 (MQ=255)
                                                                                                                         tACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGc                                                                          >  4:101974/1‑103 (MQ=255)
                                                                                                                             cTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCt                        >  4:1667598/1‑149 (MQ=255)
                                                                                                                              tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgc                                           <  3:373823/129‑1 (MQ=255)
                                                                                                                              tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgc                                           >  4:373823/1‑129 (MQ=255)
                                                                                                                                cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAgtgt                                                                                         <  2:128413/81‑1 (MQ=255)
                                                                                                                                cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAgtgt                                                                                         >  1:128413/1‑81 (MQ=255)
                                                                                                                                cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCaaa                                                                       <  1:369622/99‑1 (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                                                            <  1:1044744/107‑1 (MQ=255)
                                                                                                                                      cAGCATCAGGCAGTCGTCCGTCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaaa                                                   >  1:458715/1‑113 (MQ=255)
                                                                                                                                      cAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaaa                                                   <  2:458715/113‑1 (MQ=255)
                                                                                                                                         cATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaa                                               <  2:949953/114‑1 (MQ=255)
                                                                                                                                         cATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaa                                               >  1:949953/1‑114 (MQ=255)
                                                                                                                                          aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                            <  3:645770/84‑1 (MQ=255)
                                                                                                                                          aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                            <  1:964441/84‑1 (MQ=255)
                                                                                                                                          aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                            >  4:645770/1‑84 (MQ=255)
                                                                                                                                               gCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGt      >  2:634119/1‑149 (MQ=255)
                                                                                                                                                 aGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                            >  3:719008/1‑77 (MQ=255)
                                                                                                                                                 aGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                            <  4:719008/77‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGA  >  NC_000913/4398416‑4398711

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