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

GGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTC  >  NC_000913/1822052‑1822345                                                                                                                                                    |                                                                                                                                                     ggCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAAccc                                                                                                                                                   <  2:53941/149‑1 (MQ=255) ggCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCg                                                                                                                                                  <  1:260763/150‑1 (MQ=255)   cAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCa                                                                                                                                               <  2:807542/151‑1 (MQ=255)   cAAGTACCGCAAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCa                                                                                                                                               <  1:1046656/151‑1 (MQ=255)    aaGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGcc                                                                                                                                                >  2:949838/1‑149 (MQ=255)          cgcgAGCATGGATAATGCTCACTTCAGACGAAGGTTTACCCGCAATGTGTGCAACCTGCGCCCGGC‑TCATGCTTTTTTTGACGTCTTTCCCCACAGGCTGTACAAACTGGTTTTTATTACGATCATAAC‑CCGTACAACCCGCCAGc                                                                                                                                             <  1:543228/146‑1 (MQ=255)            cgAGCATGTATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAACCTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTa                                                                                                                                      <  1:628352/151‑1 (MQ=255)               gCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAAt                                                                                                                                    >  1:141498/1‑150 (MQ=255)                cATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATTATAAG‑CCGTACAACCCGCCAGCATGGTTAATa                                                                                                                                   <  1:324853/150‑1 (MQ=255)                 aTGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATAcc                                                                                                                                 >  1:224278/1‑151 (MQ=255)                   ggATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATAccg                                                                                                                                <  1:288933/150‑1 (MQ=255)                    gATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATAccgcc                                                                                                                              >  2:1030359/1‑151 (MQ=255)                          gCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGc                                                                                                                         >  1:973751/1‑150 (MQ=255)                          gCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGc                                                                                                                         >  1:806608/1‑150 (MQ=255)                              aCTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTc                                                                                                                     >  1:580483/1‑150 (MQ=255)                               cTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTc                                                                                                                     >  1:1036121/1‑149 (MQ=255)                                      cGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCt                                                                                                            >  2:125274/1‑151 (MQ=255)                                        aaGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGc                                                                                                          >  1:219382/1‑151 (MQ=255)                                        aaGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGc                                                                                                          >  1:792593/1‑151 (MQ=255)                                        aaGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTAAAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCt                                                                                                            <  1:648596/149‑1 (MQ=255)                                          gttttCCCCGCATTTTGTGCACCCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTCCCCCCCAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCa                                                                                                        <  1:65418/149‑1 (MQ=255)                                            tttACCCGCAATTTGTGCACCCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCat                                                                                                       <  2:428161/150‑1 (MQ=255)                                               ccccGCATTTTTTGCACCTGGCGCCCGCC‑TCATGCCTTTTTTGGCTTTCTCCCCCACAGCTGGAACAAACTGTTCTTTGGTCCGACAATAAC‑CCTTCCAACCCGCCAGCAGGTTTAATACCGCCGCGGCACTCAGAATTCCTGCCATAtt                                                                                                    <  1:584864/149‑1 (MQ=255)                                                cccGCAATTTGTGCACCCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCCCCCCAGGCTGAACAAACTGGTCTTTGGTACGATAAAAAG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTc                                                                                                   <  1:522232/150‑1 (MQ=255)                                                 cccaattttGTCAAACCTGC‑CCCGGCTTCTTGCCTTTTTTGACGTCTTTCACCACAGGCTTTAAAAACGGTTCTTTGGTCCGATCATAAG‑CCGTACACCCCGCCATCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCatat                                                                                                     <  1:810132/141‑1 (MQ=255)                                                 ccGCAATTTGTGCACCCTGCGCCCGGC‑TCATGCTTTTTTTGCCGTCTTTCCCCACAGGCTGTACAAACGGGTCTTTGGTACGATCATAAG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCtt                                                                                                 <  2:731396/151‑1 (MQ=255)                                                 ccGCAATTTGTGCACCCTGCGCCCGGC‑TCATGCTTTTTTTGACGTCTTCCCCCACAGGCTGTAAAAACTGGTCTTTGGTCCGATCATAAG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCtt                                                                                                 <  1:368721/151‑1 (MQ=255)                                                    cAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGtt                                                                                              >  2:851045/1‑151 (MQ=255)                                                     aaTTTGTCAACCCTGCGCCCGTC‑TCTTGCCTTTTTGGCCGTCTTTCCCCAAAGGCTGTACAAACTGGTCTTTGTTACGATCAAACG‑CCGTACACCCCGCCAGCTTGGTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTg                                                                                                <  1:238061/148‑1 (MQ=255)                                                        ttGTGCACCCTGCGCCCGCC‑TCATGCCTTTTTTGACGTCTTTCCCCACAGGCTGTACAACCGGGTTTTTGGTACGATCATAAG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCat                                                                                           <  2:372105/150‑1 (MQ=255)                                                           tgCACCCGGCCCCCGGC‑TAATCCCTTTTTTGCCGCTTTTCACCACAGGCTGTAAAACCTGGTTTTTGGTACGATCATAAG‑CCTTCCACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCata                                                                                          <  2:156875/148‑1 (MQ=255)                                                           tgCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTc                                                                                       >  1:895258/1‑151 (MQ=255)                                                              acccggcgcCCGGC‑TCTTCCTTTTTTGGACGTCTTCCCCCACAGCCTGTCAAAACGGTTTTTGGGTCGGATAATAAC‑CCGAACACCCCGCCAGCGTGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTc                                                                                       <  1:736949/143‑1 (MQ=255)                                                              acccggcccccGCC‑TCTTCCCTTTTTTGCCTTTTTTCCCCACAGGTTGTAAAAACTGGTTTTTGGTACGATCATAGG‑CCGTACACCCCCCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCg                                                                                     <  1:677861/142‑1 (MQ=255)                                                               aCCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTc                                                                                   >  1:474672/1‑151 (MQ=255)                                                                    cgcCCGGC‑TCATCCCTTTTTTGCCGCCTTCCCCCACAGGCTGTACAACCTGGTCTTGGGTACGATCATAAG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCtt                                                                              <  1:415204/151‑1 (MQ=255)                                                                      cccGGC‑TCTTCCTTTTTTTGACGTCTTCCCCCACGGGTTGAACAACCTGTTTTTGGGTCCGATAATAAG‑CCGTCCACCCCGCCAGTATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTg                                                                             <  1:110329/150‑1 (MQ=255)                                                                             tCATGCTTTTTTGGCCGTCTTTCCCCACAGGCTGTAAAACCGGGTTTTGGGTACGATCAAAAG‑CCGTACACCCCGCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTAt                                                                        <  1:1013011/149‑1 (MQ=255)                                                                             cccttccTTTTTTAGCTTCTTAACCCTCAGGGTGACCAACCTGGTTTTTGTTACGACCATAAG‑CCGTACACCCCGCCAACATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCctc                                                                                <  2:747952/136‑1 (MQ=255)                                                                              cATGCTTTTTTTGACTTCTTCCACCAAGGCCTGTACAAACTGGTTTTTGGTCCGATCAAAAG‑CCTTACACCCCGCCAGCATGTTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATc                                                                       <  1:271722/149‑1 (MQ=255)                                                                                  ctttttttGCCGTCTTTCCCCACAGGCTGTACAACCTGGTTTTTGGTACGATCAAAAC‑CCGTAAACCCCGCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAg                                                                     <  2:535864/145‑1 (MQ=255)                                                                                   cTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCgtgt                                                                >  2:705145/1‑151 (MQ=255)                                                                                    ttttttGACTTCTTTCACCACAGGCTGTACAACCGGGTTTTTGGTACGATCAAAAG‑CCGTACACCCCCCCAGCATGTTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCg                                                                   <  2:270634/147‑1 (MQ=255)                                                                                       tttGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGtt                                                               >  2:10790/1‑148 (MQ=255)                                                                                              ttttcccccacAGGCTTTAAAACCTGTTTTTGGGTCCGATCAAACG‑CCTACCACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGAt                                                           <  1:573242/138‑1 (MQ=255)                                                                                              tCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCt                                                     >  1:144598/1‑151 (MQ=255)                                                                                              tCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCt                                                     >  1:870463/1‑151 (MQ=255)                                                                                                tttCACCACAGCTTGTACAACCGGGTTTTTGGTACGACAATAAG‑CCTTAAACCCCGCCAGAATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGAt                                                           <  2:196875/143‑1 (MQ=255)                                                                                                 tcccccacAGCCTGTACAACCTGGTTTTTGGTACGATCATAGC‑CCGTACACCCCGCCAGCTTGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTgg                                                   <  1:257078/146‑1 (MQ=255)                                                                                                 tccaccacAGGCTGTACAACCTGTTTTTTGGTACGATCATAAG‑CCTTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCt                                                     <  2:17733/146‑1 (MQ=255)                                                                                                  tCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGa                                                  >  2:1048925/1‑150 (MQ=255)                                                                                                          ggCTGTACAAACTGGTTTTTGGTACGATCATAAG‑CCGTACACCCCGCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTg                                          <  2:926192/150‑1 (MQ=255)                                                                                                            cTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGgc                                       >  1:186099/1‑151 (MQ=255)                                                                                                            cTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGgc                                       >  1:935695/1‑151 (MQ=255)                                                                                                                aCAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCtt                                   >  2:93874/1‑151 (MQ=255)                                                                                                                aCAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCtt                                   >  1:907181/1‑151 (MQ=255)                                                                                                                       gtttttgggTACG‑TCAAAAGCCCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAg                            <  1:992008/144‑1 (MQ=255)                                                                                                                       ggTTTTTGGTACGTTAATAAG‑CCTTACATCCCGCCGGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTc                               <  2:208956/148‑1 (MQ=255)                                                                                                                       ggTCTTTGGTACTATAAAAAG‑CCTTACACCCCGCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTc                               <  2:421538/148‑1 (MQ=255)                                                                                                                           tttGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCg                        >  2:816734/1‑151 (MQ=255)                                                                                                                           tttGGTACGATCAAAAC‑CCGTCCACCCCCCCAGCTTGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAg                            <  1:457820/147‑1 (MQ=255)                                                                                                                            ttGGTACGATCATAAG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGcc                          <  2:914796/148‑1 (MQ=255)                                                                                                                                tACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAAc                   >  2:111109/1‑151 (MQ=255)                                                                                                                                 aCGATCATAAG‑CCTTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAAc                   <  2:410978/150‑1 (MQ=255)                                                                                                                                      tcaaaaccggTACACCCCGCCAGCATGTTTAATCCCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGaaa               <  2:660147/141‑1 (MQ=255)                                                                                                                                       aTAAG‑CCTTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACgg            <  2:1051151/151‑1 (MQ=255)                                                                                                                                       aTAAG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGaaa               <  2:623131/148‑1 (MQ=255)                                                                                                                                        tAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGc           >  1:136437/1‑151 (MQ=255)                                                                                                                                         aaG‑CCGTACACCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCt          <  2:578506/151‑1 (MQ=255)                                                                                                                                             ccaaacacCCCGCCAGCATTGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTcc        <  1:949700/146‑1 (MQ=255)                                                                                                                                             ccGTACACCCCGCCAGCTTGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCg       <  1:435629/151‑1 (MQ=255)                                                                                                                                              cGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGc      >  1:285758/1‑151 (MQ=255)                                                                                                                                                 acacCCCGCCAGCATGTTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCt     <  2:948958/149‑1 (MQ=255)                                                                                                                                                   aaCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTc  >  2:501746/1‑150 (MQ=255)                                                                                                                                                   aaCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTc  >  1:576666/1‑150 (MQ=255)                                                                                                                                                   aaCCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTc  >  1:278516/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGC‑TCATGCCTTTTTTGACGTCTTTCACCACAGGCTGTACAAACTGGTCTTTGGTACGATCATAAG‑CCGTACAACCCGCCAGCATGGTTAATACCGCCGCTGCACTCAGAATTCCTGCCATATTCTTGTTCATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTC  >  NC_000913/1822052‑1822345

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: ‑