Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,398,580 G→A 11.5% W390* (TGG→TAG)  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,5800GA11.5% 252.3 / 18.9 87W390* (TGG→TAG) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (41/36);  new base A (5/5);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTA  >  NC_000913/4398437‑4398725
                                                                                                                                               |                                                                                                                                                 
tGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGa                                                                                                                                              >  1:189422/1‑149 (MQ=255)
 gAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGaa                                                                                                                                             <  2:671627/149‑1 (MQ=255)
  aaCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAAt                                                                                                                                            <  2:408652/149‑1 (MQ=255)
  aaCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAAt                                                                                                                                            <  2:840731/149‑1 (MQ=255)
  aaCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAAt                                                                                                                                            <  1:779526/149‑1 (MQ=255)
   aCCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATg                                                                                                                                           >  1:395770/1‑149 (MQ=255)
       cAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATgcgc                                                                                                                                        >  4:470093/1‑148 (MQ=255)
       cAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCa                                                                                                                                       <  2:395770/149‑1 (MQ=255)
        aaCCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAg                                                                                                                                      >  1:474036/1‑149 (MQ=255)
            tGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCa                                                                                                                                   >  4:88461/1‑148 (MQ=255)
                   cgTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTAcc                                                                                                                            <  3:871535/148‑1 (MQ=255)
                    gTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTGCCAg                                                                                                                          >  1:22202/1‑149 (MQ=255)
                    gTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAg                                                                                                                          >  2:825884/1‑149 (MQ=255)
                    gTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAg                                                                                                                          >  2:985092/1‑149 (MQ=255)
                      tCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGa                                                                                                                         >  3:881060/1‑148 (MQ=255)
                        cATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCAAATGCGCAGCCAGGCTACCAGAAAc                                                                                                                      <  2:889397/149‑1 (MQ=255)
                         aTTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACa                                                                                                                     >  1:839319/1‑149 (MQ=255)
                           tCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAg                                                                                                                    >  4:313444/1‑148 (MQ=255)
                            ccGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCa                                                                                                                  >  2:726823/1‑149 (MQ=255)
                             cGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCa                                                                                                                  <  3:313444/148‑1 (MQ=255)
                              ggAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAg                                                                                                                <  2:643180/149‑1 (MQ=255)
                              ggAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGACGAATGCGCAGCCAGGCTACCAGAAACAGCAAg                                                                                                                >  2:980841/1‑149 (MQ=255)
                               gAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAg                                                                                                                >  4:302834/1‑148 (MQ=255)
                                aaaaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCaa                                                                                                                 >  3:986064/1‑146 (MQ=255)
                                aaaaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCaa                                                                                                                 <  4:986064/146‑1 (MQ=255)
                                aaaaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGt                                                                                                              <  1:395875/149‑1 (MQ=255)
                                 aaaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGt                                                                                                              <  4:873763/148‑1 (MQ=255)
                                  aaCCGCGTGGCGGCGGGGCGCAATCACTTTGCAGATCCGGAAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTg                                                                                                             <  4:689295/148‑1 (MQ=255)
                                    ccGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGaa                                                                                                           <  3:470093/148‑1 (MQ=255)
                                     cgcgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAgt                                                                                                         <  2:839319/149‑1 (MQ=255)
                                       cgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAgtg                                                                                                        >  3:105326/1‑148 (MQ=255)
                                             ggcggGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGc                                                                                                  <  3:1053202/148‑1 (MQ=255)
                                                ggggCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGc                                                                                              >  2:633704/1‑149 (MQ=255)
                                                   gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCttt                                                                                           >  1:197460/1‑149 (MQ=255)
                                                    cgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                                          <  2:222840/149‑1 (MQ=255)
                                                     gcAATCACTTTGCAGAACCGGCCGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtttt                                                                                          <  3:649956/148‑1 (MQ=255)
                                                       aaTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCa                                                                                       >  1:91636/1‑149 (MQ=255)
                                                        aTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCa                                                                                       >  3:592423/1‑148 (MQ=255)
                                                         tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCaaa                                                                                     <  2:189422/149‑1 (MQ=255)
                                                              ttGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGcc                                                                                 >  3:458104/1‑148 (MQ=255)
                                                                   gAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCgcgc                                                                            >  4:627695/1‑148 (MQ=255)
                                                                    aaCCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGcc                                                                           >  3:227188/1‑148 (MQ=255)
                                                                      ccGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCAAGCATCAGGCAGTCGTCCGGATGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGAt                                                                        <  2:992098/149‑1 (MQ=255)
                                                                         gCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCa                                                                     >  1:496512/1‑149 (MQ=255)
                                                                          cAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCa                                                                     >  3:915242/1‑148 (MQ=255)
                                                                           aGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAGCAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaa                                                                   >  1:621659/1‑149 (MQ=255)
                                                                             cTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaaa                                                                 <  2:956979/149‑1 (MQ=255)
                                                                               cGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAAt                                                                >  4:296263/1‑148 (MQ=255)
                                                                                    gagtgTAGCGCGGAGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTCTCGCCa                                                                                                <  2:762486/107‑1 (MQ=255)
                                                                                    gCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAg                                                          <  1:633704/149‑1 (MQ=255)
                                                                                      cGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgcg                                                        >  2:156169/1‑149 (MQ=255)
                                                                                       ggTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgcg                                                        <  3:34229/148‑1 (MQ=255)
                                                                                        gTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgcgc                                                       <  4:139359/148‑1 (MQ=255)
                                                                                        gTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGcc                                                      <  1:114565/149‑1 (MQ=255)
                                                                                        gTAGCCCCCCGCTACATTCCTGCGCCAGCATAAGGCAGTCGTCCGACTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgcgct                                                      <  2:874467/149‑2 (MQ=255)
                                                                                         tAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGcc                                                      >  4:681491/1‑148 (MQ=255)
                                                                                             tCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAAc                                                 >  1:916294/1‑149 (MQ=255)
                                                                                              ccGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAAc                                                 >  3:721676/1‑148 (MQ=255)
                                                                                                      cacTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGaa                                         <  3:302834/148‑1 (MQ=255)
                                                                                                      cacTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCTGAACCGCAGGaa                                         <  4:915242/148‑1 (MQ=255)
                                                                                                       acTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAAc                                        <  3:88461/148‑1 (MQ=255)
                                                                                                         tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGc                                                                                  >  2:650693/1‑104 (MQ=255)
                                                                                                         tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGc                                                                                  <  1:650693/104‑1 (MQ=255)
                                                                                                          ccTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTg                                     >  4:308794/1‑148 (MQ=255)
                                                                                                           cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCa                                   >  2:380545/1‑149 (MQ=255)
                                                                                                               gcCAGCATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACtt                                <  3:308794/148‑1 (MQ=255)
                                                                                                                ccAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTg                               <  3:296263/148‑1 (MQ=255)
                                                                                                                 cAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcg                             >  2:348413/1‑149 (MQ=255)
                                                                                                                    cATCAGGCAGTCGTCCGGCTGCCCCCTAGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCTAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcggc                           >  4:762237/1‑148 (MQ=255)
                                                                                                                     aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcggcg                          >  4:367762/1‑148 (MQ=255)
                                                                                                                     aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgc                              <  1:513841/144‑1 (MQ=255)
                                                                                                                     aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgc                              >  2:513841/1‑144 (MQ=255)
                                                                                                                       cAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                             >  4:386274/1‑127 (MQ=255)
                                                                                                                       cAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                             <  3:386274/127‑1 (MQ=255)
                                                                                                                       cAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGaa                        >  3:989851/1‑148 (MQ=255)
                                                                                                                       cAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAAc                       <  1:726823/149‑1 (MQ=255)
                                                                                                                          gCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAg                     <  4:893023/148‑1 (MQ=255)
                                                                                                                              tcgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCaga                >  2:888004/1‑149 (MQ=255)
                                                                                                                               cgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCaga                <  3:316555/148‑1 (MQ=255)
                                                                                                                                 tcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGt              >  3:852796/1‑148 (MQ=255)
                                                                                                                                   cGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGttt            <  4:881060/148‑1 (MQ=255)
                                                                                                                                      cTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTgg         >  4:119756/1‑148 (MQ=255)
                                                                                                                                      cTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTgg         <  3:55632/148‑1 (MQ=255)
                                                                                                                                        gCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCg      >  1:381879/1‑149 (MQ=255)
                                                                                                                                            ccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCttt                                                                                           <  4:965175/60‑1 (MQ=255)
                                                                                                                                            ccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCttt                                                                                           >  3:965175/1‑60 (MQ=255)
                                                                                                                                             cTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa  >  3:8228/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTA  >  NC_000913/4398437‑4398725

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: ‑