Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,005 G→C P276P (CCG→CCC)  creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,0050GC93.3% 128.4 / 2.9 45P276P (CCG→CCC) creDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base C (20/22);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

AACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCG  >  NC_000913/4638757‑4639243
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
aaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTCCCGGCGATCTTTCTGTGGTGACTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGCGATATTTTGCCAGCCAAACGGGCAGTTAGCGAGTCAGGTTTTCCGGAGCAATGGCAAAGCAGCTGGTTTGCTAATACTCTCGGTTAACTTTTTGCTTAAGGCTATGATACCGGCTAGGGAAACTTTCCCgg                                                                                                                                                                                                                                               >  2:236609/1‑251 (MQ=255)
   cTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCcgcggt                                                                                                                                                                                                                                           >  1:631748/1‑249 (MQ=255)
                  gCCCTGAATTTAAGCGGTACCGGCCATCTTTCTGGGGTGTCTG‑CGGTGCGAATTGGCAAAATTACCTTAACCAGCAACTGGCCGCATCCCATTTTTTTAGGTGATTTTTTGCAAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGAAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTGCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGTAGtgac                                                                                                                                                                                                                            <  2:404862/251‑4 (MQ=255)
                              aGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGGCGCccgccg                                                                                                                                                                                                                 >  1:397812/1‑250 (MQ=255)
                                    cccGGCGATCTTTCTGTGGTGCCTCCAGG‑GTGTAATAGCTAAATGCCCTTAACCAGCAACTGGCCGTATCCCAGTTTTTTAGGTGATTTTTTGCCCGCCAAACGTGAAGTTAGCGAGTCAGTTTTTCAGGCTCAATGTCAAAGCCGCTGGTTTGTTAATAATCTCGGTGAGCTTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTTGCGTCGCAGTGACGACGCCAGCCGATCAAt                                                                                                                                                                                                           <  2:524681/249‑1 (MQ=255)
                                         cGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCCACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGGAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCCGCCCGCCCGTCAATACCact                                                                                                                                                                                                     >  2:351637/1‑249 (MQ=255)
                                               ttCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTa                                                                                                                                                                                                   <  1:20705/247‑1 (MQ=255)
                                                 gtccgatctctgccgg‑gcgTAATAGCGAAATGACCTTAACCAGCAACCGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAAcgga                                                                                                                                                                                              <  2:17474/237‑4 (MQ=255)
                                                   gtggtgCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGAc                                                                                                                                                                                             <  1:420908/249‑1 (MQ=255)
                                                          ctgccgg‑gcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAAcgga                                                                                                                                                                                              >  1:17474/5‑238 (MQ=255)
                                                            gccgg‑gcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTa                                                                                                                                                                                   <  1:342657/248‑1 (MQ=255)
                                                                       acgtgttctcttccggtctcAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGGCGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCAtt                                                                                                                                                                        <  1:443572/232‑1 (MQ=255)
                                                                         aGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAAt                                                                                                                                                                                                           >  1:503138/1‑214 (MQ=255)
                                                                             aaaTGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGAAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATAACGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGAGACGACGCCAGCCGATCAATACCCATTAACGGACCGGGCGAATAAGTACGCCATTctgctg                                                                                                                                                                  >  1:622077/1‑251 (MQ=255)
                                                                                      tagccAGAAACTGGCCTCATCCCATTTTTTTAGGTGATTTTCTACCATCCAAACGGGAAGGTATCTAGTCAGGGTGTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAAAAATCGCGGTGGGCGTTTGGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCGGCTTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACt                                                                                                                                                          <  2:649176/247‑1 (MQ=255)
                                                                                          cAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTCCCACCCAACCGGGAGGTTAGCGGGTCAGGTTTTCAGGCGCATTGCCAAAGCAGCTGGTTTGATAATAATCTCGGTGAGCGTTTTGCTTCAGGCAAGGATACCGGCT‑GGGAAAACTTCCCCGCCTTTAGCGCCCCAGGGCCAACGCCACCCGATCAATACCAAATAAGGGAACGGGCGAATAAGTACGACAtataacaggaagagcgccggg                                                                                                                                                     >  2:443572/1‑231 (MQ=255)
                                                                                           aGCATCTGGCCGCAGCCCAGTTTTTGAGGTGATTTGTTGCGATCCAAACGGGAAGTGAGCGTGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGGTCGTTTTGCTTCTGGCAATGATACCGGCT‑GGGAAAACGTCCCCGGGTTTAGCGGCGCAGTGTGGACGTCAGCCGATCAACACCAATTAACGGTCCGGGCGCCTAAGGACGCCATTCTGCTGATTGCACTGACttt                                                                                                                                                    <  1:415290/251‑1 (MQ=255)
                                                                                                          cccAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTg                                                                                                                                      <  1:477872/250‑1 (MQ=255)
                                                                                                             aGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTTCTGATTGCCCTGACTTTTTTGGCGTTCTTTGtttt                                                                                                                                  >  2:220162/1‑251 (MQ=255)
                                                                                                                         gATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCa                                                                                                                       <  1:219591/250‑1 (MQ=255)
                                                                                                                                  ccAGCCAAACGGGACGTTAGCGAGTCAGGTTTTCCGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGGCGCCGCCAGCCGATCAATCCCCATTCACGGCCAGGGCGCCTAAGAACGCCATTCTGCTGATTGCACTTAATTTTATGGCGTTTTTTGttttttacaagctcacagcgacaac                                                                                                             >  2:614152/1‑231 (MQ=255)
                                                                                                                                    aGCCAGACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                                                                                                           >  1:86822/1‑251 (MQ=255)
                                                                                                                                    aGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                                                                                                           >  1:203731/1‑251 (MQ=255)
                                                                                                                                     gCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGttt                                                                                                          <  1:351637/251‑1 (MQ=255)
                                                                                                                                            gggAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGCTCAATACCAATTAACGGACCCGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGACAAGCTCACCGCGCAACGTTTACAccc                                                                                                    >  2:471254/1‑250 (MQ=255)
                                                                                                                                                     gCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCCatat                                                                                          >  1:400170/1‑251 (MQ=255)
                                                                                                                                                            aGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCtgg                                                                                   <  1:220162/251‑1 (MQ=255)
                                                                                                                                                              gTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCCATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGCCtt                                                                                                                                                     >  2:108677/1‑183 (MQ=255)
                                                                                                                                                              gTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACtt                                                                                                                                                     <  1:108677/183‑1 (MQ=255)
                                                                                                                                                                         caatGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTGCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCAtt                                                                      <  1:212603/248‑1 (MQ=255)
                                                                                                                                                                          aatGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGATCAATACCAATTAAAGGACCGGGCGGCTAAGTACTCCATTCTGCTGCTTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACCACCACCGCGCAACGTTTTAAACCAATGCAATTTTTTATGGTGGGGCttttcttg                                                                     >  2:93008/3‑246 (MQ=255)
                                                                                                                                                                           atGGCAAAGAAGCTTGTTTGCTAATAATCGCGGTGCGCGTTTTGCTCTAGGCAATGATACCGGTT‑GGTAAAACTGCCCCGCTTTTAGCGTCGCAGTGACGCCGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGTTGGGGCTTTCATTgg                                                                    <  2:510358/250‑1 (MQ=255)
                                                                                                                                                                             ggCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGGACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTg                                                                  <  1:640611/251‑1 (MQ=255)
                                                                                                                                                                                 aaGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCAGCCAGCCGATCAATACCAATTAACGGAACGGGCGACTAAATACGCCATTCTGCTGATTGCACTGCCTTTTTTGGCGTTTTTTGTTTTTGAAACGCTCAACGCGCAACGTTTACAACCAATGCAACTTTTTCTGGGGGGGCTTTTCTTGGtgttgt                                                              >  2:474330/1‑251 (MQ=255)
                                                                                                                                                                                      gcTGGTTTGCTAATAATCTCGGTGAGCGTTTTGATTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACCCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCCTTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGGTTTTGAAACCCTCACCGCGCAACGTTTACACCCAATGCCATATTTGCTTGTGGGGCTTTCATTGGTGCTGTTTTAt                                                         >  2:41578/1‑251 (MQ=255)
                                                                                                                                                                                         ggTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTAtt                                                        <  1:99315/249‑1 (MQ=255)
                                                                                                                                                                                                taataaTCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCACCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCAATTCTGCTGAATGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCCTTGGTGATGTTTTATTAGATCTTGc                                               >  1:374295/1‑251 (MQ=255)
                                                                                                                                                                                                 aataatCTCGGTGATGGTTTTGCTTCAGGCAATGATACCGGCT‑GGTAAAACTTCCCCGCTTTTAGCGTCGCTGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATAGTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCt                                              <  1:327000/251‑1 (MQ=255)
                                                                                                                                                                                                         cGGTTAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAAc                                                                                                              >  1:115832/1‑179 (MQ=255)
                                                                                                                                                                                                         cGGTTAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCGGCGTTTACCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAAc                                                                                                              <  2:115832/179‑1 (MQ=255)
                                                                                                                                                                                                         cGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCttt                                      <  1:471254/251‑1 (MQ=255)
                                                                                                                                                                                                           gTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCt                                    >  1:76486/1‑251 (MQ=255)
                                                                                                                                                                                                                       cTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCCACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGt                        >  1:641437/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aTACCGTCT‑GGGAAACCTTCCCCGCTTTTGGCTGCGCAGTGACGACTCCAGCCGATCAATACCAATTAACGGACCGGGCTACTAAGTACGCCATTCTGCTGATTGCACTGACTGTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTGACACCCAATGCAATATTTGCTGGTGGGGCTTTCAGTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCa            <  2:76486/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGCT‑GGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGa        <  1:614152/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gggAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATc   <  1:430129/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ggAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCg  >  1:638004/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
AACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCG  >  NC_000913/4638757‑4639243

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