Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,883 C→T 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,8830CT92.9% 114.9 / ‑2.3 42L236M (CTA→TTG) creDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (15/24);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCA  >  NC_000913/4638654‑4639126
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ggaatttaccggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAACCGGGAAGTTa                                                                                                                                                                                                                                 <  2:14395/242‑1 (MQ=255)
       aacggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGGAATAGCGAAATGACCTTAACCAGCCACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTc                                                                                                                                                                                                                          >  1:234905/1‑251 (MQ=255)
         cggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGaatc                                                                                                                                                                                                                          >  1:614657/1‑246 (MQ=255)
          ggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTGTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAgg                                                                                                                                                                                                                       <  1:191940/251‑1 (MQ=255)
                        cATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCaat                                                                                                                                                                                                         <  1:577543/251‑3 (MQ=255)
                        cATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCaat                                                                                                                                                                                                         <  1:440251/251‑3 (MQ=255)
                                              gagCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCtaa                                                                                                                                                                                   <  1:132353/251‑1 (MQ=255)
                                                   aGGCGGGCAGGGAGGGCATATCCCTTTACCTGAAGGGGAATGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAACGGTACCGGCGATCTTTCTGGGGTGCCTGGCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCAGCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGTAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAAtc                                                                                                                                                                              <  2:283632/251‑1 (MQ=255)
                                                     gCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAAtctc                                                                                                                                                                            <  1:386817/251‑1 (MQ=255)
                                                             ggCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCg                                                                                                                                                                    <  1:112640/251‑1 (MQ=255)
                                                                  gCATATCCCTTTACCCGAAGGGGACTGGCGTAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTg                                                                                                                                                               <  1:658643/251‑1 (MQ=255)
                                                                            ttACCCGAAGGGGACTGGCTGAAGCATAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCTATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCa                                                                                                                                                      <  1:555033/250‑1 (MQ=255)
                                                                              aCCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGCTCTTTCTGTGGTGCCTGGCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATg                                                                                                                                                   >  1:649176/1‑251 (MQ=255)
                                                                                ccgAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGAt                                                                                                                                                 <  1:503504/249‑1 (MQ=255)
                                                                                    aGGGGACTTGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGGGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCg                                                                                                                                             <  2:234905/251‑1 (MQ=255)
                                                                                          cTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑ggg                                                                                                                                      >  1:528808/1‑251 (MQ=255)
                                                                                             gCGGAAGCAGAACCTTAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGaaa                                                                                                                                   <  1:352476/251‑1 (MQ=255)
                                                                                             gCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGaa                                                                                                                                    >  1:509598/1‑250 (MQ=255)
                                                                                                 aaGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑g                                                                                                                                        <  1:365952/242‑1 (MQ=255)
                                                                                                 aaGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGATTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑g                                                                                                                                        >  2:365952/1‑242 (MQ=255)
                                                                                                       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)
                                                                                                                                                                                         aCCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTaataat                                                                                                                                                                               <  1:296238/117‑1 (MQ=255)
                                                                                                                                                                                         aCCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTaataat                                                                                                                                                                               >  2:296238/1‑117 (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)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
TGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGG‑GCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCT‑GGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCA  >  NC_000913/4638654‑4639126

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