Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,494	0	T	G	46.8%	101.8 / 14.7	62	intergenic (‑512/+91)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/32);  new base G (29/0);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 7.32e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCA  >  NC_000913/987346‑987643                                                                                                                                                     |                                                                                                                                                      tcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGt                                                                                                                                                       <  1:707915/149‑1 (MQ=255)      cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCaaa                                                                                                                                                 >  1:59163/1‑150 (MQ=255)       aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATa                                                                                                                                               <  1:154378/151‑1 (MQ=255)            tCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCaaaaa                                                                                                                                               >  2:359416/1‑144 (MQ=255)              cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAAttttt                                                                                                                                        >  2:443769/1‑148 (MQ=255)                       aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaac                                                                                                                               <  1:442621/151‑1 (MQ=255)                       aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAAAAAATTTCTTTACttta                                                                                                                                >  1:181309/1‑148 (MQ=255)                        aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGATTGTTTGGCttta                                                                                                                                >  1:470200/1‑147 (MQ=255)                          aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAATTGCtttaatttaaaa                                                                                                                             >  2:347990/1‑140 (MQ=255)                             tAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATGCTTTACTTaacaaaaa                                                                                                                          >  1:437328/1‑148 (MQ=255)                                 ctGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGTACAAAGATTGCTGGATTTTTCAAAAAGa                                                                                                                        >  1:390456/2‑148 (MQ=255)                                 aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATaa                                                                                                                     <  2:325929/151‑1 (MQ=255)                                 aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATa                                                                                                                      <  2:570822/150‑1 (MQ=255)                                   gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGCAAAAAGAATGGGTGAGTTAAAAAAAAGaaaaa                                                                                                                    >  1:902266/1‑149 (MQ=255)                                    aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTCCAAATAGATTTCTTTTATTAACAAACAAATaaaat                                                                                                                  >  1:706699/1‑148 (MQ=255)                                     aaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTg                                                                                                                 <  1:164613/151‑1 (MQ=255)                                          gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAg                                                                                                             <  2:306032/150‑1 (MQ=255)                                          gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAg                                                                                                             <  2:718720/150‑1 (MQ=255)                                          gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGa                                                                                                            <  1:348017/151‑1 (MQ=255)                                          gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGa                                                                                                            <  1:295118/151‑1 (MQ=255)                                          gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCGAATAAATTGATTTATTTTAATAACAAATAAGtt                                                                                                                  >  1:475519/1‑145 (MQ=255)                                            gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaa                                                                                                           <  2:719072/150‑1 (MQ=255)                                                 tAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAt                                                                                                      <  2:49251/150‑1 (MQ=255)                                                   gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTATGTGTTGCGCGTAAATGTTGCAATGTAATATTCTTTTACATTATTATCGTTATTTCTTTTTTTCAATTATTTTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaac                                                                                                   <  1:868574/151‑2 (MQ=255)                                                   gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGAAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaa                                                                                                    <  2:105841/150‑1 (MQ=255)                                                     tAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGTAATAAAATGTTTGGCTTAAAAAAAAAAAAACATGAGAGAAAAAATaaaa                                                                                                  >  2:418309/1‑150 (MQ=255)                                                       aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTCTTTTCTTTAAAAAAAAAAAACTTTAAAGAAAAAATaaaaaa                                                                                                >  2:271482/1‑150 (MQ=255)                                                         aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAaga                                                                                             <  1:499172/151‑1 (MQ=255)                                                          aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTTTTTACTTTAAAAACAAAAAAATTTTCAGaaaaaaaaaaaaaaaa                                                                                             >  2:252257/1‑148 (MQ=255)                                                            aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTCTTTACTTTACAAAAAAAAAAAATTACAGAAAAAAAAAAAAAaaagac                                                                                          >  1:412032/1‑149 (MQ=255)                                                                   tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGGAAAAAAATTTGTTATTTTAAAAAAAAAAAAAATTGACTGAAAAAAAAAAAAAAAAACCa                                                                                         >  2:181493/1‑145 (MQ=255)                                                                    ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                   <  1:143002/150‑1 (MQ=255)                                                                              aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                        <  1:800144/151‑1 (MQ=255)                                                                              aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                        <  2:191155/151‑1 (MQ=255)                                                                              aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGAATGCTTTGCTTAACAAAAAAATAACCTGATAGAAAAGAAAAAAAAAAAGCAAACGGGGGCTgg                                                                              >  1:80131/1‑145 (MQ=255)                                                                              aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTCTTTACTTAAAAAAAAAAAAAAATGACAGAAAAAAAAAAAAAAAAAACAACGGGGGCTgggcctt                                                                         >  2:281716/1‑145 (MQ=255)                                                                                  cGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                    <  2:610673/151‑1 (MQ=255)                                                                                    gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACg                                                                  <  1:148830/151‑1 (MQ=255)                                                                                    gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATGCTTTAATTAAAAAAAAAAAAAATTGACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGCTCTTGAAGa                                                                    >  2:713183/1‑149 (MQ=255)                                                                                     aTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGc                                                                 <  1:444545/151‑1 (MQ=255)                                                                                      ttAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAGATTGGTTGACTTAACAAAAAAAAAACTTTACAGAAAAAAAAAAAAAAAAACCAGCGGGGGGTGGGTCtt                                                                         >  1:815621/1‑142 (MQ=255)                                                                                        aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                               <  1:6771/150‑1 (MQ=255)                                                                                           ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                            <  1:27453/150‑1 (MQ=255)                                                                                              tGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAg                                                         <  2:896523/150‑1 (MQ=255)                                                                                              tGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTGTTTACTTTAAAAAAAAAAAACTTTACAGAAAAAATAAAAAAAAAAACAGCGGGGGCGGGCTCTTTGAGAAGCTaaa                                                             >  1:392139/1‑146 (MQ=255)                                                                                               gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATTCTTTACTTTAAAAAAAGAAAAATTTGCAAAAAAAATAAAAAAAAAAACAAGGGGGGCCGGCTCCTGAAGACGCCAAAATa                                                          >  1:365800/1‑148 (MQ=255)                                                                                                     aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                 <  2:155409/151‑1 (MQ=255)                                                                                                         gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                             <  1:660904/151‑1 (MQ=255)                                                                                                         gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                             <  1:359445/151‑1 (MQ=255)                                                                                                         gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGCTTGACTTTAAAAACAAAAAAAATTACAGAAAAAAAAAAAAAAAAAACAGCGGGGGCCGGCTCTTGTAGAAGCGAAAATAAAAACTgggggt                                              >  1:509744/1‑146 (MQ=255)                                                                                                           tGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAc                                            <  2:458266/150‑1 (MQ=255)                                                                                                              aaTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGTGAAAAAATTTGTTTACTTTAAAAAAAAATAACTTGAAAGGAAAAAAAAAAAAAAAAACAGCGGGCGGCGGCCCTTTAAGAAGCTAAAATAAAAAAAGGCGttt                                             >  2:726685/1‑145 (MQ=255)                                                                                                               aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTTTTGACTTAAAAAAAAAATAAATTGACAGAAAAGAAAAAAAAAAAGGCAACGGGGGCGGGCTCTTTAAGACGCAAAATTAAAAACTTGCGtttc                                            >  1:135367/1‑144 (MQ=255)                                                                                                                  tAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCCAAAAGATTGGTTGACTTAAAAAACAGATAAATTTACAGAAAAAATAAAAAAAAAACCAGCGGGGGCTGGGTCTTTAAGACGCCAAAATAGAAGCTGGCGGTAcgc                                          >  1:963559/1‑145 (MQ=255)                                                                                                                    atatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                  <  1:141672/151‑1 (MQ=255)                                                                                                                       ttCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGAATGCTTTACTTAAAAAAAAAAAAAATTTACAGAAAAGATAAAAAAAAAAACAGCGGGGGCTGGCTCTTTAAAACGCCAAAATAGAAACTGGCGGTACGCCGAGTACGTgg                               >  1:335923/1‑151 (MQ=255)                                                                                                                         ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                              <  1:978161/150‑1 (MQ=255)                                                                                                                           gTAACAGGATGATCGTTATGGCGGGGGTGCAAATAGATTGGTTGACTTAACAAAAAGATAACTTTACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTCGCTCTTTAAGAAGCTCAAATAAAAACTGGGGCTACACGGAGGACGTggg                              >  2:480680/1‑148 (MQ=255)                                                                                                                             aaCAGGATGATCGTTATGGCGGGGGTTCCAATAAATTGCTTGACTTTAAAAACAAAAAAATTTACAGAAAAAAAAAAAAAAAAACCAGCGGGGGCTGGCCCTTTTAAAAGCTAAAATAAAAACTGGGGGTACGCGGGGGACGgggggaagg                         >  1:1000575/1‑146 (MQ=255)                                                                                                                             aaCAGGATGATCGTTATGGCGGGGGTGGAAAAAGATTTGTTTACTTTAAAAAAAAAAAAATTGGAAGAAAAAAAAAAAAAAAAGACAACGGGGGGGGGCGCTTTTAGAAGCTAAAATAAAAAATGGGGTTTTGCGGAGGAAGTGGGGACgg                         >  2:839241/1‑151 (MQ=255)                                                                                                                              aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGa                        <  1:729473/151‑1 (MQ=255)                                                                                                                                   atgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcac                   <  1:658388/151‑1 (MQ=255)                                                                                                                                                    gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa  <  2:171072/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCA  >  NC_000913/987346‑987643

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