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

CTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCA  >  NC_000913/987347‑987643                                                                                                                                                    |                                                                                                                                                        cttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTACAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTg                                                                                                                                                       >  2:51498/1‑150 (MQ=255)   ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑c                                                                                                                                                     <  2:956958/149‑1 (MQ=255)     cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑Gaaa                                                                                                                                                  >  2:1118735/1‑150 (MQ=255)      aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGtttgaaa                                                                                                                                                   >  4:55542/1‑144 (MQ=255)       gACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAg                                                                                                                                               <  3:9538/151‑1 (MQ=255)       gACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑Caaaaa                                                                                                                                                >  3:167452/1‑148 (MQ=255)              aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGAGTTG‑Caaaaa                                                                                                                                                >  1:277973/1‑141 (MQ=255)                    ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CGAATAAATTGCTTGGCtt                                                                                                                                   >  1:497508/1‑150 (MQ=255)                       aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAAAAGATTGCTTGACtt                                                                                                                                   >  1:485433/1‑147 (MQ=255)                        aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTaacaa                                                                                                                              <  2:333704/151‑1 (MQ=255)                        aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CAAAAAAATTTGTTTACTTa                                                                                                                                  >  2:1132653/1‑147 (MQ=255)                             aGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAGATTTTTTGTCTTtacaaaaa                                                                                                                           >  4:24735/1‑147 (MQ=255)                                 tGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑TAAAAAAAATTCTTTACTTAAAAAAAAAATaa                                                                                                                      >  3:146219/1‑150 (MQ=255)                                  gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAAc                                                                                                                     <  2:951422/150‑1 (MQ=255)                                  gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACt                                                                                                                    <  1:539444/151‑1 (MQ=255)                                  gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CAAATAAAATTTTTGAATTaaaaaaaaaaaaa                                                                                                                      >  3:219731/1‑142 (MQ=255)                                   aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACtt                                                                                                                   <  1:1063856/151‑1 (MQ=255)                                     aaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGa                                                                                                                 <  1:342560/151‑1 (MQ=255)                                     aaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAATTTCTTGACTTTAAAAAAAAATaa                                                                                                                      >  3:270562/1‑146 (MQ=255)                                      aTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAGATTGCTTTACTTaaaaaaaaaaaaa                                                                                                                      >  1:1061635/1‑138 (MQ=255)                                        tgtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACa                                                                                                               <  1:1234769/150‑1 (MQ=255)                                         gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGa                                                                                                             <  1:716684/151‑1 (MQ=255)                                         gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAAATTCCTTAATTTAA‑AAACAAATAACTTgta                                                                                                                >  2:21446/1‑146 (MQ=255)                                           gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaa                                                                                                            <  1:658009/150‑1 (MQ=255)                                            cTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAATTTCTTGACTTTATAAACAAATAAATTgaaagaaa                                                                                                           >  2:982651/1‑150 (MQ=255)                                              aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGa                                                                                                        <  1:732021/151‑1 (MQ=255)                                               ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CAAAAAGATTTCTTTACTTTACAAAAAAATAACTTGACAGaaaa                                                                                                          >  4:75637/1‑148 (MQ=255)                                                   cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GTAAAAAATTGTTTTTCTTTACAAAAAAATAAATTTACAGaaaaaaa‑aaa                                                                                                   >  2:897642/1‑150 (MQ=255)                                                   cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAAATTCTTTACTTTACAAACAAATAACTTGACAGaaaaaaa‑aaa                                                                                                   >  3:162539/1‑150 (MQ=255)                                                      aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑aaaaaaa                                                                                               <  1:6175/151‑1 (MQ=255)                                                        aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAaga                                                                                             <  2:382171/151‑1 (MQ=255)                                                        aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAATTTCTTTACTTAAAAAAAAAAAAAACTTACAGaaaaaaa‑aaaaaaaa                                                                                              >  2:348831/1‑149 (MQ=255)                                                         aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAAATTCTTTACTTTAAAAAAAAAAAAATTTTCAGaaaaaaa‑aaaaaaaaa                                                                                             >  2:1039136/1‑148 (MQ=255)                                                             aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAg                                                                                        <  2:1150150/151‑1 (MQ=255)                                                                  tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CAAATAGATTGCTTTGCTTTACAAACAGAAAAATTTACAGAAAAGAa‑aaaaaaaaaaacaacggg                                                                                    >  1:490703/1‑139 (MQ=255)                                                                       aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTgg                                                                              <  2:670/151‑1 (MQ=255)                                                                          aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CAAAAAAATTTTTTGGCTTAACAAAAAAAAAAATTTACAGAAAAAAA‑AAAAAAAAAAACAGCggg                                                                                    >  1:740279/1‑142 (MQ=255)                                                                              aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                       <  2:84010/151‑1 (MQ=255)                                                                                gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAATTGCTTGACTTAACAAAAAAAAAACTTTACAGAAAGGAA‑AAAAAAAAAGACAGCGGGGGCCgg                                                                              >  1:1210279/1‑142 (MQ=255)                                                                                  ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                                   <  1:562425/151‑1 (MQ=255)                                                                                  ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAATTTTTTGACTTAACAAAAAAAAAAATTTACAGAAAAAAA‑AAAAAAAAAAACAGCGGGCGCGGGCTCtt                                                                         >  2:439197/1‑145 (MQ=255)                                                                                           tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTT‑TAAAAAGATTTCTTTTCTTAAAAAAAAAATAAATTTAAAGAAAAAAT‑AAAAAAAAAAACAACGGGGGCCGGGGCTtgatgaaactcaa                                                             >  1:214514/1‑141 (MQ=255)                                                                                           tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTT‑CAAAAAGAATTCTTTTCTTAATAAAAAAAAAAATCTACATAAAAAAA‑AAAAAAAAAAACAGCGGGGGCCGGCTCTTGTAAACGCCaaa                                                             >  1:575301/1‑148 (MQ=255)                                                                                           tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                               <  4:234047/146‑1 (MQ=255)                                                                                             tGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CAAAAAAAATTCTTGACTTTACAAACAAAAAACATGAAAGAAAAAAT‑AAAAAAAAAACCAACGGGGGCGGGCTCTTTAAGAAGCCaaa                                                             >  1:902257/1‑146 (MQ=255)                                                                                              gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                       <  1:327446/151‑1 (MQ=255)                                                                                               cgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CGAAAAGATTTCTTTACTTTACAAAAAAATAACATGACAGAAAAAAT‑AAAAAAAAAAACAACGGGGGCTGGGTCTTGTAGAAGCCAAAATaaaa                                                       >  4:436684/1‑150 (MQ=255)                                                                                               cgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑CCAATAGATTGCTTGTTTTTAAAAAAATAAAAATTGACAGAAAAAAT‑AAAAAAAAAGACAGCGGGGGCTGGCTCTTCAAAACGCGCAAATAGaa                                                       >  1:324818/1‑150 (MQ=255)                                                                                                  gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTg                                                   <  2:522954/151‑1 (MQ=255)                                                                                                  gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAAAATTTTGTGCTAAAAAAAAAAAAAATTGACGGGAAAGAT‑AAAAAAAAAAACACCCGGGGCTGGCCCTTTTATACACTaaaaaaaaag                                                      >  3:199063/1‑141 (MQ=255)                                                                                                       tGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAAATGCTTGACTTTACAAACAAAAAAATTGGCAGAAAAGAA‑AAAAAAAAAGCCCGCGGGCGCTGGCTCTTTGAGAAGCCAAATTAGAAGGTGGCggt                                              >  4:15188/1‑149 (MQ=255)                                                                                                               tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGgc                                                                                   <  2:1085150/106‑1 (MQ=255)                                                                                                               tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAATAAATTTCTTTGCTTTAAAAAAAGAAAAAATTACAGAAAAGATAAAAAAAAAAGCCA‑CGGGGGCCGGGTCTTTGAGAAGCCAAACTAAAAACTGGCGGTAc                                            >  1:1222172/1‑145 (MQ=255)                                                                                                                  aaTATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGt                                     <  1:1080748/150‑1 (MQ=255)                                                                                                                  aaTATTCCGTAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAAAATGTTTTTCTTAAAAAAAAAAAAAAATTACAGAAAAAAA‑AAAAAAAAAAACAGAGGGGGCCGGCTCTTTGAGACGCCAAAATAAAAACAGGCGGTAagcg                                         >  2:152987/1‑145 (MQ=255)                                                                                                                     atTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                <  2:383818/151‑1 (MQ=255)                                                                                                                       tCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                               <  2:531029/150‑1 (MQ=255)                                                                                                                          gTAACAGGATGATCGTTATGGCGGGGGTG‑GAAATAGATTGCTTGGCTTTACAAAAAAAAAAAATGACAGAAAAAAA‑AAAAAAAAAACCAGCGGGGGCCGGCTCTTTAAGAAGCCAAATTAAAAGATGGCGGTTCGCGGAggaagggggg                             >  4:363913/1‑140 (MQ=255)                                                                                                                           tAACAGGATGATCGTTATGGCGGGGGTG‑GAAAAAGAATGGTTTTCTTAACAAAAAAAAAAATTGGCCGGAAAAAA‑AAAAAAAAAAACAACGGGCGCCGGGCCTTGAAAAAGCTAAAATAAAAACTAGGGGTAcgccgggcgagggggga                            >  2:550949/1‑135 (MQ=255)                                                                                                                            aaCAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACg                          <  2:901457/150‑1 (MQ=255)                                                                                                                               aGGATGATCGTTATGGCGGGGGTG‑TGAATAAATTTCTTTTCTTTACAAACAGATAAATTGGCAGAAAAAAT‑AAAAAAAAAAACAGCGGGGGCTGGGTCTTTTAGAAGCTAAATTAAAAAATGGCGGTACGCGGAGGACGGGGGGACgg                         >  2:119596/1‑148 (MQ=255)                                                                                                                                       cGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatcac              <  2:277933/151‑1 (MQ=255)                                                                                                                                           aTGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAATGGAATCACATCacgtac          <  1:1037874/151‑1 (MQ=255)                                                                                                                                               cGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt      >  4:145835/1‑151 (MQ=255)                                                                                                                                                ggggTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTTCGTGGGAACGGAATCACATCACGTACGttt      <  1:785556/150‑1 (MQ=255)                                                                                                                                                ggggTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt     >  1:1136949/1‑151 (MQ=255)                                                                                                                                                   gTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa  <  1:243689/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTG‑CAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAT‑AAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCA  >  NC_000913/987347‑987643

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