Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,503	0	G	A	43.2%	87.4 / 43.6	82	intergenic (‑521/+82)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (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 = 1.00e+00

ATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAC  >  NC_000913/987356‑987653                                                                                                                                                       |                                                                                                                                                        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                                                                                                                                        <  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: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:805486/1‑151 (MQ=255)                                                                                                                                        ggggTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt               >  2:1706332/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:62500/1‑151 (MQ=255)                                                                                                                                               gC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc         >  1:1874449/1‑151 (MQ=255)                                                                                                                                                       gATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc  <  2:268617/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        ATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTT‑GC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAC  >  NC_000913/987356‑987653

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