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→GGGlntapolipoprotein 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.

GATAATGCACTGATGCCGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCATTTGCC  >  NC_000913/689272‑689754
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
gATAATGCACTGATGCCGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCgg                                                                                                                                                                                                                                           >  1:462442/1‑251 (MQ=255)
                     ccGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGcc                                                                                                                                                                                                                     >  1:371654/1‑251 (MQ=255)
                           aGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCCGGGGTAGGCGTCACGTTAGTGGATAACACCTCGCGGGTGAACTGCGGGATCATCGACTGAATCTCACCCTACGGGCCAATCACCGCCGTAATGCCgttgtg                                                                                                                                                                                                               >  2:115108/1‑250 (MQ=255)
                              gtgcgaaaaTATGGGATGTTTACGGGCACGATAAGAAGGGATTTTTTCCGTCGCGAACGAGGCCTCATCAACACAGCAGCAAACCCAAACAATGCCGTCAGCCCCCCCGGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGCCGGCGTCACGTTAGTGGTTACCACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGt                                                                                                                                                                                                            <  2:75137/245‑1 (MQ=255)
                              gTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGt                                                                                                                                                                                                            >  1:524383/1‑251 (MQ=255)
                                        aggggATGTATCCCGGCCCGAAAAGAAGGTATTATTTCCGTCGCTTCCGCGGACCCATCACCACGGCAGCAAAACCAACCAACGCCGTCGGCCCCCACAGCGGCCAGGTCCCGGGCCGTGCGTTTGGTTTGATTCCGGTGGTCGGCGTCACGTTAGTGGTTACCACCTCGCGGTTGAACTGCGGGATCATCGCCTGATTCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCaa                                                                                                                                                                                                   <  1:1346/248‑1 (MQ=255)
                                               gtattcgcgCACAATAAGAGGTGATTGTTTACGTCGCTGCCGCAGACTAATCAACACAGAAGCAAAACAAACCAATGCCGCCAGCCCCCACAGCGGCCATTTGCCGGTACGTGCGTAGGGTGTGAGTCCGGTGGTCGGCGTCACTTTAGTGGTTACCCCCTCGGGGGTAAACGGGGGGATCTCCGCCTGAATCCCCCCCTGCGGGCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgc                                                                                                                                                                                           <  2:550350/243‑1 (MQ=255)
                                                 aTTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACTCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgc                                                                                                                                                                                         >  1:285983/1‑251 (MQ=255)
                                                       gcccGATAAGAAGGGATTATTTCCGTCGCTGACGCAGACTCACCAACACACCAGCAAAACCAAACAATGCCGGCAGCCCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGACCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGc                                                                                                                                                                                   <  1:53934/248‑1 (MQ=255)
                                                          gaaaaaaaagGGATTTTTCACGTCGCGCCCGCACACTCATAAACACAGCAGCAAAACAAAACAATGCCGACAGCCCCCAAGGAGGCGGTTTGCGGGTACGTGCGTGTGGTGTGATTCGGGTGGTCGGCGCCACGTTAGTGGTTAACCCCGCGGGGGGTACCGGTGGGATCCTCGCCGAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTcc                                                                                                                                                                                <  2:333892/244‑1 (MQ=255)
                                                             aaagaagGTATTATTAACGTGGCTGACGCAGACTCACCAACACAGCAGCAAACCCAAACAAGGCCGTCGGCCCCCACAGCGGCCAGTTGCCGGTACGGGCTTATGGTGTGAGGCCGGTGGTCGGCGTCACGTTTTGGGTTAACACCTCGCGGGTGAACTGCGGGAACACCGCCTAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGGGGGCGCGCCAGCTCCAgc                                                                                                                                                                             <  1:437736/250‑1 (MQ=255)
                                                                                   gtGACGCACACTCAACAACACAGCCGCAAAACCAAACAAGGCCGTCAGCCCCCACGGCGGCCAGTTGCCGTTATGTGGGTCTGGTGTTAGTCCGTTGGCCGGCGTCACTTTAGTGGTTAACCCCTCGGGGGTGAACTGGGGGATCATCACCTGATTCTCCCCTTTGGGGCAATAAGCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTTGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                                                                                                                       <  1:255377/250‑1 (MQ=255)
                                                                                   cggACGCAGCCTCATCAACACAGAAGAAAACCCAAAAAAGGCCGTCAGCACCCACAGGGGCCAGTTGCCGTTCTGTTCGTATGGTGTGAGTCCGGGGGTCGGCGTCACGTTAGTGTTTAACCCCGCGCGGGTGAACTGCGGGATCTTCGCCTGAATCTCCCCCTGCGGACCAATCACCGCCGTAATGCCTGTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                                                                                                                       <  1:528797/249‑1 (MQ=255)
                                                                                     aaCGCAGACCCACCAACACAGCAGCAAAACCAAACAAGGCGGTCATACCCCACAGCGGCCGGTGGCCGGTAGGTGCGTAGGGTGTGAGTCCGGGGGTCGGCGTCACGTTAGGGGTTAACCCCTCGCGGGGGACCGGTGGGATCATCGCCTGAATTTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAg                                                                                                                                                     <  1:105409/250‑1 (MQ=255)
                                                                                                  caaaaaacaacaaaaaCCAAACAGTGCCGTCAGCCCCCCCAGGGCCCAGTTGCGGGTAGTTGCGTAGGGTGTGAGCCCGGGTGCGGGTGTCACGTTAGGGGTTAACCCCCCGCGGGTGAACGGCGGGATCATCGCCGGACCCCCCCCCGGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAc                                                                                                                                        <  2:156545/239‑1 (MQ=255)
                                                                                                  caACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAGGTCTCGGGGGCCGCGGCATCATAGAGAATAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTTCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCCTGGAc                                                                                                                                        >  2:422213/1‑251 (MQ=255)
                                                                                                     cacaGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCa                                                                                                                                      >  1:91678/1‑250 (MQ=255)
                                                                                                        agcagcAAAGCCAAACAGGGCCGTCACCCCCCACGGCGGCAGTTTGCGGGTCGGTGCGTGTGGTGTGGGTCCGGGGGTCGGCGTCCCTTGGGGGGTAACCCCCCCGGGGGGGACGTGCGGGATAAGCCCCTGAACCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAg                                                                                                                                  <  2:397933/251‑1 (MQ=255)
                                                                                                             caAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAGATCCCGGTGGCCGCCGCATCATAGAGAATAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGACATCTGGAAGTGTTGCCATGGACCAATAGAttt                                                                                                                              >  2:485226/1‑250 (MQ=255)
                                                                                                             caAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTCGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGGCGGCGGCACGTTAGTGGTAAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGGTGGGGCTGCGCAACAGTGGGCGCGACAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGATGCCATGGACCAATAGAttt                                                                                                                              >  2:114884/1‑250 (MQ=255)
                                                                                                                    aaaCAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAAATCTCGGTGTCGGCGTTATGATAGTGGTTAACACCTCGAGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGGGGGCGCGCCCGCTCCAGCGCACGAATTCGCGCCATCTGGAAGTGTTGCCATGGGCCAATAGAGTTACAAAAcc                                                                                                                      >  2:284430/1‑251 (MQ=255)
                                                                                                                         aTGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGGATCTCGGGGGCGGCCGTATGATTAAAAAAAACACATCGCGGGTGAACTGCGGGATCATCGCATGAATATAACCCTGCGGGCCACTCACCGCCGGAATGCCGTTGTTGGGGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCAGTCGCGCCATCTGGAAGGGTTGCACTGGACCAATAGATTTACCAAACCACGcc                                                                                                                 >  2:437599/1‑250 (MQ=255)
                                                                                                                              gTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGtt                                                                                                            >  1:464804/1‑251 (MQ=255)
                                                                                                                                 aGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAAGTCCGGGGGGCCGCGTCACGATAGTAGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGAGTTGCCCTGGACCAATAGATTTCCCAAACCACGCATCGATgg                                                                                                          >  2:261487/1‑250 (MQ=255)
                                                                                                                                    aCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCa                                                                                                                                                                                     >  1:170032/1‑172 (MQ=255)
                                                                                                                                    aCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg                                                                                                          >  1:343538/1‑247 (MQ=255)
                                                                                                                                      ccacaagcggccgttggccgttaggTGCGTATGGTGTGGGTCCGGTGGTGGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGACCTGGGGGATCACCGCCTGAATCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATa                                                                                                     <  2:546914/227‑1 (MQ=255)
                                                                                                                                         caGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGGCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTc                                                                                                  >  2:36484/1‑250 (MQ=255)
                                                                                                                                                    tGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCg                                                                                      >  1:82717/1‑251 (MQ=255)
                                                                                                                                                                gggTATGGTGGGAGGCGGGTGGTGGGCGTCACGTTAGGGGTTAACCCCCCGCGGGTGAGCTGCGGGTTCGTCCCCGTACTCTCCCCCGGCGGCCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGc                                                                            <  2:372924/248‑1 (MQ=255)
                                                                                                                                                                   gtatgtgttgtgTCCGGGGGTCGGCTTCCCGTTGGTGGTAACCACCTCGCGGGTGAACTGGGGGATACTCGCCAGAATCTCACCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGa                                                                         <  1:261487/240‑1 (MQ=255)
                                                                                                                                                                     tgtggtgAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGGGAACTGGGGGATCATCGCCTGAATCTCACCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAg                                                                       <  1:263029/246‑1 (MQ=255)
                                                                                                                                                                        tgtgAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTAt                                                                   >  1:370170/1‑251 (MQ=255)
                                                                                                                                                                            gagcccGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAACCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATcgc                                                                <  1:66853/247‑1 (MQ=255)
                                                                                                                                                                             attCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgct                                                                                                                                                                                                         >  2:279279/3‑113 (MQ=255)
                                                                                                                                                                              ggtcgggtgtgcGGCGTCAGGTTGGTGGTAAACCCCTCGGGTGTGAGCTGGGGGATCCTCGCCTGAACCGCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATcgcg                                                               <  2:343538/240‑1 (MQ=255)
                                                                                                                                                                                   ggtggtCGGCGTCAGGTTAGTTTTTACCCCCTCGCGGGTCAACTGCGGAACCATCGCCTCACTCCCCCCCTGCGGCCCACTCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACtt                                                          <  2:524383/251‑1 (MQ=255)
                                                                                                                                                                                   ggaggtCGGCGTCACGTTCGTGGTTAACACCTCGCGGGTGAACTGCGAGATCATCGCCTGAATCTCACCCTGCGGGCAAATCACCCCCGTACTGCCGTTGTTGGTGCTGCGCAACAGCGGGCGCGCCGGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGTACCAATAAATTTAACAAACCACGCATCGTTGGAGATAGTCAGCCGATAGTCGGTATCCGGTCGGCAGTTCTCTCGCACtt                                                          >  2:27382/4‑251 (MQ=255)
                                                                                                                                                                                   agtggtCGTCGTCACCTTAGTGGGTAACACCTCGCGGGTGAAATGCGGGATCATCGCCTGAATCTCAACCTGCGGGCCAATCACCAACGTAATGCCGGTGGTGGTGCTGCTCAACAGGGGGCGCGACGGCTCCAGCGCACGCAGTCGAGCCATCTGGAAGTGTTGCAAGGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATGGTCGGTATCCGAGCGGAAGTTATCGTGCACtt                                                          >  2:400633/2‑251 (MQ=255)
                                                                                                                                                                                      ggtCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCt                                                       >  1:374762/1‑251 (MQ=255)
                                                                                                                                                                                        ggtcggcgccaggtgagggtttAACCCCGCGCGGGTGAACTGGGGGATCATCGCCTGAATCCCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGc                                                        <  1:134558/230‑1 (MQ=255)
                                                                                                                                                                                          ggCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGcc                                                   >  1:509221/1‑251 (MQ=255)
                                                                                                                                                                                            gcgccgcGTTGGTGGTTACCCCCTGGCGGTTGAACTGCGGGCTCAGCGCCTGAACCCCCCCCTGCGGGCAAATCGCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCg                                                  <  1:115108/245‑1 (MQ=255)
                                                                                                                                                                                                  cGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATg                                            >  1:504581/1‑251 (MQ=255)
                                                                                                                                                                                                      tAGTGTTTAACCCCGCGCGGGTGAACTGCGGGTTCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAt                                          <  1:117906/250‑1 (MQ=255)
                                                                                                                                                                                                       ggTGGTTACCCCCTCGGGGGTGAACTGCGGTATCGCCCCCTAAACCTCACCCTGGGGGCCAACCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAtct                                        <  2:298275/250‑1 (MQ=255)
                                                                                                                                                                                                        gTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGCTTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAtctc                                       >  2:92734/1‑251 (MQ=255)
                                                                                                                                                                                                        gTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGAtctc                                       >  1:387432/1‑251 (MQ=255)
                                                                                                                                                                                                         gggTTAACCCCTCGCGGGTGAACTGGGGGATCAGCGCCTGAATCCCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCg                                      <  1:538793/250‑1 (MQ=255)
                                                                                                                                                                                                             tAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGc                                  >  1:341357/1‑251 (MQ=255)
                                                                                                                                                                                                                             gAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTa                    >  1:442372/1‑250 (MQ=255)
                                                                                                                                                                                                                              atcTGCGGGATCATCGTCTGAATCTCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTaa                   <  2:134636/248‑1 (MQ=255)
                                                                                                                                                                                                                               aCTGCGGGTTCATCGCCTTATTTCCCCCCTGCGGGCCAATCGCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt                                                                                                      <  1:281803/166‑1 (MQ=255)
                                                                                                                                                                                                                               aCTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt                                                                                                      >  2:281803/1‑166 (MQ=255)
                                                                                                                                                                                                                                  gctggttcatccctggAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCt                <  1:18366/234‑1 (MQ=255)
                                                                                                                                                                                                                                    gggATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAAt            >  2:218946/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gggATCACCCCCTAAATCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAAt            <  2:341357/251‑1 (MQ=255)
                                                                                                                                                                                                                                       atccccgcgtgaaccccACCCGGCGGGCAAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGAt                                                                                           <  2:100270/153‑1 (MQ=255)
                                                                                                                                                                                                                                       atcatcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATACCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGAt                                                                                           >  1:100270/1‑169 (MQ=255)
                                                                                                                                                                                                                                          atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCAtt      >  2:292076/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCATTTGc   >  2:274708/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ccTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGCAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGTGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCATTTGcc  >  2:280989/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
GATAATGCACTGATGCCGCATCCGGTGAGCGTGCCGAAATATGGGATGTATTCCGGCACGATAAGAAGGGATTATTTACGTCGCTGACGCAGACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGCAGATAGTCGGTATCCGGGCGGAAGTTATCGCGCACTTGCTCGCCGAGAATGATCTCGTAGCAAATAGCCGCAGTAAGCTCAATACCATTTGCC  >  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: