Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913689,5140AC34.4% 114.3 / 11.3 62G456G (GGT→GGG) lntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/6);  new base C (0/21);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 2.43e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATAATGCACTGATGCCGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCATTTGCC  >  NC_000913/689272‑689754
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
gATAATGCACTGATGCCGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCgg                                                                                                                                                                                                                                           >  1:462442/1‑251 (MQ=255)
                     ccGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGcc                                                                                                                                                                                                                     >  1:371654/1‑251 (MQ=255)
                           aGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CCGGGGTAGGC‑GTCACG‑TTAGTGGATAACACCTC‑GCGGGTGAACTGCGGGATCATCGACTGAATCTCACCCTACGGGCCAATCA‑CCGCCGTAATGCCgttgtg                                                                                                                                                                                                               >  2:115108/1‑250 (MQ=255)
                              gtgcgaaaaTATGGGATGTTTACGGGCACGATAAGAAGGGATTTTTTCCGTCGCGAACGAGGCCTCATCAACACAGCAGCAAACCCAAACAATGCCGTCAGCCCCCCCGGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGCCGGC‑GTCACG‑TTAGTGGTTACCACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGGGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGt                                                                                                                                                                                                            <  2:75137/245‑1 (MQ=255)
                              gTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGt                                                                                                                                                                                                            >  1:524383/1‑251 (MQ=255)
                                        aggggATGTATCCCGGCCCGAAAAGAAGGTATTATTTCCGTCGCTTCCGCGGACCCATCACCACGGCAGCAAAACCAACCAACGCCGTCGGCCCCCACAGCGGCCAGGT‑CCCGGGCCGTGCGTTTGGTTT‑GATTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTACCACCTC‑GCGGTTGAACTGCGGGATCATCGCCTGATTCTCCCCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCaa                                                                                                                                                                                                   <  1:1346/248‑1 (MQ=255)
                                               gtattcgcgCACAATAAGAGGTGATTGTTTACGTCGCTGCCGCAGACTAATCAACACAGAAGCAAAACAAACCAATGCCGCCAGCCCCCACAGCGGCCATTT‑GCCGGTACGTGCGTAGGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACT‑TTAGTGGTTACCCCCTC‑GGGGGTAAACGGGGGGATCTCCGCCTGAATCCCCCCCTGCGGGCCAATCC‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgc                                                                                                                                                                                           <  2:550350/243‑1 (MQ=255)
                                                 aTTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACTCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgc                                                                                                                                                                                         >  1:285983/1‑251 (MQ=255)
                                                       gcccGATAAGAAGGGATTATTTCCGTCGCTGACGCAGACTCACCAACACACCAGCAAAACCAAACAATGCCGGCAGCCCCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACCCCTC‑GCGGGTGAACTGCGGGACCATCGCCTGAATCTCCCCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGc                                                                                                                                                                                   <  1:53934/248‑1 (MQ=255)
                                                          gaaaaaaaagGGATTTTTCACGTCGCGCCCGCACACTCATAAACACAGCAGCAAAACAAAACAATGCCGACAGCCCCCAAGGAGGCGGTTT‑GCGGGTACGTGCGTGTGGTGT‑GATTC‑GGGTGGTCGGC‑GCCACG‑TTAGTGGTTAACCCCGC‑GGGGGGTACCGGTGGGATCCTCGCCGAAATCTCCCCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTcc                                                                                                                                                                                <  2:333892/244‑1 (MQ=255)
                                                             aaagaagGTATTATTAACGTGGCTGACGCAGACTCACCAACACAGCAGCAAACCCAAACAAGGCCGTCGGCCCCCACAGCGGCCAGTT‑GCCGGTACGGGCTTATGGTGT‑GAGGC‑CGGTGGTCGGC‑GTCACG‑TTTTGGGTTAACACCTC‑GCGGGTGAACTGCGGGAACACCGCCTAAATCTCCCCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGGGGGCGCGCCAGCTCCAgc                                                                                                                                                                             <  1:437736/250‑1 (MQ=255)
                                                                                   gtGACGCACACTCAACAACACAGCCGCAAAACCAAACAAGGCCGTCAGCCCCCACGGCGGCCAGTT‑GCCGTTATGTGGGTCTGGTGT‑TAGTC‑CGTTGGCCGGC‑GTCACT‑TTAGTGGTTAACCCCTC‑GGGGGTGAACTGGGGGATCATCACCTGATTCTC‑CCCTTTGGGGCAATAAGCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTTGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                                                                                                                       <  1:255377/250‑1 (MQ=255)
                                                                                   cggACGCAGCCTCATCAACACAGAAGAAAACCCAAAAAAGGCCGTCAGCACCCACAGGGGCCAGTT‑GCCGTTCTGTTCGTATGGTGT‑GAGTC‑CGGGGGTCGGC‑GTCACG‑TTAGTGTTTAACCCCGC‑GCGGGTGAACTGCGGGATCTTCGCCTGAATCTCCCCCTGCGGACCAATCA‑CCGCCGTAATGCCTGTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                                                                                                                       <  1:528797/249‑1 (MQ=255)
                                                                                     aaCGCAGACCCACCAACACAGCAGCAAAACCAAACAAGGCGGTCATACCCCACAGCGGCCGGTG‑GCCGGTAGGTGCGTAGGGTGT‑GAGTC‑CGGGGGTCGGC‑GTCACG‑TTAGGGGTTAACCCCTC‑GCGGGGGACCGGTGGGATCATCGCCTGAATTTCCCCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAg                                                                                                                                                     <  1:105409/250‑1 (MQ=255)
                                                                                                  caaaaaacaacaaaaaCCAAACAGTGCCGTCAGCCCCCCCAGGGCCCAGTT‑GCGGGTAGTTGCGTAGGGTGT‑GAGCC‑CGGGTGCGGGT‑GTCACG‑TTAGGGGTTAACCCCCC‑GCGGGTGAACGGCGGGATCATCGCCGGACCCCCCCCCGGCGGCCCAACCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAc                                                                                                                                        <  2:156545/239‑1 (MQ=255)
                                                                                                  caACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑AGGTCTCGG‑GGGCCGC‑GGCATC‑ATAGAGAATAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTTCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCCTGGAc                                                                                                                                        >  2:422213/1‑251 (MQ=255)
                                                                                                     cacaGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCa                                                                                                                                      >  1:91678/1‑250 (MQ=255)
                                                                                                        agcagcAAAGCCAAACAGGGCCGTCACCCCCCACGGCGGC‑AGTTTGCGGGTCGGTGCGTGTGGTGT‑GGGTC‑CGGGGGTCGGC‑GTC‑CC‑TTGGGGGGTAACCCCCCCGGGGGGGACGTGCGGGATAAGCCCCTGAACCCCCCCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAg                                                                                                                                  <  2:397933/251‑1 (MQ=255)
                                                                                                             caAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGTAGATCC‑CGGTGGCCGCC‑G‑CATC‑ATAGAGAATAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGACATCTGGAAGTGTTGCCATGGACCAATAGAttt                                                                                                                              >  2:485226/1‑250 (MQ=255)
                                                                                                             caAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTC‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGGCGGC‑GGCACG‑TTAGTGGTAAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGGTGGGGCTGCGCAACAGTGGGCGCGACAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGATGCCATGGACCAATAGAttt                                                                                                                              >  2:114884/1‑250 (MQ=255)
                                                                                                                    aaaCAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑AAATCTCGGTG‑TCGGC‑GTTATG‑ATAGTGGTTAACACCTC‑GAGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGGGGGCGCGCCCGCTCCAGCGCACGAATTCGCGCCATCTGGAAGTGTTGCCATGGGCCAATAGAGTTACAAAAcc                                                                                                                      >  2:284430/1‑251 (MQ=255)
                                                                                                                         aTGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GGATCTCGG‑GGGCGGCCGT‑ATGATTA‑AAAAAAACACATC‑GCGGGTGAACTGCGGGATCATCGCATGAATATAACCCTGCGGGCCACTCA‑CCGCCGGAATGCCGTTGTTGGGGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCAGTCGCGCCATCTGGAAGGGTTGCACTGGACCAATAGATTTACCAAACCACGcc                                                                                                                 >  2:437599/1‑250 (MQ=255)
                                                                                                                              gTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGtt                                                                                                            >  1:464804/1‑251 (MQ=255)
                                                                                                                                 aGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑AAGTC‑CGGGGGGCCGC‑GTCACG‑ATAGTAGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGAGTTGCCCTGGACCAATAGATTTCCCAAACCACGCATCGATgg                                                                                                          >  2:261487/1‑250 (MQ=255)
                                                                                                                                    aCCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCa                                                                                                                                                                                     >  1:170032/1‑172 (MQ=255)
                                                                                                                                    aCCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg                                                                                                          >  1:343538/1‑247 (MQ=255)
                                                                                                                                      ccacaagcggccgttggccgttaggTGCGTATGGTGT‑GGGTC‑CGGTGGTGGGC‑GTCACG‑TTAGTGGTTAACCCCTC‑GCGGGTGACCTGGGGGATCACCGCCTGAATCCCCCCCGGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATa                                                                                                     <  2:546914/227‑1 (MQ=255)
                                                                                                                                         caGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGGCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTc                                                                                                  >  2:36484/1‑250 (MQ=255)
                                                                                                                                                    t‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCg                                                                                      >  1:82717/1‑251 (MQ=255)
                                                                                                                                                                gggTATGGTGG‑GAGGC‑GGGTGGTGGGC‑GTCACG‑TTAGGGGTTAACCCCCC‑GCGGGTGAGCTGCGGGTTCGTCCCCGTACTCTCCCCCGGCGGCCCAATCC‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGc                                                                            <  2:372924/248‑1 (MQ=255)
                                                                                                                                                                   gtatgtgttgtgTC‑CGGGGGTCGGC‑TTCCCG‑TTGGTGGTAACCACCTC‑GCGGGTGAACTGGGGGATACTCGCCAGAATCTCACCCTGCGGGCAAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGa                                                                         <  1:261487/240‑1 (MQ=255)
                                                                                                                                                                     tgtggt‑gAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACCCCTC‑GCGGGGGAACTGGGGGATCATCGCCTGAATCTCACCCGGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAg                                                                       <  1:263029/246‑1 (MQ=255)
                                                                                                                                                                        tgt‑gAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTAt                                                                   >  1:370170/1‑251 (MQ=255)
                                                                                                                                                                            gagcc‑cGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACCCCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAACCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATcgc                                                                <  1:66853/247‑1 (MQ=255)
                                                                                                                                                                             attC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGtgct                                                                                                                                                                                                         >  2:279279/3‑113 (MQ=255)
                                                                                                                                                                              ggtcgggtgtgcGGC‑GTCAGG‑TTGGTGGTAAACCCCTC‑GGGTGTGAGCTGGGGGATCCTCGCCTGAACCGCCCCCTGCGGGCAAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATcgcg                                                               <  2:343538/240‑1 (MQ=255)
                                                                                                                                                                                   ggtggtCGGC‑GTCAGG‑TTAGTTTTTACCCCCTC‑GCGGGTCAACTGCGGAACCATCGCCTCACTCCCCCCCTGCGGCCCACTCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACtt                                                          <  2:524383/251‑1 (MQ=255)
                                                                                                                                                                                   ggaggtCGGC‑GTCACG‑TTCGTGGTTAACACCTC‑GCGGGTGAACTGCGAGATCATCGCCTGAATCTCACCCTGCGGGCAAATCA‑CCCCCGTACTGCCGTTGTTGGTGCTGCGCAACAGCGGGCGCGCCGGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGTACCAATAAATTTAACAAACCACGCATCGTTGGAGATAGTCAGCCGATAGTCGGTATCCGGTCGGCAGTTCTCTCGCACtt                                                          >  2:27382/4‑251 (MQ=255)
                                                                                                                                                                                   agtggtCGTC‑GTCACC‑TTAGTGGGTAACACCTC‑GCGGGTGAAATGCGGGATCATCGCCTGAATCTCAACCTGCGGGCCAATCA‑CCAACGTAATGCCGGTGGTGGTGCTGCTCAACAGGGGGCGCGACGGCTCCAGCGCACGCAGTCGAGCCATCTGGAAGTGTTGCAAGGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATGGTCGGTATCCGAGCGGAAGTTATCGTGCACtt                                                          >  2:400633/2‑251 (MQ=255)
                                                                                                                                                                                      ggtCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCt                                                       >  1:374762/1‑251 (MQ=255)
                                                                                                                                                                                        ggtcggcgccaggtgagggtttAACCCCGC‑GCGGGTGAACTGGGGGATCATCGCCTGAATCCCCCCCTGCGGGCCAACCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGc                                                        <  1:134558/230‑1 (MQ=255)
                                                                                                                                                                                          ggC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGcc                                                   >  1:509221/1‑251 (MQ=255)
                                                                                                                                                                                            gcgccgcG‑TTGGTGGTTACCCCCTG‑GCGGTTGAACTGCGGGCTCAGCGCCTGAACCCCCCCCTGCGGGCAAATCG‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCg                                                  <  1:115108/245‑1 (MQ=255)
                                                                                                                                                                                                  cG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATg                                            >  1:504581/1‑251 (MQ=255)
                                                                                                                                                                                                      tAGTGTTTAACCCCGC‑GCGGGTGAACTGCGGGTTCATCGCCTGAATCTCCCCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAt                                          <  1:117906/250‑1 (MQ=255)
                                                                                                                                                                                                       ggTGGTTACCCCCTC‑GGGGGTGAACTGCGGTATCGCCCCCTAAACCTCACCCTGGGGGCCAACCC‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAtct                                        <  2:298275/250‑1 (MQ=255)
                                                                                                                                                                                                        gTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGCTTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAtctc                                       >  2:92734/1‑251 (MQ=255)
                                                                                                                                                                                                        gTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAtctc                                       >  1:387432/1‑251 (MQ=255)
                                                                                                                                                                                                         gggTTAACCCCTC‑GCGGGTGAACTGGGGGATCAGCGCCTGAATCCCCCCCTGCGGGCAAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCg                                      <  1:538793/250‑1 (MQ=255)
                                                                                                                                                                                                             tAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGc                                  >  1:341357/1‑251 (MQ=255)
                                                                                                                                                                                                                             gAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTa                    >  1:442372/1‑250 (MQ=255)
                                                                                                                                                                                                                              atcTGCGGGATCATCGTCTGAATCTCCCCCTGCGGCCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTaa                   <  2:134636/248‑1 (MQ=255)
                                                                                                                                                                                                                               aCTGCGGGTTCATCGCCTTATTTCCCCCCTGCGGGCCAATCG‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt                                                                                                      <  1:281803/166‑1 (MQ=255)
                                                                                                                                                                                                                               aCTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt                                                                                                      >  2:281803/1‑166 (MQ=255)
                                                                                                                                                                                                                                  gctggttcatccctggAATCTCCCCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCt                <  1:18366/234‑1 (MQ=255)
                                                                                                                                                                                                                                    gggATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAAt            >  2:218946/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gggATCACCCCCTAAATCTCCCCCGGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAAt            <  2:341357/251‑1 (MQ=255)
                                                                                                                                                                                                                                       atccccgcgtgaaccccACCCGGCGGGCAAACCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGAt                                                                                           <  2:100270/153‑1 (MQ=255)
                                                                                                                                                                                                                                       atcatcGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATACCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGAt                                                                                           >  1:100270/1‑169 (MQ=255)
                                                                                                                                                                                                                                          atcGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCAtt      >  2:292076/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCATTTGc   >  2:274708/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ccTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGCAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGTGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCATTTGcc  >  2:280989/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
GATAATGCACTGATGCCGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTT‑GCCGGTACGTGCGTATGGTGT‑GAGTC‑CGGTGGTCGGC‑GTCACG‑TTAGTGGTTAACACCTC‑GCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCA‑CCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCATTTGCC  >  NC_000913/689272‑689754

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