Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,822,197	0	A	C	43.1%	101.7 / 14.9	65	G20G (GGT→GGG)	osmE	osmotically‑inducible lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (22/15);  new base C (0/28);  total (22/43)
Fisher's exact test for biased strand distribution p-value = 9.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCT  >  NC_000913/1822047‑1822342                                                                                                                                                        |                                                                                                                                                  ggTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACaa                                                                                                                                                   <  2:165345/151‑1 (MQ=255)     tGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCg                                                                                                                                               >  1:542739/1‑151 (MQ=255)      ggCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAAccc                                                                                                                                                <  1:198985/149‑1 (MQ=255)      ggCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCg                                                                                                                                               <  1:734214/150‑1 (MQ=255)      ggCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGc                                                                                                                                              <  1:475676/151‑1 (MQ=255)        cAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCa                                                                                                                                            <  1:1020992/151‑1 (MQ=255)           gTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCa                                                                                                                                         <  2:37970/151‑1 (MQ=255)               cgcgAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGt                                                                                                                                     >  2:653272/1‑151 (MQ=255)               cgcgAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGt                                                                                                                                     >  2:371631/1‑151 (MQ=255)                gcgAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGtt                                                                                                                                    <  1:1107531/151‑1 (MQ=255)                  gAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTa                                                                                                                                   >  1:171285/1‑150 (MQ=255)                   aGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAAt                                                                                                                                 <  1:790093/151‑1 (MQ=255)                   aGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAAt                                                                                                                                 >  2:142121/1‑151 (MQ=255)                     cATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATAc                                                                                                                               <  1:80462/151‑1 (MQ=255)                      aTGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATAcc                                                                                                                              >  1:565884/1‑151 (MQ=255)                        ggATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATAccgc                                                                                                                            <  2:605600/151‑1 (MQ=255)                          aTCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATAccgcc                                                                                                                           <  1:649986/150‑1 (MQ=255)                                    cTTCAGACGAAGGTTTCCCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCCCCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCa                                                                                                                 <  1:93136/150‑1 (MQ=255)                                       cAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCCCAGGCTGTACAAACTGGTTTTTGGTACGATCATAAGCCTT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAAt                                                                                                             <  1:627117/151‑1 (MQ=255)                                             aaGGTTTCCCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACTTCTTTCACCACAGGCTGTACAAACGGGTTTTTGGTACGATCATAACCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTcc                                                                                                          <  2:139222/148‑1 (MQ=255)                                                  ttccccGCATTTTGTGCACCCTGCGCCCGGCTCATGCTTTTTTTGCCGTCTTCCCCCACGGTCGGTAAAAACTGGTTTTTGGTCCGATCATAAGCCGT‑ACACCCCGCAAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCata                                                                                                   <  2:1102887/147‑1 (MQ=255)                                                  ttACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCatat                                                                                                  >  1:732384/1‑151 (MQ=255)                                                  ttACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCatat                                                                                                  >  1:50794/1‑151 (MQ=255)                                                       cGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTg                                                                                             >  2:507664/1‑151 (MQ=255)                                                             ttGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGCCGTCTTTCACCACAGGCTGTACAAACTGGTTTTTGGTACGATCATAAGCCTT‑AAACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCat                                                                                        <  1:442543/150‑1 (MQ=255)                                                                tgCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTc                                                                                    >  2:828090/1‑151 (MQ=255)                                                                  cAACCTGCGCCCGGCTCATGCTTTTTTGGACGTTTTTCCCCACAGGCTGTACAAACGGGTTTTTGGTACGATCATAAGCCGT‑ACACCCCGCCGGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCg                                                                                  <  1:728652/151‑1 (MQ=255)                                                                  cAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCg                                                                                  >  1:966558/1‑151 (MQ=255)                                                                   agccTGCGCCCGGCTCATGCCTTTTTTGACGTTTTCCACCAAAGGCTGTAAAAACTGGTCTTTGGTACGATCAAAAGCCGT‑ACAACCCGCCGGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGt                                                                                 <  1:295627/149‑1 (MQ=255)                                                                   acccTGCGCCCGCCTCTTCCTTTTTTTGACGTTTTCCCCCACAGGCTGTAAAAACTGTTTTTTGGTCCAATCAAAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTcc                                                                                   <  2:254127/147‑1 (MQ=255)                                                                   aaCCTGCCCCCGGCTCACGCTTTTTTTGACTTTTTCCACCACAGGTTGTAAAACCGGTTTTTTGGTCCGATAATAACCCTA‑AAACCCCCCCAGCTTGGTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTc                                                                                    <  2:682613/148‑1 (MQ=255)                                                                     ccgccccccGGCTCTTGCCTTTTTAGCCGTCTTTCCCCACGGCCTATAAAACCTGTTTTTCGTCACGA‑CAAAAGCCATCAAACCCCGCAACAATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTcc                                                                                   <  2:493264/141‑1 (MQ=255)                                                                     ccTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTcc                                                                               >  1:1069110/1‑151 (MQ=255)                                                                         cgcCCGGCTCATGCTTTTTTTGACGTCTTCCCCCACAGGCTGTACAAACTGGTTTTTGGTACGATCATAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCctc                                                                             <  1:503072/149‑1 (MQ=255)                                                                         cgcCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCtt                                                                           >  2:1007095/1‑151 (MQ=255)                                                                           cccgtttCATGCTTTTTTTGCCGTCTTTCCCCCCAGGCTGTACAACCTGGTCTTTGGTACGATCATAAGCCTT‑CCAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGt                                                                         <  1:857379/145‑1 (MQ=255)                                                                           cccGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCtt                                                                           <  1:252992/149‑1 (MQ=255)                                                                           cccGCCTCATCCTTTTTTTGCCTTTTTCCCCCACGGGCTGAAAAAACTGGTTTTTGTTCCGATCAAAACCCGT‑CCACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCatat                                                                                      <  2:606348/138‑1 (MQ=255)                                                                                       cTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCgtgt                                                             >  1:391570/1‑151 (MQ=255)                                                                                        ttttttGACGTCTTTCCCCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCgtgt                                                             <  1:343051/150‑1 (MQ=255)                                                                                         tttttGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGtt                                                            >  2:585267/1‑150 (MQ=255)                                                                                          ttttGACGTTTTCCCCCACGGCTTTTACAAATTGTTTTTTGTTACTACCATAGGCCTT‑ACACCCCGCCAGCATGGTTAATCCCGCCGCTCCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTa                                                           <  2:1074511/150‑1 (MQ=255)                                                                                          ttttGACGTCTTTCACCACAGGCTGTACAACCGGGTCTTTGGTACGATCATAAGCCTT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTa                                                           <  1:33524/150‑1 (MQ=255)                                                                                           tttGACGTCTTTCCCCACAGGCTGTAAAAACTGGTTTTTGGTACGATCATAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAg                                                          <  1:465014/150‑1 (MQ=255)                                                                                           tttGACGTCTTTCACCACAGGCGGTACAAACTGGTTTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAg                                                          <  1:98846/150‑1 (MQ=255)                                                                                            ttgccttCTTTCACCGCAGGTTGTACAAACTGGTTTTTGGTACGATCATAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGa                                                         <  1:31838/144‑1 (MQ=255)                                                                                             tGACGTCTTTCACCACAGGCTGTACAAACTGGTTTTTGGTACGATCATAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGAt                                                        <  1:46443/150‑1 (MQ=255)                                                                                               aCGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATa                                                       >  2:879641/1‑149 (MQ=255)                                                                                                    ttcccccacAGGCTGTAAAACCTGGTTTTGGGTACTATCAAAACCCGA‑ACACCCCCCAAGCATGTTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATaa                                                      <  2:893041/140‑1 (MQ=255)                                                                                                       caccacAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAAt                                             >  2:898862/1‑151 (MQ=255)                                                                                                            caGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACAtt                                        >  2:32389/1‑151 (MQ=255)                                                                                                             aGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTg                                       >  2:146484/1‑151 (MQ=255)                                                                                                              gtctgaaaaaccTGGTTTTTGGTACGTTCAAAAGCCGT‑ACACCCCGCCAGCATGGTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGaa                                              <  2:602034/132‑1 (MQ=255)                                                                                                                      aaccggggttttGTTACCATCAAAGGCAAT‑ACAACCCCCAACCATGGTGAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGgc                                    <  1:946300/136‑1 (MQ=255)                                                                                                                      aaccTGGTCTTTGGTACGATCATAACCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGgcg                                   <  1:317400/143‑1 (MQ=255)                                                                                                                      aaattgtttttgggTACGATCAAAAGCCGA‑ACACCCCGCAAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCttt                               <  2:466570/138‑1 (MQ=255)                                                                                                                          tGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTc                            >  1:558112/1‑149 (MQ=255)                                                                                                                             tCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGcc                       >  2:352884/1‑151 (MQ=255)                                                                                                                               tttgttcgattCAAAACCCGT‑CCACCCCGCCAGCATGTTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGccc                      <  2:1118513/140‑1 (MQ=255)                                                                                                                               tttGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCg                     >  1:859452/1‑151 (MQ=255)                                                                                                                               ttggttACGATCAAACGCCTT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGccc                      <  1:389704/145‑1 (MQ=255)                                                                                                                                  ggTACGATAAAAACCGTT‑ACACCCCGCCACCATGTTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCg                     <  2:229917/148‑1 (MQ=255)                                                                                                                                     aCGATCATAACCCGT‑ACACCCCGCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGaa                 <  2:161851/149‑1 (MQ=255)                                                                                                                                      cGATCAAAAGCCGT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAAc                <  1:765371/149‑1 (MQ=255)                                                                                                                                           aaaaGCCTT‑ACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACg               <  2:504274/143‑1 (MQ=255)                                                                                                                                            aaaccgtt‑ACACCCCCCCAGAATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGc        <  2:774167/144‑1 (MQ=255)                                                                                                                                                     acacCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCt  <  1:355185/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  GGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAGCCGT‑ACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCT  >  NC_000913/1822047‑1822342

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