Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,927 C→A H250Q (CAC→CAA creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,9270CA95.2% 117.5 / ‑3.7 43H250Q (CAC→CAAcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (20/20);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTC  >  NC_000913/4638690‑4639172
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
cACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCt                                                                                                                                                                                                                                          <  1:103881/251‑1 (MQ=255)
  ccGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGATCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTgg                                                                                                                                                                                                                                        <  1:338871/251‑1 (MQ=255)
                   gggCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCgg                                                                                                                                                                                                                       >  1:98773/1‑251 (MQ=255)
                       agggcggtcATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAg                                                                                                                                                                                                                   <  2:285231/243‑1 (MQ=255)
                       aGGGCGTGCATATCACTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCAGCATCCAAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGAGTCATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTg                                                                                                                                                                                                                     >  1:118254/1‑249 (MQ=255)
                            gTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGttt                                                                                                                                                                                                              <  1:24101/251‑1 (MQ=255)
                             tGCATATCCCTTTACCCGAAGGGGACTGGCGGCAGCAGAACATGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGGGGTGCCTGCCGGGCGGACTAGCGAAATGGCCTTAACCAGCAACTGGGCGAATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCACACGGGAAGTTAGGGAGATCGGGTTTCAGGCGGAATGGGAAAGAAGCTGGTTTGCtattaaatccggtgggggtttt                                                                                                                                                                                                             >  2:58897/1‑235 (MQ=255)
                              gCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCCGCATATGGACCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTg                                                                                                                                                                                                            >  1:174892/1‑251 (MQ=255)
                                  atCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTc                                                                                                                                                                                                        >  1:390238/1‑251 (MQ=255)
                                            ccgAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATAAGGGCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCCGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGGGAGCGTTTTGCTTCAGGCAATGa                                                                                                                                                                                               >  1:377795/3‑250 (MQ=255)
                                             cgAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGAt                                                                                                                                                                                              >  1:283157/2‑250 (MQ=255)
                                              gAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGGATTTAAGCGGTACCGGCGATCTTTCTTTGGTGCCTGCCGGGCGGAATATCGAAATGACCTTAACCCGCAACTGGCCGCATCCCCGTTTTTTAGGAGATTTTTTGCCCGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGACGCAATAGCAACACCGCTGGTTTGCTAATAATCTCGGTGAGCGGTTTGCTTCAGGCAATgtga                                                                                                                                                                                             >  2:128913/1‑247 (MQ=255)
                                                   ggACTGGCGGAAGCGGAACTTGAAGCTGAATATGGCCCTGAAATTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCTAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCt                                                                                                                                                                                       <  2:51116/251‑1 (MQ=255)
                                                                 acaaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGTCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTcctc                                                                                                                                                                         <  2:200790/249‑3 (MQ=255)
                                                                 aGAACCTGAAGCTGAATATGGCCCTGAATTTACGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGAACTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCATAGCAGCTGGTTTGCTAATAATCTCGATGAGCGATTTGCTTCAGGCAATGATAACGGCTGGGAAAACtttccc                                                                                                                                                                         >  1:146923/1‑247 (MQ=255)
                                                                  gAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGGTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAAATTccccg                                                                                                                                                                        >  1:406393/1‑249 (MQ=255)
                                                                    aCCTGAAACTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAAATGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAATTTTGTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCcgcg                                                                                                                                                                      <  2:377795/251‑1 (MQ=255)
                                                                                          gAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCCGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAACGCAGCTGGTTTGCTAATAATCTCTGTGGGCGGTTTGCTTCAGGCAATGATCCCGGCTGGGAAAAATTACCCGCGTTTAAGGTCGCAgagacaaccg                                                                                                                                                >  1:269500/1‑242 (MQ=255)
                                                                                            aTTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTGCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGGTTAGCGTCGCAGTGACGACGCCa                                                                                                                                              <  1:378572/251‑1 (MQ=255)
                                                                                                aaGCGGTCCCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGcc                                                                                                                                           <  2:283157/250‑1 (MQ=255)
                                                                                                 aGCGGTCCCGGCGATCTTGCTGTGGTGCCGGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTGTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGa                                                                                                                                         <  2:445126/251‑1 (MQ=255)
                                                                                                     gttccGGCGATCTTTCTGGGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGGCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATTCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCaa                                                                                                                                     <  1:418483/248‑1 (MQ=255)
                                                                                                                  ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTTGCGTCGCCGTGACGCCGGCAGCCGGTCAAttcacattaa                                                                                                                           >  2:165486/1‑239 (MQ=255)
                                                                                                                  ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTaa                                                                                                                           <  1:165486/248‑1 (MQ=255)
                                                                                                                       gggggCCTGCCGGGCGTAAAAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGATTCAGGTTTTCAGGCGCAATGGAAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGGGACGACGCCAGCCGAGCAATACCAATTAAcggtccgg                                                                                                                   <  2:367211/247‑8 (MQ=255)
                                                                                                                                   gggggaaaaGCGAAATGACCCTAACCAGCAACTGGGCGAATCCCAGGTTTTTAGGTGAATTTTTGCTAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAACGTCAAAGCAGCTGGTTTGCTGATAACCTCGGTGAGCGTTTTGTTTCGGGCAATGATACGGGCTGGGAAAACGTCCCCGCGTTTAGCGGCGCAGTGACGAGGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTAc                                                                                                       <  2:45351/243‑1 (MQ=255)
                                                                                                                                              gAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCg                                                                                                                                                                                          >  2:70724/1‑157 (MQ=255)
                                                                                                                                              gAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCg                                                                                                                                                                                          <  1:70724/157‑1 (MQ=255)
                                                                                                                                                           ccAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGa                                                                                 >  1:25410/1‑249 (MQ=255)
                                                                                                                                                           ccAGCAACTGGCCGCAGCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGa                                                                                 <  2:25410/249‑1 (MQ=255)
                                                                                                                                                                       cGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGt                                                                                                         >  1:355533/1‑213 (MQ=255)
                                                                                                                                                                       cGCATCCCAGTTTTTTAGGTGAGTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTGCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATGCCAATGAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTc                                                                   <  1:223283/251‑1 (MQ=255)
                                                                                                                                                                        gCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTAAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGCCTTTTATGGCGTTCt                                                                  >  1:476684/1‑251 (MQ=255)
                                                                                                                                                                            cccAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGt                                                              >  1:495646/1‑251 (MQ=255)
                                                                                                                                                                                       aGGTGATTTTTTGCTAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAAGGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTATCGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACg                                                    <  1:465693/250‑1 (MQ=255)
                                                                                                                                                                                        ggTGATTTTTTGCCAGCCACACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTGGCGTAGCCGTGGCGACGCCAGCCGATCCATACCAAATCACGGAACGGGCGAATTAATACCACACTTTGCTGATTTCACTGGCTTTTCTGGCGTTCTTTTTTTTTGAAAcacc                                                  >  2:61603/1‑248 (MQ=255)
                                                                                                                                                                                         gtgagttgttgccAGCCAACCGGGAAGTTAGCGAGTCAGGTTTGCCGGCGCAAGGGCAAAGCAGCTGGTTGGTTAATAATTGCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTGGAAACGCTc                                                 <  2:149178/240‑1 (MQ=255)
                                                                                                                                                                                                  tgccAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGGCTAAGTACGCCCTTCTGCTGATTTCCATGGCTTTTATGGCGTTCTTTTTTTTTGAAACGCTCACCTCGCa                                         >  2:229749/3‑250 (MQ=255)
                                                                                                                                                                                                      aGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                                    <  1:211303/251‑1 (MQ=255)
                                                                                                                                                                                                                gAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAAt                          >  1:40980/1‑251 (MQ=255)
                                                                                                                                                                                                                     tAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAat                     >  1:331507/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tCAGGCGTAATGGCAAAGCAGCTGGTTTGCTAATAATTTCGGTGGGCGTTTTGCTACAGGCAAGGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTgggg       <  2:274168/250‑1 (MQ=255)
                                                                                                                                                                                                                                          cgcAATGGCAAAGCAGCTGGTTTGCTACTAATCTCGGTGAGCGTTTTGCTTCAGGCACTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCCTTCTGCTGATTGCACTGCCTTTTATGGCGTTCTTTGTTTTTTAAACGCTCACCGCGCAACGTTTAAACCCAATGCAATATTTGCTGGTGGGGCTTTc  >  1:163111/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
CACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTC  >  NC_000913/4638690‑4639172

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