Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 21,750 A→T H115L (CAT→CTT)  ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091321,7500AT98.0% 148.6 / ‑6.2 50H115L (CAT→CTT) ribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (26/23);  minor base C (0/1);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTC  >  NC_000913/21505‑21994
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
cgcgcTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTgc                                                                                                                                                                                                                                                 <  1:34065/251‑1 (MQ=255)
  ggttGTTAAAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTTATGCTCTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTTCCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTgcgc                                                                                                                                                                                                                                               <  2:226190/248‑1 (MQ=255)
      gTTACAGGGCTTTCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTa                                                                                                                                                                                                                                            <  1:337459/250‑1 (MQ=255)
       ttACAGTGCTTGCAGGAAGAAAGGCGAAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGCACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGATTCGCGGAGTTAACCGCGCAAAATTTCCTCAGCGATCGTCTGGTGAAGCTTTTGCGCGTaa                                                                                                                                                                                                                                           >  2:5496/1‑250 (MQ=255)
          cAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAAtt                                                                                                                                                                                                                                       >  1:51371/1‑251 (MQ=255)
                 tGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGcc                                                                                                                                                                                                                                >  1:418352/1‑251 (MQ=255)
                  gCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCACGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGcc                                                                                                                                                                                                                                <  2:360689/250‑1 (MQ=255)
                       aagaagGGCGCAAGGTCAACTTACCGGTGATGGTGATGCTTCTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCCTTTGCGCGTAAAATTTCTTGCCGTAGGt                                                                                                                                                                                                                          <  1:256695/251‑1 (MQ=255)
                          aagGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGtgat                                                                                                                                                                                                                       <  2:3475/251‑1 (MQ=255)
                                     cgcaaCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGCGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCTGCtt                                                                                                                                                                                                            >  2:100152/1‑251 (MQ=255)
                                     cgcaaCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGATGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGCTCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGCTTTCCGCtt                                                                                                                                                                                                            >  1:99656/1‑251 (MQ=255)
                                      gcaaCTTACCGGTGATGGAGATGCATTTTGAACCTCAACCCCTGGAACTGTTTGCTCCCGATAAAGCCCCGGAACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCCGAGTGTGGCGTTGATTACGTGCTGTTCGTGCGTTTCGACAGGCGTTACGCGGCGTTAACCCCGCAAAATTTCATCAGCGATCTTCTGTTGAAGCTTTTGCGCGTAAAATTTCTTGCCGAAGGTGATGTTTTCCGCttt                                                                                                                                                                                                           >  2:473276/1‑251 (MQ=255)
                                                ggtgatggtgatGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTc                                                                                                                                                                                                 <  1:20350/251‑1 (MQ=255)
                                                           tgCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGAt                                                                                                                                                                                      >  1:184577/1‑251 (MQ=255)
                                                                          aaCCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGCAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTAc                                                                                                                                                                         >  2:166444/1‑249 (MQ=255)
                                                                            ccACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGa                                                                                                                                                                      >  1:339586/1‑250 (MQ=255)
                                                                                       gTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATgg                                                                                                                                                          >  1:148316/1‑251 (MQ=255)
                                                                                         ttGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGa                                                                                                                                                         >  2:459645/1‑250 (MQ=255)
                                                                                           gCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCGTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATa                                                                                                                                                      <  2:399210/251‑1 (MQ=255)
                                                                                             tACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGCTGATTTCCGCTTTGTCGCTGGACGAGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGACTACg                                                                                                                                                    >  2:136244/1‑251 (MQ=255)
                                                                                              aCCGATGAAGCCCCGGAACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTGGCAGAGTGGGGCGTTGATTACGTGCTGGTCGTGCGTGTCGACAGGCGTTTCGCGGCGTTACCCGCGCAAAATTTCATCAGCGATCTGCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATAc                                                                                                                                                     <  2:122324/249‑1 (MQ=255)
                                                                                                  aTAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGACGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTc                                                                                                                                               >  1:351369/1‑251 (MQ=255)
                                                                                                     aaGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGat                                                                                                                                            >  1:373676/1‑251 (MQ=255)
                                                                                                      aGCCCCGACACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGAGGTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGCAGCTTTTGCGCGGAAAATTTCTTGCCGAAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCGTGCATGGAATACGGCTTCGaa                                                                                                                                            >  1:302410/1‑249 (MQ=255)
                                                                                                                cggcTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGGGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTAc                                                                                                                                 <  2:51371/248‑1 (MQ=255)
                                                                                                                         cGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACttt                                                                                                                        <  1:216234/251‑1 (MQ=255)
                                                                                                                               cGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGCCAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGCTCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGACGGCGATTTCTTGTTATTCCAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAACCTTTTTGCg                                                                                                                   >  2:296916/1‑250 (MQ=255)
                                                                                                                                  gAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAgg                                                                                                               <  2:373676/251‑1 (MQ=255)
                                                                                                                                     aaaCTGCGTTACCTTGCAGAGTGTGGCGTTGATTGCGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTgg                                                                                                            <  1:10816/251‑1 (MQ=255)
                                                                                                                                               aCCTTGCAGAGTGGGGCGTTGATTTCGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTgcgtat                                                                                                   <  2:408275/250‑4 (MQ=255)
                                                                                                                                                cgctgccgagTGTGGCGTTGATTACGGGATGTGCGTTCGTTTCGACAGGCGTTTCGCGGCGGTTACCGCGCAAAATGTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGGTATGACAGAAAGCTGGCATGGAATACGGCTTCGAGATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATca                                                                                                 <  2:111922/244‑1 (MQ=255)
                                                                                                                                                ccTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGAtt                                                                                                                                                                                     >  2:390790/1‑167 (MQ=255)
                                                                                                                                                ccTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGGCGTGAAGGCGAtt                                                                                                                                                                                     <  1:390790/167‑1 (MQ=255)
                                                                                                                                                 cTTGCAGGGTGGGGCGCTGATTACGTGCTGTGCGTACGTTTCGCCAGGCGTTTCGTGGCGTTAACCGCGCATAATTTCACGAGTGGTCTTCTGGTGAAGCTTTTGCGCGTCAAAATTCTTGCTGTAGGCGATGCTTGCCGCTTTCGCGGTGGACGTGAAGGCGATTTCTTGCTATTACAGAAGGCTGGCATGGAATACGGCTTCGATATCACCAAACCGCAAACTTTTTGAGCAGGTGGCggccgtatc                                                                                                  >  2:288173/1‑241 (MQ=255)
                                                                                                                                                    gCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATcagc                                                                                               >  1:222747/1‑249 (MQ=255)
                                                                                                                                                          tgtgGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCAccgc                                                                                       >  1:62753/1‑251 (MQ=255)
                                                                                                                                                               cGTTGATTACGTGCTGGGCGGGCGTTTCGACAGGCTTTTCGCGGCGGTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGc                                                                                  <  2:125727/251‑1 (MQ=255)
                                                                                                                                                                           gCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGcgctgg                                                                      <  2:338276/251‑6 (MQ=255)
                                                                                                                                                                            cTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAGGGGGCCGTGCGTATCAGCGGCACCGCTGGGCGTCAGGcgctgcg                                                                     >  2:31845/1‑245 (MQ=255)
                                                                                                                                                                             tgtgCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGcgctggcg                                                                    >  2:1357/1‑244 (MQ=255)
                                                                                                                                                                                gcgtgcgtTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGAt                                                                 <  1:98245/251‑1 (MQ=255)
                                                                                                                                                                                               ggCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCgggtagatcgggagggcgt                                                  >  2:295660/1‑234 (MQ=255)
                                                                                                                                                                                               ggCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGAt                                                                 <  1:295660/236‑1 (MQ=255)
                                                                                                                                                                                                      cgcgGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCagaaa                                           <  1:203060/251‑3 (MQ=255)
                                                                                                                                                                                                        ctgcgTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCagaaa                                           <  2:222747/247‑3 (MQ=255)
                                                                                                                                                                                                           cgTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTa                                      <  1:1357/251‑1 (MQ=255)
                                                                                                                                                                                                             ttAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGAGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGCTGACAATCTGGCTCTGGCAGAAAGTTTACt                                    >  2:417646/1‑251 (MQ=255)
                                                                                                                                                                                                                aCCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGTCAATCTGGCTCTGGCAGAAAGTTTACTggg                                 >  2:76787/1‑251 (MQ=255)
                                                                                                                                                                                                                               agcaGCGATCTTCTGGTGAAGCTTTTGCGCGTACAATTTCTTGCCGTAGGTGATTATTTCCGCTTTGGCGCTGGTCGTGAATGCTATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTGCGATATCACCAGTACGCAAACTTTTTTCGAAGGTGGCGTGCGGATCAGCATCTCCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGTCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCAt                  <  1:39665/249‑1 (MQ=255)
                                                                                                                                                                                                                                               tGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTc  >  1:356996/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTC  >  NC_000913/21505‑21994

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