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

GGGCAACAAACAGGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGC  >  NC_000913/4379388‑4379886
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                
gggCAACAAACAGAATTACGATGGTGGCATCCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGTTCCCATTTTTTCGTCTGTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGGTTTGTTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGAGCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGGTTCCGTGCCATTTGTCATGTCAAACAGCCCg                                                                                                                                                                                                                                                                  <  2:288908/251‑1 (MQ=255)
       aaacaGGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATc                                                                                                                                                                                                                                                            >  1:177588/1‑250 (MQ=255)
               ttACGATGGCAGCAACTACAGTTACCGCCCCGAGACTTTTAATAATTGGCTCTGATCCCATTTTTTCGTCTTTTACACTGATATTGGCCGCTTTCGGTGCCCGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGAAAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCCTTTTTCAGGGCCAACCGCCCGAAAATAAGTTCAATg                                                                                                                                                                                                                                                 >  2:212188/1‑251 (MQ=255)
                tACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATg                                                                                                                                                                                                                                                 <  1:324160/250‑1 (MQ=255)
                   gATGGTGGCAACCACAATTACCGCCCAAAGACTTTTGATAATTGGCTCTGGACCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCCTTTTTCAGGGCAAACAGCACGACAATCAGTTCAATGCTg                                                                                                                                                                                                                                              >  2:44686/1‑250 (MQ=255)
                         ggCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACcac                                                                                                                                                                                                                                       <  1:224976/250‑1 (MQ=255)
                                      aCCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAAccgc                                                                                                                                                                                                                    >  2:137240/1‑251 (MQ=255)
                                      aCCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGGCTGTTACAAGTATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTGTGTAAAAATTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAATCAGCCTGAAAATCAGTTCAATGCTGAACCACACAGCCGGAAccgc                                                                                                                                                                                                                    <  1:428175/251‑1 (MQ=255)
                                         gCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGt                                                                                                                                                                                                                 >  2:410238/1‑251 (MQ=255)
                                                       gATAATTGGCTCTGGTGCCATTTTTTTGTCTTTTACAATGATATTGGCAGCTTTCGGTGAAAGCTTAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCTCCAGATCGATCAGGTTAATGATCACGATAACCGGGTTTTGCAAAAAGTCGAAGAATCCCGCACAGGCTTCCGTCCCATTTTTAAGGGCAAACAGCACGAAAATAATTTCGATGCTGAACCACACATACGCAAccgccgcgacctttcccctct                                                                                                                                                                                                   >  2:260515/1‑236 (MQ=255)
                                                                          aTTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACgg                                                                                                                                                                                >  2:181347/1‑251 (MQ=255)
                                                                             ttttCGTCTTTTACACAGATATTTGCCGCCTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGATGCCGCCAGAGTGATCAGGTTAATGATCACGATCACCGGGTTTTGTAACAAGTCGACGAATACCGCCCAGGCTACCGGGCCATTTTTCAGGGCAAACAGCCCGATAATCAGTTAAATGCTGAACCACACCGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAACCGGCa                                                                                                                                                                              >  2:292992/1‑250 (MQ=255)
                                                                                                    gccggCTTCCGGTGCCAGTTAAACCCGGTTTTGGGGGGCAAGATGGGCTGCCGCAAGAGTGATCAGTTAAATGATCCCGATAGCCGGTTTTTCAAAAAAGTCGACGAACCCCGCCCGTCTTTCGGGGCCTTTTTTCAGGGCAAACAGCCCAAAAAGCAGTAAAAGGCTAACCCCCCCACCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCa                                                                                                                                                      <  1:285009/247‑1 (MQ=255)
                                                                                                              ggTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTAc                                                                                                                                            >  2:28761/1‑251 (MQ=255)
                                                                                                                       tCAAACCAGGGTTTTGGTGGGCAGAGGGGCTGCCGCAAGAGTGAGAAGTTAAATGATCACAATACCGGGGTTTTAAAAAAAGCCGACAACCCCCCCCCGGCCTCCCGGGCATTTTTTCGGGCCAAACGGCCCAAAAATAAGTTAAAGGCTAACCCCCACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGGCGTCATTGGCCGTACATACGGt                                                                                                                                     <  1:137240/250‑1 (MQ=255)
                                                                                                                          aaCCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACg                                                                                                                                >  2:62269/1‑251 (MQ=255)
                                                                                                                            cccggTTTTGGGTGTGAAGAAGAGCTGCCGCCAGAGTGACCAGTTAAAGGACCACAAAACCCGGTTTTTGAAAAAAGCCGCCAAACCCCCCCCAGCCTTCCGGCCATTTTTCCAGGGAAAACAGCCCAAAAATCAGTTCAATCCTGACCCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGt                                                                                                                               <  1:148231/248‑1 (MQ=255)
                                                                                                                              aGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCAACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGttt                                                                                                                             >  2:130245/1‑250 (MQ=255)
                                                                                                                                  ttttGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAACCAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtc                                                                                                                         >  2:169763/1‑251 (MQ=255)
                                                                                                                                     tGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcgtc                                                                                                                      >  1:371607/1‑251 (MQ=255)
                                                                                                                                     tGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcgtc                                                                                                                      >  1:307432/1‑251 (MQ=255)
                                                                                                                                      ggTGTGCAGCAGAGCTGAAGCCAGAGTGATCAGGCTAATGATAGCGATAACCGGGTTTTGTAAAAAGTCGAGGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGACAATCAGTTCAATGCTGAACTACACAGCCTGAACCGCCGTGCATTCGCGCAGCCTGTAAAAGCGATAAAACGGCACTTGTTTCCACCAGTTGTACGTGATTGGCCGTACATACGGTTTACGTTTAgtcgtc                                                                                                                      >  2:76207/1‑250 (MQ=255)
                                                                                                                                       gggtgtcaagaagggctcccgccAGAGTGATAAGGTAAATGACCACAATACCGGGGTTTTGAAAAAAGCCGACGAACCCCGCCCAGGCTGCCGGGCCTTTTTCAGGGGAAAACACCCCAAAAATAATTTCAATCCGAACCCCCACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAgtcgt                                                                                                                       <  1:169763/232‑1 (MQ=255)
                                                                                                                                                  agagCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGACCCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTc                                                                                                           >  2:435235/1‑251 (MQ=255)
                                                                                                                                                       tGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCtt                                                                                                        >  2:334406/1‑249 (MQ=255)
                                                                                                                                                                agTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTc                                                                                             >  1:166243/1‑251 (MQ=255)
                                                                                                                                                                    aTCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAggg                                                                                         >  2:11710/1‑251 (MQ=255)
                                                                                                                                                                      aaGGTTAATGATCACGATAACGGGTTTTTGAAAAAATGCGACGAACCCCCCCCAGCCTCCCGGCCCTTTTTTCAGGCCAAACAGCCCAAAAAGCAGTACAATGCTGAACCCCACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTc                                                                                       <  1:435235/250‑1 (MQ=255)
                                                                                                                                                                             aTGACCACGATAACGGGGTTTTGAAAAAAGTCGACGAACCCCGCCCAGGCTTCCGGCCAATTTTCCAGGGCAAACAGCCAAAAAACCATTTAAATGCGGAACCACACGGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATaa                                                                                <  1:181347/251‑1 (MQ=255)
                                                                                                                                                                              tGATCACGCTAACCCGGTTTTGTAAAAAGTCGCCGACTCCCGCCAAGCCTTCAGGGCCATTTTTCAGGGCAAACCGCCCGAACCTCCGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTCAAAGCGATAAAACGGCAATTTTTTCCACCAGGAGGACGTCATTGGCCGTCCATACGGTTTAAGTTTAGTCGTCATGTTGTACTCCTTAGCGTGGTTTCAGGGTCGCGATa                                                                                 >  2:182416/1‑249 (MQ=255)
                                                                                                                                                                                aTCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGaa                                                                                                                                                                                                                         <  2:343620/112‑1 (MQ=255)
                                                                                                                                                                                        aaCCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTc                                                                     >  2:179365/1‑251 (MQ=255)
                                                                                                                                                                                              ttgttaaaaaagcccgacgaACCCCGCCCAGGTTCCGGGACCATTTCCCGGGGCAACCAGCCCAAAAATAAATTCAATGTGAACCACCCCCGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACt                                                                <  2:112167/236‑1 (MQ=255)
                                                                                                                                                                                                attgaaaaaaGGCGACTACTCCCGCCCGGGTGTCCGGGCCTTTTTTCAGGCCAACGAGCCAAAAAATAATTTAAATGCTGACCCCAAAGACCGGAACCGTCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAAGTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACttt                                                              <  1:228086/245‑1 (MQ=255)
                                                                                                                                                                                                  aaaaaaaaGTCGACAACCCCCGCCCAGCCTCCCGGGCTGTTTTACAGGCCAAACACCCCAAAAATCCTTACAATGATAACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCTATAAGAAAGTCTTTCGAACTTTCt                                                            <  2:305956/247‑1 (MQ=255)
                                                                                                                                                                                                    aaaaaagccaccaaaccccGCCCAGGCTTCCGGCCATTTTTCCAGGGAAACCACCCCAAAAACCATTCAAAGGCTAAACCACACAGCCGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACt                                                         <  1:28761/235‑1 (MQ=255)
                                                                                                                                                                                                              acgaATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct                                               >  1:216785/1‑251 (MQ=255)
                                                                                                                                                                                                               ccgacccccgccAGGCTGCCGGGCCATTTTTACGGGCAAACCCCCCAAAAAAGAGTTACAATGCGGACCACACCGGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctg                                                 <  2:236253/239‑1 (MQ=255)
                                                                                                                                                                                                                 aaTCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGa                                                                   >  1:16742/1‑228 (MQ=255)
                                                                                                                                                                                                                 aaTCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGa                                             >  1:375806/1‑250 (MQ=255)
                                                                                                                                                                                                                    cccccccccgggctCCGGGCCTTTTTTCAAGGCCAAACGCCCCAAAAATCAGTTAAAGGTTAACCCACACAGCGGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAAt                                           <  1:11710/238‑1 (MQ=255)
                                                                                                                                                                                                                         ccAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAACTCAGTTCAATGCTGACCCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACt                                                                >  2:469832/1‑223 (MQ=255)
                                                                                                                                                                                                                                    cGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCg                           >  2:261704/1‑251 (MQ=255)
                                                                                                                                                                                                                                             atttgagcaggCAAACACCCCCAAAAACAGTTCAAAGGCGACACCACACGCCGGAACCCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTGCATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACgtg                      <  2:255334/242‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGc  >  2:9186/1‑249 (MQ=255)
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                
GGGCAACAAACAGGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTTGATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGGTGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTTGTAAAAAGTCGACGAATCCCGCCCAGGCTTCCGGGCCATTTTTCAGGGCAAACAGCCCGAAAATCAGTTCAATGCTGAACCACACAGCCGGAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAGACTTCGGAGC  >  NC_000913/4379388‑4379886

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