Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,413 T→G 8.7% intergenic (‑12/+102) yihR ← / ← yihS putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4130TG 8.7% 299.0 / 21.4 92intergenic (‑12/+102)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (44/40);  new base G (4/4);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-02

GATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCC  >  NC_000913/4070143‑4070675
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                       
gATATTCCGGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     <  4:334757/275‑1 (MQ=255)
         gCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCGTAATTCACG‑ACATAAAGCCTCCCGTGAGGCGG‑ATGTc                                                                                                                                                                                                                                                            >  2:421451/1‑274 (MQ=255)
                                               ggCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                      <  2:151523/274‑1 (MQ=255)
                                                   aCGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                                                                                  <  4:53692/274‑1 (MQ=255)
                                                    cGGAATAAGCACTTTACCCAGATCCCCCAGCGGCATCTCTTCGGGTTTGTCCGGGATTACTAAATGACCCCCCTGAAAGGTTAATTCAGCCAGCCCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGAACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  1:75077/275‑1 (MQ=255)
                                                     ggAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgag                                                                                                                                                                                                               <  2:131988/275‑1 (MQ=255)
                                                                 ttACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGCTTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTT‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGCGAGAGGGTTAGGga                                                                                                                                                                                                    >  2:72682/1‑274 (MQ=255)
                                                                            cgcCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                                                                                                                                          >  1:154838/1‑209 (MQ=255)
                                                                            cgcCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  4:237616/1‑247 (MQ=255)
                                                                                aGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAccc                                                                                                                                                                                                                                                      <  3:154838/209‑1 (MQ=255)
                                                                                aGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  2:237617/247‑1 (MQ=255)
                                                                                                   tgtgCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                               >  4:538582/1‑181 (MQ=255)
                                                                                                     tgCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAACACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTa                                                                                                                                                                >  2:524348/1‑274 (MQ=255)
                                                                                                       cGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCa                                                                                                                                                                                                                                                           <  2:538583/181‑1 (MQ=255)
                                                                                                          gATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                                                                                                                                                          >  2:87625/1‑275 (MQ=255)
                                                                                                                        aaCCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  4:429662/1‑203 (MQ=255)
                                                                                                                          cccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCACTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGATAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                                                                                                                                           >  2:445036/1‑274 (MQ=255)
                                                                                                                            cTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  2:429663/203‑1 (MQ=255)
                                                                                                                             tGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGa                                                                                                                                                                                                                                                                       <  5:1076039/148‑1 (MQ=255)
                                                                                                                                    ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGc                                                                                                                                                                                                                                                                    <  8:110245/144‑1 (MQ=255)
                                                                                                                                    ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGc                                                                                                                                                                                                                                                                    >  7:110245/1‑144 (MQ=255)
                                                                                                                                     tAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                                                                                                                                                                                         >  2:314276/1‑217 (MQ=255)
                                                                                                                                       aTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAg                                                                                                                                                                                                                                                                      >  6:996749/1‑139 (MQ=255)
                                                                                                                                       aTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAg                                                                                                                                                                                                                                                                      <  5:996749/139‑1 (MQ=255)
                                                                                                                                             ccAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCCTATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGTTCCCCTCGCCCCTTTGGGGACAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCCGATCAGTCATTTCACATTAATATCCACCAGACCCGCTGCAACCCCTGGAg                                                                                                                       >  1:147008/1‑275 (MQ=255)
                                                                                                                                                   cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACAT‑CAGCCTCACGGGAGGCTgtatg                                                                                                                                                                                                                                                              >  5:490588/1‑130 (MQ=255)
                                                                                                                                                   cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACAT‑CAGCCTCACGGGAGGCTgtatg                                                                                                                                                                                                                                                              <  6:490588/134‑5 (MQ=255)
                                                                                                                                                    aGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACt                                                                                                                                                                                                                                               >  8:44140/1‑148 (MQ=255)
                                                                                                                                                           aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGCGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATTCAGAAGCACCGATGCACCCGCTGGCAGCATTCCTGgcgcc                                                                                                         >  2:537731/1‑274 (MQ=255)
                                                                                                                                                                 acaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTc                                                                                                                                                                                                                                  >  8:96848/1‑148 (MQ=255)
                                                                                                                                                                  caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaaa                                                                                                                                                                                                                                                         >  6:784275/1‑124 (MQ=255)
                                                                                                                                                                  caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaaa                                                                                                                                                                                                                                                         <  5:784275/124‑1 (MQ=255)
                                                                                                                                                                           cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                                                                                                                                          <  5:1126903/114‑1 (MQ=255)
                                                                                                                                                                           cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                                                                                                                                          >  6:1126903/1‑114 (MQ=255)
                                                                                                                                                                            gTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                                                                                                                                                                                       <  7:44140/148‑1 (MQ=255)
                                                                                                                                                                             ggATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                      <  7:96848/147‑1 (MQ=255)
                                                                                                                                                                                tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                               <  6:465230/104‑1 (MQ=255)
                                                                                                                                                                                tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                               >  5:465230/1‑104 (MQ=255)
                                                                                                                                                                                tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  6:923484/147‑1 (MQ=255)
                                                                                                                                                                                tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  5:923484/1‑147 (MQ=255)
                                                                                                                                                                                 aaTCCCCCGGGGCGATAGAACCGTACTGTCCAATACAATGCAAATGGGTGATATGCATCAGACCCACCTTAATAAACA‑ACAGACATCCTCCCGTGAGGCGT‑GTGTAAAACCCAAAAACTGCACCGACCGGTCCCCTCCCCCTTTTGGGGAGAGGGTTGGGGAGGGGGAAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCa                                                                                   <  4:445035/275‑1 (MQ=255)
                                                                                                                                                                                    ccccctggcaatgTGAAACCGTTCGTCCACTACAGTGACTTTTGGTGTTTTGCATACAATGCCCCTTAATTAACG‑CCTTACAGCCCCCCGTGAGGCTG‑TTGTCAAACCCAAAACCTGCCCCGGACGGCCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGGGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGATCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAgg                                                                                 <  4:87624/262‑1 (MQ=255)
                                                                                                                                                                                    cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  8:876529/143‑1 (MQ=255)
                                                                                                                                                                                    cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  7:876529/1‑143 (MQ=255)
                                                                                                                                                                                       ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCACGGGAGGCt                                                                                                                                                                                                                                                                   >  7:435844/1‑94 (MQ=255)
                                                                                                                                                                                       ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCACGGGAGGCt                                                                                                                                                                                                                                                                   <  8:435844/94‑1 (MQ=255)
                                                                                                                                                                                       ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCACGGGAGGCTgta                                                                                                                                                                                                                                                                <  6:239365/97‑3 (MQ=255)
                                                                                                                                                                                       ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCACGGGAGGCTgta                                                                                                                                                                                                                                                                >  5:239365/1‑95 (MQ=255)
                                                                                                                                                                                        cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                                                                                                                                                                                       <  5:914830/104‑1 (MQ=255)
                                                                                                                                                                                        cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                                                                                                                                                                                       >  6:914830/1‑104 (MQ=255)
                                                                                                                                                                                            tgagtagAAACGGTGTGCCCACAGAAAAGAAGAGTGGTTAGTTGGTTGCAAAGCACCTGAATTCAGGCACA‑GCAGCACCCCGTGGGCCGG‑AGGTGAAACCAAGAACCGGCGCCGGCCGGGCCCCGGGCCCCTTGAGGGGGAGGGTTAGGGGGACGGGAAGACCCGCACCCTATTTTCAGGTCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCATCCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcag                                                                         <  4:281740/268‑1 (MQ=255)
                                                                                                                                                                                              aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGCGAGAGGGTTAGGCAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCATCCCCCCTTCAACCGCTGGACCCATGCCGTGCGCTAACGGGATACCCGTGCTCCCCCGGCAAtcccg                                                                         >  4:250076/1‑268 (MQ=255)
                                                                                                                                                                                                           tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                                                                                  >  4:88582/1‑122 (MQ=255)
                                                                                                                                                                                                               cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  2:88583/122‑1 (MQ=255)
                                                                                                                                                                                                                 cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     <  5:563514/66‑1 (MQ=255)
                                                                                                                                                                                                                 cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     >  6:563514/1‑66 (MQ=255)
                                                                                                                                                                                                                 cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                                                                                                                                     >  8:191727/1‑129 (MQ=255)
                                                                                                                                                                                                                 cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                                                                                                                                     <  7:191727/129‑1 (MQ=255)
                                                                                                                                                                                                                       tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACt                                                                                                                                                                                                                                               >  1:27585/1‑81 (MQ=255)
                                                                                                                                                                                                                       tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  7:444499/108‑1 (MQ=255)
                                                                                                                                                                                                                       tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  8:444499/1‑108 (MQ=255)
                                                                                                                                                                                                                        gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATg                                                                                                                                                                                                                                                              <  6:1086383/65‑1 (MQ=255)
                                                                                                                                                                                                                        gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATg                                                                                                                                                                                                                                                              >  5:1086383/1‑65 (MQ=255)
                                                                                                                                                                                                                        gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAATCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                                                                                                                                           <  5:893759/116‑1 (MQ=255)
                                                                                                                                                                                                                        acATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                                                                                                                                           >  6:893759/2‑116 (MQ=255)
                                                                                                                                                                                                                         cATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  7:200494/106‑1 (MQ=255)
                                                                                                                                                                                                                         cATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  8:200494/1‑106 (MQ=255)
                                                                                                                                                                                                                           tatTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCAc                                                                                                                                                                                                                                           <  3:27585/81‑1 (MQ=255)
                                                                                                                                                                                                                                  gATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                                                                                                                                                   <  7:658136/98‑1 (MQ=255)
                                                                                                                                                                                                                                  gATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                                                                                                                                                   >  8:658136/1‑98 (MQ=255)
                                                                                                                                                                                                                                   aTTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑a                                                                                                                                                                                                                                                                >  6:885603/1‑52 (MQ=255)
                                                                                                                                                                                                                                   aTTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑a                                                                                                                                                                                                                                                                <  5:885603/52‑1 (MQ=255)
                                                                                                                                                                                                                                        cATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGACGGGAAAACCCGCACCCTATTTTCAGATCNGTGATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCGAACGGGAGACGCGGGNACACCAGGCAATGCACCAGGGGATAAAGATNNTGTTTCGTGTCCCACACTTTGgtggtg                              >  3:64913/1‑273 (MQ=255)
                                                                                                                                                                                                                                              aGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              >  7:140545/1‑91 (MQ=255)
                                                                                                                                                                                                                                              aGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  8:140545/91‑1 (MQ=255)
                                                                                                                                                                                                                                              aGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATNTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATACCCNGGCAATGCAGACGGGGATAAATATCCCGTTNNCCATCCCAGACTTTGGCGGTTActgtc                        >  3:58543/1‑269 (MQ=255)
                                                                                                                                                                                                                                                   ccTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              >  8:602876/1‑86 (MQ=255)
                                                                                                                                                                                                                                                   ccTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  7:602876/86‑1 (MQ=255)
                                                                                                                                                                                                                                                      taattaaCG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCa                                                                                                                                                                                                                                            <  6:258963/53‑1 (MQ=255)
                                                                                                                                                                                                                                                      taattaaCG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCa                                                                                                                                                                                                                                            >  5:258963/1‑53 (MQ=255)
                                                                                                                                                                                                                                                      taattaaCG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGccc                                                                                                                                                                                                                         >  7:202905/1‑72 (MQ=255)
                                                                                                                                                                                                                                                      taattaaCG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGccc                                                                                                                                                                                                                         <  8:202905/72‑1 (MQ=255)
                                                                                                                                                                                                                                                              g‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAGCAGTTATTTCGCATTAATATCCAGCAGATCCGCTGCACCCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATACTGTTTTCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCgatttcatctc       >  1:22945/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                   tACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGcc  >  2:523476/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     cAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                <  6:859602/83‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                >  5:859602/1‑83 (MQ=255)
                                                                                                                                                                                                                                                                       gCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtttt                                                                                                                                                                                                                    <  6:892725/62‑2 (MQ=25)
                                                                                                                                                                                                                                                                       gCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  5:892725/1‑62 (MQ=37)
                                                                                                                                                                                                                                                                             cGTGAGGCTG‑ATGTCAAACCCCAAAACTGCACGGGACGGTCCCCTCGCCCTTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAtt                                                                                                                                                              <  8:264309/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  6:372699/59‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                >  5:372699/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                             cGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAtt                                                                                                                                                              >  7:264309/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                       
GATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACG‑ACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCC  >  NC_000913/4070143‑4070675

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑