Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,885 A→G L236M (CTA→TTG)  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,8850AG93.3% 114.5 / ‑2.2 45L236M (CTA→TTG) creDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (20/22);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

GTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCT  >  NC_000913/4638640‑4639124
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
gtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGc                                                                                                                                                                                                                                              >  1:210857/1‑251 (MQ=255)
gtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGCGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAAAACCTGAAGCAGAAAATGGACCTGAATTTAAGCGGTCCCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTCAAAGCGAAATGACCTTAACCAGGAACTGGCCGCATGCCAGGTTTTTAGGGGATTTTTTGCCAGc                                                                                                                                                                                                                                              >  1:22394/1‑251 (MQ=255)
               gAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCCAACGGGAAGtt                                                                                                                                                                                                                                 >  1:285231/1‑249 (MQ=255)
                    tAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAg                                                                                                                                                                                                                           <  1:319578/250‑1 (MQ=255)
                      acggaactgCGCCGACCATTGAACCCGGCCTCGGGTTTGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGTGACTGGCGTAATCAGAACCTGAAGCTGAATATGGCCCGGAATTTAAGCGGTGCCGTCGATCTTTCTGTGGGGCCTGCCGGGCGTAATGGCGAAAGGACCTTAACCAGCAACTGGGCGCATCCCAGTTTTTTGGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAggca                                                                                                                                                                                                                        <  2:284386/251‑4 (MQ=255)
                             ggCGCTGACAATTGAACCCGGCACCGGGATAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCTAAATGACCTTAACCAGCAACTGGCCGAATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTc                                                                                                                                                                                                                <  2:34728/251‑1 (MQ=255)
                               cgcTGACCATTGAACCCGGCACCGGGCTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTc                                                                                                                                                                                                                >  1:33660/1‑249 (MQ=255)
                                gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAgg                                                                                                                                                                                                             <  1:288207/251‑1 (MQ=255)
                                   gACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGgc                                                                                                                                                                                                            <  1:4354/249‑1 (MQ=255)
                                                  cACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCt                                                                                                                                                                                          <  1:103881/251‑1 (MQ=255)
                                                    ccGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGATCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTgg                                                                                                                                                                                        <  1:338871/251‑1 (MQ=255)
                                                                     gggCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCgg                                                                                                                                                                       >  1:98773/1‑251 (MQ=255)
                                                                         agggcggtcATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAg                                                                                                                                                                   <  2:285231/243‑1 (MQ=255)
                                                                         aGGGCGTGCATATCACTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCAGCATCCAAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGAGTCA‑TGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTg                                                                                                                                                                     >  1:118254/1‑249 (MQ=255)
                                                                              gTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGttt                                                                                                                                                              <  1:24101/251‑1 (MQ=255)
                                                                               tGCATATCCCTTTACCCGAAGGGGACTGGCGGCAGCAGAACATGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGGGGTGCCTGCCGGGCGGACTAGCGAAATGGCCTTAACCAGCAACTGGGCGAATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCACACGGGAAGTTAGGGAGATC‑GGGTTTCAGGCG‑GAATGGGAAAGAAGCTGGTTTGCtattaaatccggtgggggtttt                                                                                                                                                             >  2:58897/1‑235 (MQ=255)
                                                                                gCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCCGCATATGGACCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTg                                                                                                                                                            >  1:174892/1‑251 (MQ=255)
                                                                                    atCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTc                                                                                                                                                        >  1:390238/1‑251 (MQ=255)
                                                                                              ccgAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATAAGGGCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCCGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGGGAGCGTTTTGCTTCAGGCAATGa                                                                                                                                               >  1:377795/3‑250 (MQ=255)
                                                                                               cgAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGAt                                                                                                                                              >  1:283157/2‑250 (MQ=255)
                                                                                                gAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGGATTTAAGCGGTACCGGCGATCTTTCTTTGGTGCCTGCCGGGCGGAATATCGAAATGACCTTAACCCGCAACTGGCCGCATCCCCGTTTTTTAGGAGATTTTTTGCCCGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGACG‑CAATAGCAACACCGCTGGTTTGCTAATAATCTCGGTGAGCGGTTTGCTTCAGGCAATgtga                                                                                                                                             >  2:128913/1‑247 (MQ=255)
                                                                                                     ggACTGGCGGAAGCGGAACTTGAAGCTGAATATGGCCCTGAAATTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCTAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCt                                                                                                                                       <  2:51116/251‑1 (MQ=255)
                                                                                                                   acaaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGTCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTcctc                                                                                                                         <  2:200790/249‑3 (MQ=255)
                                                                                                                   aGAACCTGAAGCTGAATATGGCCCTGAATTTACGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGAACTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCATAGCAGCTGGTTTGCTAATAATCTCGATGAGCGATTTGCTTCAGGCAATGATAACGGCTGGGAAAACtttccc                                                                                                                         >  1:146923/1‑247 (MQ=255)
                                                                                                                    gAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGGTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAAATTccccg                                                                                                                        >  1:406393/1‑249 (MQ=255)
                                                                                                                      aCCTGAAACTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAAATGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAATTTTGTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCcgcg                                                                                                                      <  2:377795/251‑1 (MQ=255)
                                                                                                                                            gAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCCGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAACGCAGCTGGTTTGCTAATAATCTCTGTGGGCGGTTTGCTTCAGGCAATGATCCCGGCTGGGAAAAATTACCCGCGTTTAAGGTCGCAgagacaaccg                                                                                                >  1:269500/1‑242 (MQ=255)
                                                                                                                                              aTTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTGCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGGTTAGCGTCGCAGTGACGACGCCa                                                                                              <  1:378572/251‑1 (MQ=255)
                                                                                                                                                  aaGCGGTCCCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGcc                                                                                           <  2:283157/250‑1 (MQ=255)
                                                                                                                                                   aGCGGTCCCGGCGATCTTGCTGTGGTGCCGGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTGTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGa                                                                                         <  2:445126/251‑1 (MQ=255)
                                                                                                                                                       gttccGGCGATCTTTCTGGGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑GCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATTCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCaa                                                                                     <  1:418483/248‑1 (MQ=255)
                                                                                                                                                                    ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTTGCGTCGCCGTGACGCCGGCAGCCGGTCAAttcacattaa                                                                           >  2:165486/1‑239 (MQ=255)
                                                                                                                                                                    ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTaa                                                                           <  1:165486/248‑1 (MQ=255)
                                                                                                                                                                         gggggCCTGCCGGGCGTAAAAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAT‑TCAGGTTTTCAGGCG‑CAATGGAAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGGGACGACGCCAGCCGAGCAATACCAATTAAcggtccgg                                                                   <  2:367211/247‑8 (MQ=255)
                                                                                                                                                                                     gggggaaaaGCGAAATGACCCTAACCAGCAACTGGGCGAATCCCAGGTTTTTAGGTGAATTTTTGCTAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAACGTCAAAGCAGCTGGTTTGCTGATAACCTCGGTGAGCGTTTTGTTTCGGGCAATGATACGGGCTGGGAAAACGTCCCCGCGTTTAGCGGCGCAGTGACGAGGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTAc                                                       <  2:45351/243‑1 (MQ=255)
                                                                                                                                                                                                gAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCg                                                                                                                                          <  1:70724/157‑1 (MQ=255)
                                                                                                                                                                                                gAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCg                                                                                                                                          >  2:70724/1‑157 (MQ=255)
                                                                                                                                                                                                             ccAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGa                                 >  1:25410/1‑249 (MQ=255)
                                                                                                                                                                                                             ccAGCAACTGGCCGCAGCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGa                                 <  2:25410/249‑1 (MQ=255)
                                                                                                                                                                                                                         cGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGt                                                         >  1:355533/1‑213 (MQ=255)
                                                                                                                                                                                                                         cGCATCCCAGTTTTTTAGGTGAGTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTGCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATGCCAATGAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTc                   <  1:223283/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TAAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGCCTTTTATGGCGTTCt                  >  1:476684/1‑251 (MQ=255)
                                                                                                                                                                                                                              cccAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGt              >  1:495646/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aGGTGATTTTTTGCTAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAAGGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTATCGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACg    <  1:465693/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTGATTTTTTGCCAGCCACACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTGGCGTAGCCGTGGCGACGCCAGCCGATCCATACCAAATCACGGAACGGGCGAATTAATACCACACTTTGCTGATTTCACTGGCTTTTCTGGCGTTCTTTTTTTTTGAAAcacc  >  2:61603/1‑248 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
GTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCT  >  NC_000913/4638640‑4639124

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: ‑