Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,332,8230AG31.7% 74.4 / 12.9 42pseudogene (1580/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (2/26);  new base G (11/2);  total (14/28)
Fisher's exact test for biased strand distribution p-value = 1.69e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCATCC‑GGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCA‑GTTTTGAACTCTTTTTTGGCGAG  >  NC_000913/2332604‑2333043
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      
aaCCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGCGACGGTTAAAAAATAGCGGCTCGGCTTATC‑GGC‑GGCC‑GGGGGGAtttcat                                                                                                                                                                                                                                >  2:1746470/1‑226 (MQ=255)
            gTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑c                                                                                                                                                                                                                  <  1:1337797/229‑1 (MQ=255)
            gTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATC‑GGC‑GG‑C‑GGGGGGGTTTAAAGCcacctaaacgc                                                                                                                                                                                                                   >  1:1355782/1‑221 (MQ=255)
                          ttttGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑CGGGGGGGATTTACGCCCACATAACCC‑CGGG‑GGCCGGaac                                                                                                                                                                                                     >  1:560857/1‑226 (MQ=255)
                               aaGATTCCAGCGGGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTCCTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGttt                                                                                                                                                                                               <  1:1553681/228‑1 (MQ=255)
                                                        ccaacccGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGa‑ca                                                                                                                                                                            <  1:1835337/216‑1 (MQ=255)
                                                        cacaACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTTTCGGGC‑GG‑G‑GGGGGGATTTAACCCCCCATAACAC‑CTGG‑GGCGGGACCCCGTTTGCCTGCCTTCCAGCGC‑AAC‑GGGGa                                                                                                                                                                     >  2:1832215/1‑228 (MQ=255)
                                                           aaCCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGCGGG‑G‑GGGGGGGTTTAAC‑CCACACGACAC‑CGTG‑CGGCCGGCACCACGTTGCTGGCTTCCATCGA‑CCC‑gggggtaagt                                                                                                                                                                >  2:735400/1‑223 (MQ=255)
                                                           aaCCGACTCGCCGCTATGACTGTATTGTGCGGCAGTACTTAATGAGGAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGaa                                                                                                                                                                    <  1:736690/225‑1 (MQ=255)
                                                              cGACTCGCCGCTATTACTGTATTGTGCGTCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACAGGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CCC‑GGGGAAACGTc                                                                                                                                                               <  2:1633118/227‑1 (MQ=255)
                                                                                 tATTGTGCGGCAGTACTTAATGGGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑c                                                                                                                                         <  1:1082438/229‑1 (MQ=255)
                                                                                    tgtgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGATATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGa                                                                                                                                      <  1:1831451/229‑1 (MQ=255)
                                                                                       gCGGCAGTACTTAAATAGTAATGCTCTAGACCGCGTTCAATGATGATCTCTTTGGTGGTGGTGCCGCGATGCGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCGCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCa                                                                                                                                   <  2:884591/229‑1 (MQ=255)
                                                                                           cAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGaa                                                                                                                                                                    <  1:1306761/193‑1 (MQ=255)
                                                                                           cAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCCCGGGTTATC‑GGC‑GG‑C‑GGGGGGGTTTTAAGCCACACGAACG‑CGGG‑CGGGGGAGACCGCGTTGCTGCCTCCCACCGA‑GAC‑cgggga                                                                                                                                                                    >  2:1304506/1‑191 (MQ=255)
                                                                                                   aaTGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑Cg                                                                                                                                                                                                            <  1:1295406/147‑1 (MQ=255)
                                                                                                        gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTa                                                                                                                                           <  1:421297/205‑1 (MQ=255)
                                                                                                           aTGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAACCGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAgcgcgc                                                                                                              <  1:681291/229‑1 (MQ=255)
                                                                                                               gcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAc                                                                                                          <  1:818997/229‑1 (MQ=255)
                                                                                                                        cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGGGGTTTTAACCCCCATTACCGCCCGG‑GGCGGGAACCCGTTTGCGGGCCTTCCAGGGA‑CCA‑GGGGAAAGGCCCGGCATACACCAAAACGTTA‑CAGGCAAATGAAAA‑TTGTGGAGCGCGCAAAcatcacctt                                                                                                  >  2:1056572/1‑221 (MQ=255)
                                                                                                                        cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAgcgc                                                                                                                <  1:1058436/214‑1 (MQ=255)
                                                                                                                                   gATCGCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCggg                                                                                          <  1:1455112/224‑1 (MQ=255)
                                                                                                                                         tttGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATC‑GGCGGG‑G‑GGGGGGATTTACAGCC‑CAATACCG‑CTGGCCGCGGGAAACGGTTTGCCGGCTTTCCGGCG‑‑CACGGGGGAAACG‑CCGGGACAAAGCAAAACGCTATCCAACAATGGATAA‑TAGCTGGGGGC‑CAAACGTAAAACACACGCGGGCGTTTTTTa                                                                                 >  2:1441955/1‑228 (MQ=255)
                                                                                                                                                   gtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGGGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCAt                                                                    <  1:472287/229‑1 (MQ=255)
                                                                                                                                                   gtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGGGAATTAAAGCCACCTGACCC‑CCCG‑CGGCCGGCACCCGTTCGCTTCCTTCCCGCGACCCC‑GGGGAAACG‑CCCGCACAAACCCAATCGGTA‑CCGACCCTTGCGAA‑TTGCTGTGCGCCCCAAC‑CTAACCTCACCCGGGCCTTTTTTTCCGGGTCG‑ccc                                                                     >  2:471458/1‑227 (MQ=255)
                                                                                                                                                   gtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAGCGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCAt                                                                    <  1:675513/229‑1 (MQ=255)
                                                                                                                                                   gtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCAt                                                                    <  2:231386/229‑1 (MQ=255)
                                                                                                                                                               gccgcgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCATCC‑GGGTAGAgc                                                        <  1:1139279/228‑1 (MQ=255)
                                                                                                                                                                                  agaCGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTTACGCCACATGACCG‑CTGG‑CGTCGGaca                                                                                                                                                                                                     >  1:586388/1‑75 (MQ=255)
                                                                                                                                                                                  agaCGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGaca                                                                                                                                                                                                     <  2:585364/75‑1 (MQ=255)
                                                                                                                                                                                           aCAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCATCC‑GGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGc                            <  1:387897/228‑1 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGGGGTTTAACGCCCCCTGACCG‑CCGG‑CGTCGGGCACCGTTTTGCTGCCTTTCAGCGA‑CAC‑GGGGGAACGTCCGGCCTAACGCCAATCGTTT‑CCCCCCCTTGATAA‑GTGCTGTGCGCGCAAAC‑TTAACTGCACCCGGGCATTTTTTTCCGGCTCG‑CCCTCC‑CGGGCGAGCGGTTTCCGTTTGCCGCTTAACCCCTCGcccttgttt                    >  2:271115/1‑222 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGGGATTTAAAGCCACATTACCG‑c                                                                                                                                                                                                                  >  1:531797/1‑47 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑c                                                                                                                                                                                                                  <  2:530854/47‑1 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCC‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑cc                                                                                                                                        <  1:962276/117‑1 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGAACG‑CTGG‑GGGCGGGCACCGTTTCGCTGCCCTTCCGCGA‑CAC‑GGGGAAAAGTCCGGGATAACGGAAATCGTTA‑CCGCCCAATGATAATTTGCT‑AGCGCGCAAAC‑TTTAACGCCCCCGGGCCTTTTTTTACGGCGCGCCCCTCCGGGGTAG‑GAGGTGTCAGTTGTCCGGTTACCGCTTCCCCA‑Gtttt                    >  2:699247/1‑229 (MQ=255)
                                                                                                                                                                                                   gggCTCGGTTTATC‑GGC‑GG‑C‑GGGGGAATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCTCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCATCC‑GGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCa                          <  2:967804/221‑1 (MQ=255)
                                                                                                                                                                                                      cTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGTCCG‑CTGG‑CGTCGGACACCGTCTCGCTGCCTTCCATCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCAtt                                                                                                                                 <  1:1187920/120‑1 (MQ=255)
                                                                                                                                                                                                              aTC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCATCC‑GGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCA‑GTTTTGaa                 <  1:686231/219‑1 (MQ=255)
                                                                                                                                                                                                                  ggc‑gg‑c‑ggGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCATCC‑GGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCA‑GTTTTGAACTCTTTTTTGGc     <  1:1749447/228‑1 (MQ=255)
                                                                                                                                                                                                                   gc‑gg‑c‑ggGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCATCC‑GGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCA‑GTTTTGAACTCTTTTTTGGCg    <  2:559873/228‑1 (MQ=255)
                                                                                                                                                                                                                       g‑c‑ggGGTGATTTAACGCCACATGACCG‑CTTG‑CGTCGGGCACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAAAGCAAATCGTTA‑CCGACCATTGATTA‑TTGCTGAGCGCCCAAAA‑CTAACTGCACGCGGGCATTTTTTACCGGGCCG‑CCCACC‑GGGGCGGGCAGGTGCAGTTCTCCGCCGACCCCCTCCCCACGTTTTGAACCGTTTTTTgtgcgg  >  2:1390557/1‑223 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      
AACCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATC‑GGC‑GG‑C‑GGGGAGATTTAACGCCACATGACCG‑CTGG‑CGTCGGACACCGTTTCGCTGCCTTCCAGCGA‑CAC‑GGGGAAACGTCCGGCATAACGCAAATCGTTA‑CCGACCATTGATAA‑TTGCTGAGCGCGCAAAC‑TTAACTGCACGCGGGCATTTTTTACCGGCTCG‑CCATCC‑GGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCA‑GTTTTGAACTCTTTTTTGGCGAG  >  NC_000913/2332604‑2333043

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