Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,494	0	T	G	52.9%	114.4 / 148.1	85	intergenic (‑512/+91)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (6/34);  new base G (45/0);  total (51/34)
Fisher's exact test for biased strand distribution p-value = 6.24e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-16

TTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAC  >  NC_000913/987349‑987644                                                                                                                                                   |                                                                                                                                                          ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑g                                                                                                                                                        <  2:1788679/148‑1 (MQ=255)   cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAAt                                                                                                                                                  <  1:201051/151‑1 (MQ=255)   cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAAt                                                                                                                                                  <  1:2001451/151‑1 (MQ=255)        aTCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGAtt                                                                                                                                             <  2:276445/151‑1 (MQ=255)         tCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGttttcaa                                                                                                                                                     >  1:99679/1‑139 (MQ=255)           cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAATTTCt                                                                                                                                          >  2:356092/1‑151 (MQ=255)                tGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GT‑aaaaaaaatttttgtct                                                                                                                                     >  2:1821645/1‑137 (MQ=255)                tGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑aaaaaaaatgtttttat                                                                                                                                     >  2:1107985/1‑137 (MQ=255)                tGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAATTGCtttaat                                                                                                                                     >  1:1164593/1‑147 (MQ=255)                 gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑tt                                                                                                                                   <  2:59499/151‑1 (MQ=255)                 gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑t                                                                                                                                    <  2:1374078/150‑1 (MQ=255)                  ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAGATTTATTGTC‑tt                                                                                                                                   >  2:1351413/1‑150 (MQ=255)                   cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑aaaaaaaattctt                                                                                                                                         >  1:1389646/1‑134 (MQ=255)                     aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTaaca                                                                                                                               <  2:1389622/151‑1 (MQ=255)                     aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTaaca                                                                                                                               <  1:517587/151‑1 (MQ=255)                        aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAATAAAATTCTTTAC‑TTtaaaaa                                                                                                                             >  1:184016/1‑150 (MQ=255)                          tAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAAATTTTTTTTC‑TTaaaaaaaa                                                                                                                           >  2:141055/1‑148 (MQ=255)                           aGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAg                                                                                                                          <  1:1670306/150‑1 (MQ=255)                              aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATa                                                                                                                       <  2:1715526/150‑1 (MQ=255)                              aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATa                                                                                                                       <  1:613356/150‑1 (MQ=255)                              aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATa                                                                                                                       <  1:962060/150‑1 (MQ=255)                              aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAAAATGCTTTAC‑TTTACAAAAAAATa                                                                                                                       >  1:807725/1‑150 (MQ=255)                                  aaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTg                                                                                                                  <  1:2151460/151‑1 (MQ=255)                                   aaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGa                                                                                                                 <  1:127669/151‑1 (MQ=255)                                      tgtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAAATTTTTTAC‑TTaaaaaaaaaaaaaatttaaa                                                                                                               >  2:691231/1‑136 (MQ=255)                                         gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GT‑AAAAAAATTTCTTTAC‑TTTAAAAAAAAAAAACTTTACAGaa                                                                                                            >  1:68754/1‑150 (MQ=255)                                          cTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑GAAAAAATTGTTTTAC‑TTTACAAAAAAAAAAATTTACAGaaa                                                                                                           >  2:1498230/1‑150 (MQ=255)                                            aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAg                                                                                                         <  2:1302178/150‑1 (MQ=255)                                             ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GG‑AAAAAAATTGTTTGAC‑TTAAAAAAAAAAAAAAATTACAGaaaaaa                                                                                                        >  2:542575/1‑148 (MQ=255)                                             ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑GAAAAAATTTCTTGAccttataaa                                                                                                                              >  1:2108698/1‑121 (MQ=255)                                               aGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATa                                                                                                      <  1:497525/150‑1 (MQ=255)                                                gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATaaa                                                                                                    <  1:1080572/151‑1 (MQ=255)                                                 cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTT‑TC‑AAA‑ATATTGTTTTGA‑TTTACAAAAAGAAAAATTGCGAGAAAAAATaaa                                                                                                    >  1:1726866/1‑150 (MQ=255)                                                 cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAATAAATTTTTTGAC‑TTAAAAAAAAAATAAATTGAAAGAAAAAATaaa                                                                                                    >  2:283464/1‑150 (MQ=255)                                                 cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAAATTGTTGAC‑TTTAAAAAAAAAAAAACTGACAGaaaaaaaaaa                                                                                                    >  1:1531839/1‑150 (MQ=255)                                                   aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATaaaaaa                                                                                                 <  2:1909979/151‑1 (MQ=255)                                                    aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑TG‑AAATAGATTGCTTGAC‑CTAAAAAAAAAAAAAATTCACAGCAAAAATaaaaaa                                                                                                 >  1:1327102/1‑150 (MQ=255)                                                    aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAGATTTGGTTAA‑TTAAAAAACAAAAAAATTGACAGAaaaaaaaaaaaa                                                                                                 >  2:347934/1‑150 (MQ=255)                                                     gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAag                                                                                               <  2:864862/151‑1 (MQ=255)                                                     gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAATTGGTTGAC‑TTTACAAAAAAATAAAATTAAAAAAAAGAAtaaaaaa                                                                                                >  1:1772362/1‑150 (MQ=255)                                                      aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAaga                                                                                              <  1:1754970/151‑1 (MQ=255)                                                        aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAATAAAATGCTTTAC‑TTTAAAAAAAAAAAAATTTACAGaaaaaaaaaaaaaaaa                                                                                              >  2:895830/1‑147 (MQ=255)                                                         aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAAATGGTTGTT‑TTTACAAAAAAAAAAAATTACAGaaaaaaaaaaaaaaaa                                                                                              >  2:1696335/1‑146 (MQ=255)                                                         aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑GAATAGATTGCTTTGC‑TTAACAAACAAAAAAATTTACAGaaaaaaaaaaaaaaaa                                                                                              >  1:943278/1‑146 (MQ=255)                                                          aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAGATTGTTTTTC‑TTAAAAAAAAGATAACTTGACAGAAAAGATaaaaaaaaa                                                                                              >  1:1663271/1‑145 (MQ=255)                                                             tAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAGAATTCTTTAC‑TTAACAAACAAAAAAATTTACAGaaaaaaaaaaaaaaaa                                                                                              >  1:606963/1‑142 (MQ=255)                                                              aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGc                                                                                        <  1:1242943/149‑1 (MQ=255)                                                              aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGc                                                                                        <  1:270847/149‑1 (MQ=255)                                                              aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                                      <  2:833465/151‑1 (MQ=255)                                                              aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                                      <  1:295525/151‑1 (MQ=255)                                                               aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAAATTGCTTTAC‑TTTACAAAAAAAAAAATTTACAGAAAAAATAAAAAAAAAGCCCGCggg                                                                                     >  1:418192/1‑151 (MQ=255)                                                                 ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                    <  2:619138/150‑1 (MQ=255)                                                                 ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAATAAATTGGTTTAC‑TTAACAAAAAAAAAAATTTACAGAAAAAATaaaaaaaaaaaccacggg                                                                                     >  2:2171670/1‑138 (MQ=255)                                                                  cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAAATGCTTGGC‑TTAACAAACAGAAAAATTTACAGAAAAAAAAAAAAAAAAGACAACggg                                                                                     >  1:452026/1‑148 (MQ=255)                                                                   gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAAATTGTTTAA‑TTAAAAAAAAAAAAAAATGACAGAAAAAAAAAAAAAAAAAACAGCggggg                                                                                   >  1:930334/1‑147 (MQ=255)                                                                    cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAAATTATTGTC‑TTAAAAAAAAAAAAACATTACAGAAAAAAAAAAAAAAAAAACAGCGGgcgc                                                                                  >  2:569127/1‑149 (MQ=255)                                                                     aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAAATTGCTTGTC‑TTTAAAAAAAAAAAAATTTAAAGAAAAAATAAAAAAAAAACCAGGGGGGGCCgg                                                                               >  1:933998/1‑151 (MQ=255)                                                                         aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctct                                                                           <  2:428940/151‑1 (MQ=255)                                                                          gAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GT‑AAAAAAAATGTTTTAC‑TTTAAAAAAAGAAAAATTGACAGAAAAAATAAAAAAAAAAACAGCGGGCGCGGGctcct                                                                          >  2:2017948/1‑149 (MQ=255)                                                                             tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGG‑GG‑AAAAAGAATTATTTTG‑TTAAAAAAAAAATAACTTTACAGAAAAGAAAAAAAAAAAAAAAACGGGCGCTgg                                                                               >  1:1011359/1‑143 (MQ=255)                                                                              gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAg                                                                      <  1:550614/151‑1 (MQ=255)                                                                                ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAAATTGTTTAC‑TTAAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAGACAGCGGGGGCTGGCTCTTTAAg                                                                      >  1:1043541/1‑149 (MQ=255)                                                                                  aTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAATTTCTTTGT‑TTTACAAACAAAAAAATTTAAAGAAAAAAAAAAAAAAAAAACAACGGGGGCTGGCCCtttta                                                                       >  2:1657465/1‑143 (MQ=255)                                                                                      gTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑TT‑AAAAAAATTGCTTGAT‑TTAATAAAAAGATAAATTTACAAAAAAAAAAAAAAAAAAAACCACGGGCGCTGGCCCTTTAAGACGCTaa                                                               >  1:1354989/1‑150 (MQ=255)                                                                                       tttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGGG‑AAAAAAATTG‑TTGAA‑TTAAAAAAAAGAAAAATTGAAAGAAAAAAAAAAAAAAAAAAAAACGGGGGGTGGGTCTTGAAAACGCCaa                                                               >  1:1581104/1‑149 (MQ=255)                                                                                           tGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGG‑GG‑AAAAAAAATTTTTTAC‑TTAAAAAAAAAAAAAAATTACAGAAAAAAAAAAAAAAAAAACAACGGGCGCTGGGCCTTGAAAACGCCaaa                                                              >  2:258861/1‑146 (MQ=255)                                                                                              gcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                                      <  1:1789706/151‑1 (MQ=255)                                                                                                 tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAATTGTTTTTT‑TTTAAAAAAAAAAAAATTGTCAGAAAAAAAAAAAAAAAAACCAGCGGGCGCCGGCTCTTTAAGACGCCAAATTAGaaagtgg                                                   >  1:256473/1‑146 (MQ=255)                                                                                                 tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAAAATTCTTTAC‑TTAAAAAAAAAATAAATTTAAAAAAAAAATAAAAAAAAAAACAGCCGGGGCCGGCTCTTTAAGACGCCAAATTaaaaaatgg                                                   >  1:1699110/1‑146 (MQ=255)                                                                                                   aaTGTTGCAATGTAATATTCTGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAATAGAATTTTTTAC‑TTTAAAAAAAGAAAAATTTAAAGAAAAAAAAAAAAAAAAAGCAACGGGGGCTGGGTCTTGAAGAAGCtaaaataaaa                                                        >  2:994103/1‑144 (MQ=255)                                                                                                       ttGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAAATGTTTGAT‑TTAAAAAAAAAATAAATTTAAAGAAAAAAAAAAAAAAAAACCAACGGGGGCCGGGGCTTTAAGAAGCCAAAATaaaaaaagggg                                                 >  2:1959281/1‑140 (MQ=255)                                                                                                        tGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAAAATGGTTGAC‑TTAAAAAAAAAAAAAAATGACAAAAAAAAAAAAAAAAAAAACAACGGGGGCTGGCTCTTTAAGAAGCCAAAATaaaaagaggcggt                                               >  1:988127/1‑139 (MQ=255)                                                                                                          cAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GCAAAAAAAATTGTTTTA‑‑TTAAAAAAAAAAAAAAATTAAAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGCTCTTGAAGACGCTAAAATAGAAGATGGCGGTAcg                                            >  2:1522466/1‑149 (MQ=255)                                                                                                              gTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGt                                      <  2:1399618/151‑1 (MQ=255)                                                                                                                 atatTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                   <  1:1051382/151‑1 (MQ=255)                                                                                                                 atatTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑GAATAAAATTTTTGTCTTTAACAAA‑AGAAAACTTTACAGAAAAAAAAAAAAAAAAAACAGCCGGGGCTGGCCCTTTTAGAAGCTAAATTaaaaacaaggggtacgcg                                          >  2:771139/1‑130 (MQ=255)                                                                                                                 atatTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑CAAAAAATTGCTTGAC‑TTAAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAACCAGCGGGGGCCGGCCCTTTAAGACGCCAAAATAAAAACTGGCGGTTCGCCGAGTa                                     >  1:1110896/1‑149 (MQ=255)                                                                                                                   atTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GC‑AAAAAGAATGCCTTAC‑TTAAAAAAAAGAAAAATTTAAAAAAAAAAAAAAAAAAAAAACAACGGGCGGTGGCTCTTTAAAAAGCCAAATTAAAAAATGGCGTTTCGCGGAGGACggg                                 >  1:1739277/1‑149 (MQ=255)                                                                                                                      ccGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                               <  2:1600219/150‑1 (MQ=255)                                                                                                                          aaCAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAATAAATTTCTTTAC‑TTTACAAAAAAAAAACTTTAAAGAAAAAAAAAAAAAAAAGCCAGCGGGGGCCGGCTCTTGTAGACGCTAAATTAGAAACTGGCGGTTCGCGGAGGACGTgggga                             >  1:952639/1‑146 (MQ=255)                                                                                                                            cAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                        <  1:1401537/151‑1 (MQ=255)                                                                                                                                  gatCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                  <  2:922939/151‑1 (MQ=255)                                                                                                                                   atCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatc                 <  2:1122754/151‑1 (MQ=255)                                                                                                                                               ggggTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt      >  2:1706332/1‑151 (MQ=255)                                                                                                                                               ggggTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt      >  2:805486/1‑151 (MQ=255)                                                                                                                                               ggggTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTAGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt      >  2:1200840/1‑151 (MQ=255)                                                                                                                                                   tt‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc  >  1:971173/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                          TTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAC  >  NC_000913/987349‑987644

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