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

GATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCC  >  NC_000913/4070143‑4070675
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                       
gATATTCCGGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     <  4:334757/275‑1 (MQ=255)
         gCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCGTAATTCACGACATAAAGCCTCCCGTGAGGCGGATGTc                                                                                                                                                                                                                                                            >  2:421451/1‑274 (MQ=255)
                                               ggCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                      <  2:151523/274‑1 (MQ=255)
                                                   aCGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                                                                                  <  4:53692/274‑1 (MQ=255)
                                                    cGGAATAAGCACTTTACCCAGATCCCCCAGCGGCATCTCTTCGGGTTTGTCCGGGATTACTAAATGACCCCCCTGAAAGGTTAATTCAGCCAGCCCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGAACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  1:75077/275‑1 (MQ=255)
                                                     ggAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgag                                                                                                                                                                                                               <  2:131988/275‑1 (MQ=255)
                                                                 ttACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGCTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGCGAGAGGGTTAGGga                                                                                                                                                                                                    >  2:72682/1‑274 (MQ=255)
                                                                            cgcCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaa                                                                                                                                                                                                                                                          >  1:154838/1‑209 (MQ=255)
                                                                            cgcCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  4:237616/1‑247 (MQ=255)
                                                                                aGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAccc                                                                                                                                                                                                                                                      <  3:154838/209‑1 (MQ=255)
                                                                                aGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  2:237617/247‑1 (MQ=255)
                                                                                                   tgtgCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                                                                                                                                                                                               >  4:538582/1‑181 (MQ=255)
                                                                                                     tgCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAACACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTa                                                                                                                                                                >  2:524348/1‑274 (MQ=255)
                                                                                                       cGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCa                                                                                                                                                                                                                                                           <  2:538583/181‑1 (MQ=255)
                                                                                                          gATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGc                                                                                                                                                          >  2:87625/1‑275 (MQ=255)
                                                                                                                        aaCCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  4:429662/1‑203 (MQ=255)
                                                                                                                          cccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCACTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGATAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                                                                                                                                           >  2:445036/1‑274 (MQ=255)
                                                                                                                            cTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  2:429663/203‑1 (MQ=255)
                                                                                                                             tGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGa                                                                                                                                                                                                                                                                       <  5:1076039/148‑1 (MQ=255)
                                                                                                                                    ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGc                                                                                                                                                                                                                                                                    <  8:110245/144‑1 (MQ=255)
                                                                                                                                    ttAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGc                                                                                                                                                                                                                                                                    >  7:110245/1‑144 (MQ=255)
                                                                                                                                     tAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                                                                                                                                                                                         >  2:314276/1‑217 (MQ=255)
                                                                                                                                       aTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAg                                                                                                                                                                                                                                                                      >  6:996749/1‑139 (MQ=255)
                                                                                                                                       aTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAg                                                                                                                                                                                                                                                                      <  5:996749/139‑1 (MQ=255)
                                                                                                                                             ccAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCCTATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGTTCCCCTCGCCCCTTTGGGGACAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCCGATCAGTCATTTCACATTAATATCCACCAGACCCGCTGCAACCCCTGGAg                                                                                                                       >  1:147008/1‑275 (MQ=255)
                                                                                                                                                   cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATCAGCCTCACGGGAGGCTgtatg                                                                                                                                                                                                                                                              >  5:490588/1‑130 (MQ=255)
                                                                                                                                                   cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATCAGCCTCACGGGAGGCTgtatg                                                                                                                                                                                                                                                              <  6:490588/134‑5 (MQ=255)
                                                                                                                                                    aGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACt                                                                                                                                                                                                                                               >  8:44140/1‑148 (MQ=255)
                                                                                                                                                           aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGCGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATTCAGAAGCACCGATGCACCCGCTGGCAGCATTCCTGgcgcc                                                                                                         >  2:537731/1‑274 (MQ=255)
                                                                                                                                                                 acaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTc                                                                                                                                                                                                                                  >  8:96848/1‑148 (MQ=255)
                                                                                                                                                                  caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaaa                                                                                                                                                                                                                                                         >  6:784275/1‑124 (MQ=255)
                                                                                                                                                                  caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaaa                                                                                                                                                                                                                                                         <  5:784275/124‑1 (MQ=255)
                                                                                                                                                                           cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaa                                                                                                                                                                                                                                                          <  5:1126903/114‑1 (MQ=255)
                                                                                                                                                                           cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaa                                                                                                                                                                                                                                                          >  6:1126903/1‑114 (MQ=255)
                                                                                                                                                                            gTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                                                                                                                                                                                       <  7:44140/148‑1 (MQ=255)
                                                                                                                                                                             ggATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                      <  7:96848/147‑1 (MQ=255)
                                                                                                                                                                                tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                                                                                                                                                                                               <  6:465230/104‑1 (MQ=255)
                                                                                                                                                                                tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                                                                                                                                                                                               >  5:465230/1‑104 (MQ=255)
                                                                                                                                                                                tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  6:923484/147‑1 (MQ=255)
                                                                                                                                                                                tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  5:923484/1‑147 (MQ=255)
                                                                                                                                                                                 aaTCCCCCGGGGCGATAGAACCGTACTGTCCAATACAATGCAAATGGGTGATATGCATCAGACCCACCTTAATAAACAACAGACATCCTCCCGTGAGGCGTGTGTAAAACCCAAAAACTGCACCGACCGGTCCCCTCCCCCTTTTGGGGAGAGGGTTGGGGAGGGGGAAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCa                                                                                   <  4:445035/275‑1 (MQ=255)
                                                                                                                                                                                    ccccctggcaatgTGAAACCGTTCGTCCACTACAGTGACTTTTGGTGTTTTGCATACAATGCCCCTTAATTAACGCCTTACAGCCCCCCGTGAGGCTGTTGTCAAACCCAAAACCTGCCCCGGACGGCCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGGGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGATCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAgg                                                                                 <  4:87624/262‑1 (MQ=255)
                                                                                                                                                                                    cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  8:876529/143‑1 (MQ=255)
                                                                                                                                                                                    cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  7:876529/1‑143 (MQ=255)
                                                                                                                                                                                       ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCt                                                                                                                                                                                                                                                                   >  7:435844/1‑94 (MQ=255)
                                                                                                                                                                                       ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCt                                                                                                                                                                                                                                                                   <  8:435844/94‑1 (MQ=255)
                                                                                                                                                                                       ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgta                                                                                                                                                                                                                                                                <  6:239365/97‑3 (MQ=255)
                                                                                                                                                                                       ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgta                                                                                                                                                                                                                                                                >  5:239365/1‑95 (MQ=255)
                                                                                                                                                                                        cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                                                                                                                                                                                       <  5:914830/104‑1 (MQ=255)
                                                                                                                                                                                        cgcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                                                                                                                                                                                       >  6:914830/1‑104 (MQ=255)
                                                                                                                                                                                            tgagtagAAACGGTGTGCCCACAGAAAAGAAGAGTGGTTAGTTGGTTGCAAAGCACCTGAATTCAGGCACAGCAGCACCCCGTGGGCCGGAGGTGAAACCAAGAACCGGCGCCGGCCGGGCCCCGGGCCCCTTGAGGGGGAGGGTTAGGGGGACGGGAAGACCCGCACCCTATTTTCAGGTCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCATCCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcag                                                                         <  4:281740/268‑1 (MQ=255)
                                                                                                                                                                                              aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGCGAGAGGGTTAGGCAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCATCCCCCCTTCAACCGCTGGACCCATGCCGTGCGCTAACGGGATACCCGTGCTCCCCCGGCAAtcccg                                                                         >  4:250076/1‑268 (MQ=255)
                                                                                                                                                                                                           tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                                                                                  >  4:88582/1‑122 (MQ=255)
                                                                                                                                                                                                               cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  2:88583/122‑1 (MQ=255)
                                                                                                                                                                                                                 cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     <  5:563514/66‑1 (MQ=255)
                                                                                                                                                                                                                 cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAgg                                                                                                                                                                                                                                                                     >  6:563514/1‑66 (MQ=255)
                                                                                                                                                                                                                 cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                                                                                                                                     >  8:191727/1‑129 (MQ=255)
                                                                                                                                                                                                                 cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                                                                                                                                     <  7:191727/129‑1 (MQ=255)
                                                                                                                                                                                                                       tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACt                                                                                                                                                                                                                                               >  1:27585/1‑81 (MQ=255)
                                                                                                                                                                                                                       tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  7:444499/108‑1 (MQ=255)
                                                                                                                                                                                                                       tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  8:444499/1‑108 (MQ=255)
                                                                                                                                                                                                                        gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATg                                                                                                                                                                                                                                                              <  6:1086383/65‑1 (MQ=255)
                                                                                                                                                                                                                        gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATg                                                                                                                                                                                                                                                              >  5:1086383/1‑65 (MQ=255)
                                                                                                                                                                                                                        gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAATCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                                                                                                                                           <  5:893759/116‑1 (MQ=255)
                                                                                                                                                                                                                        acATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                                                                                                                                           >  6:893759/2‑116 (MQ=255)
                                                                                                                                                                                                                         cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    <  7:200494/106‑1 (MQ=255)
                                                                                                                                                                                                                         cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  8:200494/1‑106 (MQ=255)
                                                                                                                                                                                                                           tatTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCAc                                                                                                                                                                                                                                           <  3:27585/81‑1 (MQ=255)
                                                                                                                                                                                                                                  gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                                                                                                                                                   <  7:658136/98‑1 (MQ=255)
                                                                                                                                                                                                                                  gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                                                                                                                                                   >  8:658136/1‑98 (MQ=255)
                                                                                                                                                                                                                                   aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGa                                                                                                                                                                                                                                                                >  6:885603/1‑52 (MQ=255)
                                                                                                                                                                                                                                   aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGa                                                                                                                                                                                                                                                                <  5:885603/52‑1 (MQ=255)
                                                                                                                                                                                                                                        cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGACGGGAAAACCCGCACCCTATTTTCAGATCNGTGATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCGAACGGGAGACGCGGGNACACCAGGCAATGCACCAGGGGATAAAGATNNTGTTTCGTGTCCCACACTTTGgtggtg                              >  3:64913/1‑273 (MQ=255)
                                                                                                                                                                                                                                              aGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              >  7:140545/1‑91 (MQ=255)
                                                                                                                                                                                                                                              aGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  8:140545/91‑1 (MQ=255)
                                                                                                                                                                                                                                              aGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATNTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATACCCNGGCAATGCAGACGGGGATAAATATCCCGTTNNCCATCCCAGACTTTGGCGGTTActgtc                        >  3:58543/1‑269 (MQ=255)
                                                                                                                                                                                                                                                   ccTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              >  8:602876/1‑86 (MQ=255)
                                                                                                                                                                                                                                                   ccTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                              <  7:602876/86‑1 (MQ=255)
                                                                                                                                                                                                                                                      taattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCa                                                                                                                                                                                                                                            <  6:258963/53‑1 (MQ=255)
                                                                                                                                                                                                                                                      taattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCa                                                                                                                                                                                                                                            >  5:258963/1‑53 (MQ=255)
                                                                                                                                                                                                                                                      taattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGccc                                                                                                                                                                                                                         >  7:202905/1‑72 (MQ=255)
                                                                                                                                                                                                                                                      taattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGccc                                                                                                                                                                                                                         <  8:202905/72‑1 (MQ=255)
                                                                                                                                                                                                                                                              gACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAGCAGTTATTTCGCATTAATATCCAGCAGATCCGCTGCACCCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATACTGTTTTCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCgatttcatctc       >  1:22945/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                   tACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGcc  >  2:523476/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     cAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                <  6:859602/83‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                >  5:859602/1‑83 (MQ=255)
                                                                                                                                                                                                                                                                       gCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtttt                                                                                                                                                                                                                    <  6:892725/62‑2 (MQ=25)
                                                                                                                                                                                                                                                                       gCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                                                                                                                                    >  5:892725/1‑62 (MQ=37)
                                                                                                                                                                                                                                                                             cGTGAGGCTGATGTCAAACCCCAAAACTGCACGGGACGGTCCCCTCGCCCTTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAtt                                                                                                                                                              <  8:264309/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                <  6:372699/59‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                >  5:372699/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                             cGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAtt                                                                                                                                                              >  7:264309/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                       
GATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCC  >  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: