Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 788,385 G→T 5.2% T151T (ACC→ACA) ‡ galM ← aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913788,3850GT 5.2% 523.2 / 11.7 150T151T (ACC→ACA) ‡galMaldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base G (72/70);  new base T (4/4);  total (76/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

Only 100 of 156 total aligned reads displayed.

TCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCAT  >  NC_000913/788114‑788658
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                 
tCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGCTTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTAtcgtcgtt                                                                                                                                                                                                                                                                                 >  2:221440/1‑272 (MQ=255)
   cTTTGGGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATGTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                                                                                                                                             <  4:219833/275‑1 (MQ=255)
                                tCACTGGGAATGATGTTGGTGCTGCGGAAATCGAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTAATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTAc                                                                                                                                                                                                                                                 <  2:415912/274‑1 (MQ=255)
                                cgcgtgggcatgatTTTGGTGCGTGGGAAATAAAAAGCGTTGTCGGTCGCAGATTTAGGGCCTACGAGGGGAATGCCGCCTGCACCACCCGGGAGATATTCGTCCGCCAGAATCTGCAAGTTGGAATTGCCAACGTCAGACTGCTCGCCGTCAAGATGGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGAAACGCGATAAGTAGTGGAGGTACGGTGATGGGCGGTCAGACGATATTGCACCGACGCGCCGGGATTAc                                                                                                                                                                                                                                                 <  4:7237/265‑1 (MQ=255)
                                                tGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGCTTACCCGGGAAGCCCTgctg                                                                                                                                                                                                                                 >  1:9905/1‑271 (MQ=255)
                                                           aaaTAAACAGCCGTTCCGGGGNNCACATNTCCGGCCGTCGTCCGTAATGCCGCCTTCAGCACACGTCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGGCAGACTGCTCGCCGTCAAGAGTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCt                                                                                                                                                                                                                       <  3:9905/274‑1 (MQ=255)
                                                                         tCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACa                                                                                                                                                                                                       >  2:536135/1‑275 (MQ=255)
                                                                                            ggCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGCCGATCGtt                                                                                                                                                                                       >  2:261234/1‑274 (MQ=255)
                                                                                                           tGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACCATCTGc                                                                                                                                                                       >  1:49416/1‑275 (MQ=255)
                                                                                                                       caACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                                                                                                                                             <  1:138768/161‑1 (MQ=255)
                                                                                                                            ggCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGaa                                                                                                                                                      <  2:287485/275‑1 (MQ=255)
                                                                                                                                        tCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCa                                                                                                                                                                                                                                                                               <  8:268085/142‑1 (MQ=255)
                                                                                                                                        tCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCa                                                                                                                                                                                                                                                                               >  7:268085/1‑142 (MQ=255)
                                                                                                                                                           tgtgATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCCAACAGCACCTGACGATCGTTCTTGTTCACAATCTTCCAGCGACGTTTGTCGAACCCTTCAGGCCCGCCGTGCAGCTGGTTAACg                                                                                                                       >  2:486701/1‑275 (MQ=255)
                                                                                                                                                                ttGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                             >  8:803923/1‑136 (MQ=255)
                                                                                                                                                                ttGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                             <  7:803923/136‑1 (MQ=255)
                                                                                                                                                                    gcACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGatat                                                                                                                                                                                                                                                                       >  5:139219/1‑122 (MQ=255)
                                                                                                                                                                    gcACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGatat                                                                                                                                                                                                                                                                       <  6:139219/122‑1 (MQ=255)
                                                                                                                                                                     cACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAg                                                                                                                                                                                                                                                                              >  7:729073/1‑114 (MQ=255)
                                                                                                                                                                     cACGTAAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAg                                                                                                                                                                                                                                                                              >  6:492230/1‑114 (MQ=255)
                                                                                                                                                                         tCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGc                                                                                                                                                                                                                                        >  8:506607/1‑148 (MQ=255)
                                                                                                                                                                          cAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGaa                                                                                                        <  1:9076/275‑1 (MQ=255)
                                                                                                                                                                             aCTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                             <  8:513558/123‑1 (MQ=255)
                                                                                                                                                                             aCTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                             >  7:513558/1‑123 (MQ=255)
                                                                                                                                                                                   cGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTAc                                                                                                                                                                                                                                                 >  6:1102431/1‑129 (MQ=255)
                                                                                                                                                                                    gCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAGGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTCGCGAAAGCGTCCCgg                                                                                              >  2:91525/1‑275 (MQ=255)
                                                                                                                                                                                    gCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACgg                                                                                              >  4:412424/1‑275 (MQ=255)
                                                                                                                                                                                        naaaGATNANAATAGACGTGATTAGTCATATTCACCGGGCAAGGNTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGc                                                                                                                                                                   <  4:554643/200‑1 (MQ=255)
                                                                                                                                                                                         cAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcatc                                                                                                                                                                                                                           <  8:386543/145‑1 (MQ=255)
                                                                                                                                                                                         cAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcatc                                                                                                                                                                                                                           >  7:386543/1‑145 (MQ=255)
                                                                                                                                                                                           aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGataa                                                                                                                                                                                                                                                                       >  8:625661/1‑98 (MQ=255)
                                                                                                                                                                                           aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                    >  6:746327/1‑134 (MQ=255)
                                                                                                                                                                                           aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                    <  5:746327/134‑1 (MQ=255)
                                                                                                                                                                                              ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCg                                                                                                                                                                                                                                                        >  5:787790/1‑111 (MQ=255)
                                                                                                                                                                                              ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAAc                                                                                                                                                                                                                   >  5:1019730/1‑148 (MQ=255)
                                                                                                                                                                                              ttGAAATAGACGTAATTAGTCATATTCACCGGGCAAGTTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAAc                                                                                                                                                                                                                   >  7:77863/1‑148 (MQ=255)
                                                                                                                                                                                              tggAAAAAGCCGTGATTAGTGAAGTTCAACGGCACAAGTTTATGAACTGGGGCGGGAGAAGTAATAGAGAGTCAGTTATCGGCGGTCAGACGACAATGCACAGTGGGGCCGGGATTAACCGGTAAGCCCTGATCACCATCATCTGAACGCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTGCACCGTc                                                                                    <  4:361930/273‑1 (MQ=255)
                                                                                                                                                                                                taaaGAGCCGTGATTAGTCATATTCACCGGGCAAGGGTTATCAACTGTGGCGCGATAAGTAATGTGGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCaa                                                                                  <  4:79822/274‑1 (MQ=255)
                                                                                                                                                                                                gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                 <  5:959719/148‑1 (MQ=255)
                                                                                                                                                                                                gAAATAGACGTGATTAGTCATATTAACCGGGCAAGGTTTAGCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGGGATTACCCGGGAAGCCCTGATCACCATCATCTGAACGCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTGCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCaa                                                                                  <  1:115203/275‑1 (MQ=255)
                                                                                                                                                                                                    tAGACTTGATTAGTAATATTCACCGGGCAAGGTTTATCAACTGTGGGGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCAGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAgg                                                                               <  2:149437/274‑1 (MQ=255)
                                                                                                                                                                                                      gACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAAGGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCa                                                                                                                                                                                                           >  6:916314/1‑148 (MQ=255)
                                                                                                                                                                                                            aTTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCCCAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGt                                                                                                                                     >  4:329274/1‑212 (MQ=255)
                                                                                                                                                                                                                gTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAg                                                                                                                                 <  2:329275/212‑1 (MQ=255)
                                                                                                                                                                                                                           cGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGTTTAACGCCCTGACTTTGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCCATTGGCGATActg                                                        >  2:202325/1‑272 (MQ=255)
                                                                                                                                                                                                                             ggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTg                                                                                                                                                                                    <  7:553802/148‑1 (MQ=255)
                                                                                                                                                                                                                             ggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACCGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCCGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTg                                                                                                                                                                                    <  5:916314/148‑1 (MQ=255)
                                                                                                                                                                                                                               cAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                 <  6:339578/117‑1 (MQ=255)
                                                                                                                                                                                                                               cAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                 >  5:339578/1‑117 (MQ=255)
                                                                                                                                                                                                                                 aGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                                                                                                                                                <  8:212379/100‑1 (MQ=255)
                                                                                                                                                                                                                                   gTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATGTTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                                                                                                                                                >  4:538817/1‑98 (MQ=255)
                                                                                                                                                                                                                                   gTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTcac                                                                                                                                                                              >  6:963320/1‑148 (MQ=255)
                                                                                                                                                                                                                                    tttATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGccc                                                                                                                                                                                                                                      >  5:270951/1‑91 (MQ=255)
                                                                                                                                                                                                                                    tttATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGccc                                                                                                                                                                                                                                      <  6:270951/91‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGt                                                                                                                                                                                                                                                                                 <  4:526267/45‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATGTTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcat                                                                                                                                                                                                                            <  2:538818/98‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGc                                                                                                                                                                                                            >  3:29619/1‑114 (MQ=255)
                                                                                                                                                                                                                                       aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCAc                                                                                        >  4:220411/1‑230 (MQ=255)
                                                                                                                                                                                                                                        tCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCt                                                                                                                                                                         >  5:464063/1‑148 (MQ=255)
                                                                                                                                                                                                                                         cAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTg                                                                                                                                                                        >  6:411154/1‑148 (MQ=255)
                                                                                                                                                                                                                                          aaCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACt                                                                                                                                                                                                                  >  6:432566/1‑105 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCATACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcat                                                                                                                                                                                                                            <  2:560551/94‑1 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                 >  6:221957/1‑105 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                 <  5:221957/105‑1 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAAc                                                                                                                                                                                                        <  1:29619/114‑1 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGAGCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTc                                                                                    <  2:220412/230‑1 (MQ=255)
                                                                                                                                                                                                                                            cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTg                                                                                                                                                                                                                      <  5:1088476/99‑1 (MQ=255)
                                                                                                                                                                                                                                            cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTg                                                                                                                                                                                                                      >  6:1088476/1‑99 (MQ=255)
                                                                                                                                                                                                                                            cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGa                                                                                                                                                                                                                     >  8:682794/1‑100 (MQ=255)
                                                                                                                                                                                                                                            cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGa                                                                                                                                                                                                                     <  7:682794/100‑1 (MQ=255)
                                                                                                                                                                                                                                               tgGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATagag                                   <  2:465890/275‑1 (MQ=255)
                                                                                                                                                                                                                                                ggCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGa                                                                                                                                                                 <  6:464063/148‑1 (MQ=255)
                                                                                                                                                                                                                                                  cgcgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACg                                                                                                                                                               >  7:215417/1‑148 (MQ=255)
                                                                                                                                                                                                                                                   gcgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGt                                                                                                                                                              >  8:1116101/1‑148 (MQ=255)
                                                                                                                                                                                                                                                    cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCt                                                                                                                                                                                     <  7:1120775/124‑1 (MQ=255)
                                                                                                                                                                                                                                                    cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCt                                                                                                                                                                                     >  8:1120775/1‑124 (MQ=255)
                                                                                                                                                                                                                                                    cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGtt                                                                                                                                                             <  8:215417/148‑1 (MQ=255)
                                                                                                                                                                                                                                                       tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAc                                                                                                                                                                                                                               >  7:1015982/1‑79 (MQ=255)
                                                                                                                                                                                                                                                       tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAc                                                                                                                                                                                                                               <  8:1015982/79‑1 (MQ=255)
                                                                                                                                                                                                                                                       tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACg                                                                                                                                                                                             <  6:302662/113‑1 (MQ=255)
                                                                                                                                                                                                                                                       tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACg                                                                                                                                                                                             >  5:302662/1‑113 (MQ=255)
                                                                                                                                                                                                                                                       tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTTTGGTACACAATCTGCCAGCGACGTTTGt                                                                                                                                                          <  8:77863/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                          <  7:307424/127‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                          >  8:307424/1‑127 (MQ=255)
                                                                                                                                                                                                                                                             aTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAAc                                                                                                                                                                                                        >  7:195510/1‑96 (MQ=255)
                                                                                                                                                                                                                                                              tGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGcc                                                                                                                                                                      >  6:796162/1‑129 (MQ=255)
                                                                                                                                                                                                                                                               ggAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                    >  7:447382/1‑66 (MQ=255)
                                                                                                                                                                                                                                                               ggAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCt                                                                                                                                                  >  5:138326/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                   aTACGGTTATCGGCGGTCCAACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGTTCACCATCATCTGAACTCAGGTCAAACAGCACCTGACGATCGTGCTGGTTCACAAGCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcct                <  4:486700/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGc                                                                                                                                            <  6:138326/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCCTAACGACCAATCGAGGCCCCCAGAAACTCGGcctga              >  2:138027/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                     aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCACTAGAGGCCCCCAGAAACGCGGcctgat             >  2:6236/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAAAGCGGcctga              >  2:356751/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                     aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGAACACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTTGAACCTTTCCGGCCCGCCGTGCCGCTGGTTAACGCCCTGACTTGGCGAAAGCGGCACGGTTTCACCGTCAAATGTCTAACGGCTATTGGCGATACGGTTGGCag                                                 >  2:243104/1‑238 (MQ=255)
                                                                                                                                                                                                                                                                         ttATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCaa                                                                                                                                                                                                          >  5:429089/1‑82 (MQ=255)
                                                                                                                                                                                                                                                                          tATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGt                                                                                                                                                                                        >  8:758060/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                           aTCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGt                                                                                                                                                                                  <  5:488751/104‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                          >  5:53567/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                               tcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                          <  6:53567/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 ggTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCAt  >  2:378922/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                 
TCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCAT  >  NC_000913/788114‑788658

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