Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,501	0	T	A	43.2%	109.9 / 42.6	82	intergenic (‑519/+84)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (15/31);  new base A (35/0);  total (51/31)
Fisher's exact test for biased strand distribution p-value = 2.66e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGC  >  NC_000913/987351‑987646                                                                                                                                                          |                                                                                                                                                   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)         cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGT‑GG‑AAAAAAATTTCt                                                                                                                                            >  2:356092/1‑151 (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)                   aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTaaca                                                                                                                                 <  1:517587/151‑1 (MQ=255)                   aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTaaca                                                                                                                                 <  2:1389622/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:962060/150‑1 (MQ=255)                            aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATa                                                                                                                         <  1:613356/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:270847/149‑1 (MQ=255)                                                            aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGc                                                                                          <  1:1242943/149‑1 (MQ=255)                                                            aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                                        <  1:295525/151‑1 (MQ=255)                                                            aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                                        <  2:833465/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)                                                                                                                                                    gC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc  >  1:1874449/1‑151 (MQ=255)                                                                                                                                                    gC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc  >  1:62500/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                   CAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGC  >  NC_000913/987351‑987646

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