Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,379,6380GA28.3% 85.1 / 13.9 46F45F (TTC→TTT) frdCfumarate reductase membrane protein FrdC
Reads supporting (aligned to +/- strand):  ref base G (29/4);  new base A (0/13);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 2.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCAACAAACAGGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTT  >  NC_000913/4379389‑4379880
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
ggCAACAAACAGGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTATGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCg                                                                                                                                                                                                                                                      >  1:575420/1‑250 (MQ=255)
      aaacaGGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATc                                                                                                                                                                                                                                                >  2:137934/1‑250 (MQ=255)
               tACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGc                                                                                                                                                                                                                                    >  1:86714/1‑251 (MQ=255)
                                                             ggCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTaaa                                                                                                                                                                                       >  2:563853/1‑250 (MQ=255)
                                                                   ggTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGc                                                                                                                                                                                    >  1:392051/1‑247 (MQ=255)
                                                                     tCCCATTTTTTCGTCTTTTACGATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGGTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATTATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATaa                                                                                                                                                                               >  1:78328/1‑250 (MQ=255)
                                                                        cATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACg                                                                                                                                                                           >  1:445882/1‑251 (MQ=255)
                                                                                  tCTTTTACATGGATTTTGGCCGTTTCCGGTGCCAGTTAAACCCGGGTTTGGGTGGGAAGCAGGGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATACCCGGTTTTT‑TAAAAAAGTCGACGAACCCCGCCCAGCCTCCGGGGCCTTTTTT‑C‑AGGGCAAAC‑AGCCCAAAAATCA‑GTTAAA‑TGCTAACCCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAAttt                                                                                                                                                                    <  1:137934/248‑1 (MQ=255)
                                                                                          aaTGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAgg                                                                                                                                                         >  1:581351/1‑251 (MQ=255)
                                                                                              atatTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTgg                                                                                                                                                      >  2:341746/1‑250 (MQ=255)
                                                                                              atatTGGCCGCTTTCGGTGCCAGTTCAAACCAGGCTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGCAAAAGCGATAAAACGGCAATTTTTTccacccggggg                                                                                                                                                      >  2:576188/1‑244 (MQ=255)
                                                                                                          ttCGGTGCCAGTTCAACCCAGGTTTGGGTGGGCGGAAGAGCGGCCGCCAGAGGTATCAGGTTAAGGAACAAGATAACCGGTTTTT‑TAAAAAAGTCGAAAAACCCCCCCCAGGCT‑CCGGGCCATTTTTTC‑AGGGCAAAC‑CGCCAAAAAACCA‑TTTCAA‑TGCGGACCCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTgg                                                                                                                                            <  1:99174/248‑1 (MQ=255)
                                                                                                            cGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGt                                                                                                                                        >  1:595920/1‑250 (MQ=255)
                                                                                                             ggTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTAc                                                                                                                                      >  2:460688/1‑251 (MQ=255)
                                                                                                                     ttaaacccggtttttGGTGGGCAGCAGGGCGGCCGCAGGAGTGATCAGTTTAATGATCACAATAACCGGGTTTT‑GAAAAAAGTCGACGAATCCCCCCCAGGCTCCGGGGCCTTTTTT‑C‑AGGGAAAAC‑AGCCAGAAAATCA‑GTTAAA‑TGCTGACCCACACGGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGt                                                                                                                               <  1:460688/239‑1 (MQ=255)
                                                                                                                               gTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGGTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgt                                                                                                                    >  2:399789/1‑251 (MQ=255)
                                                                                                                                       gtgCAGCAGGGCGGCCGCAGGAGTGACCAGTTTAATGACCGCAATACCGGGGTTTT‑GAAAAAAGTCGACAAACCCCCCCCAGGCTTCGGGCCCTTTTTT‑C‑AGG‑CAAACCAGCCCAAAAATCA‑GTAAAA‑TGCGAAACCACCCACCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGt                                                                                                            <  2:591352/251‑1 (MQ=255)
                                                                                                                                            aaggggcggacgccAGAGTGATAGGGTTAATAATCAAGATA‑CCGGGTTTTTGAAAAAAGCCGACGACCCCCGCCCAGGCTTCCGGGC‑ATTTTT‑CCAGGACAAAC‑AGCCCAAAAATCA‑GTTAAA‑TTCTGACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCAc                                                                                                       <  2:293265/241‑1 (MQ=255)
                                                                                                                                              agagCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTc                                                                                                     >  2:546470/1‑251 (MQ=255)
                                                                                                                                                       gccAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTg                                                                                            >  1:381529/1‑251 (MQ=255)
                                                                                                                                                              tGATCGGGTTAATGATCCCGATA‑CCGGGTTTTTGAAAAAAGTCGACGAACCCCCCCCAGGCTTCCGGGCTTTTTTC‑C‑AGGACAAAC‑GCCCCAAAAATCA‑GTTCAA‑TGCTGACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATa                                                                                                                                                                                <  1:208858/160‑1 (MQ=255)
                                                                                                                                                              tGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATa                                                                                                                                                                                >  2:208858/1‑160 (MQ=255)
                                                                                                                                                               gATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAgg                                                                                    >  1:115123/1‑251 (MQ=255)
                                                                                                                                                               gATCAGGATAATGATCACGATCACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAgg                                                                                    >  1:464423/1‑251 (MQ=255)
                                                                                                                                                                   aGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTcg                                                                                >  1:523618/1‑251 (MQ=255)
                                                                                                                                                                    gtttaagaacacagATACCCGGGTTTT‑AAAAAAAGGCGACGAACCCCCCACAGTCTTCCGGGCTATTTTC‑C‑AGGGCAACC‑ACCCCGAAACTCA‑GTTAAA‑TGCTGACCCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTcgc                                                                               <  2:136900/238‑1 (MQ=255)
                                                                                                                                                                         aTGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATa                                                                           >  2:172058/1‑250 (MQ=255)
                                                                                                                                                                              cACGATACCCGGGTTTT‑GAAAAAAGTCGACGAACCCCGCCCAGGCTTCCGGGCCTTTTTT‑C‑AGGGAAAAC‑AGCCCAAAAATCA‑GTTCAA‑TGCGGACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGttt                                                                                        <  1:34442/232‑1 (MQ=255)
                                                                                                                                                                              cACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGttt                                                                                        >  2:34442/1‑232 (MQ=255)
                                                                                                                                                                                cGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTc                                                                   >  1:234987/1‑251 (MQ=255)
                                                                                                                                                                                agATACCCGGTTTTT‑GAAAAAAGTCGACGACCCCCCCCCAGGCTCCCGGGCCTTTTTC‑C‑AGGCCACAC‑ACCCAGAAAATCA‑GTTCAT‑TCATGACCCACCCAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTc                                                                   <  2:395478/250‑1 (MQ=255)
                                                                                                                                                                                  aTAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCtt                                                                 >  2:592140/1‑251 (MQ=255)
                                                                                                                                                                                     aCCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCg                                                              >  2:224325/1‑251 (MQ=255)
                                                                                                                                                                                        ggttttt‑ATAAAAAGTCGACAACCCCCGCCCAGGCTCCCGGGCCTTTTTA‑C‑GGGACAAAC‑ACCCCAAAAATCA‑GTCAAA‑TGCTAACCCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAAc                                                           <  1:245831/248‑1 (MQ=255)
                                                                                                                                                                                                   aaaatgcgcaaacccccgccc‑gGGTTCCGGGGCTTTTTT‑C‑AGGGACAAC‑AACCCAAAAAACAAGTTAAA‑TGCTAAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACttt                                                 <  2:445882/237‑1 (MQ=255)
                                                                                                                                                                                                   aaaaGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACttt                                                 >  2:504124/1‑251 (MQ=255)
                                                                                                                                                                                                        ccgaggacccccgcccgGGCTTCCGGGCCATTTTG‑C‑AGGGCAAAC‑AGCCCAAAAATAA‑GTTCAA‑TGCGGACCCCCACGGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCct                                            <  1:231016/242‑1 (MQ=255)
                                                                                                                                                                                                        ccaacgaACCCCGCCAGGGCTGCCGGGCCATTTTT‑C‑GGGGCAAAA‑CCCCCAAAAATAA‑GTTAAA‑TGCTAACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCct                                            <  2:500802/248‑1 (MQ=255)
                                                                                                                                                                                                           acgaATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GCTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct                                         >  2:134145/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccgggCTTCGGGGCCATTTTT‑C‑AGGGAAAAC‑AGCCCAAAAATCA‑GTTCAA‑TGCTGACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTc                                                                                 <  1:365829/197‑1 (MQ=255)
                                                                                                                                                                                                                      ccAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTc                                                                                 >  2:365829/1‑200 (MQ=255)
                                                                                                                                                                                                                             ccggggcttttttt‑C‑GGGGCAAAC‑AGCCCAAAAATAA‑TTTCAA‑GGCAGAACCCCCCGGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGAt                       <  1:305871/242‑1 (MQ=255)
                                                                                                                                                                                                                               ggggcttttttt‑C‑AGGGCAAAC‑AGCCCAAAAATCA‑GTTAAA‑TGCTGAACCACACAGCGGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCg                     <  1:341746/244‑1 (MQ=255)
                                                                                                                                                                                                                               cGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCg                     >  1:438213/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACg                  >  1:591772/1‑251 (MQ=255)
                                                                                                                                                                                                                                          t‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCgg          >  1:482236/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    acc‑AGCCCGAAAATCA‑‑TTCAAATGCTGACCCACCCAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACtt  <  1:290411/249‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
GGCAACAAACAGGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTT‑C‑AGGGCAAAC‑AGCCCGAAAATCA‑GTTCAA‑TGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTT  >  NC_000913/4379389‑4379880

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