Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	753,125	0	T	A	42.6%	84.6 / 33.6	61	intergenic (‑330/+60)	ybgD/gltA	putative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (29/6);  new base A (0/26);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCA  >  NC_000913/752976‑753266                                                                                                                                                       |                                                                                                                                              cATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCatat                                                                                                                                               >  1:94458/1‑150 (MQ=255) cATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCatat                                                                                                                                               >  1:817020/1‑150 (MQ=255)  aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCAtata                                                                                                                                              >  2:89515/1‑150 (MQ=255)     ttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATaaaa                                                                                                                                           <  1:253183/150‑1 (MQ=255)     ttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATaaaa                                                                                                                                           >  1:765126/1‑150 (MQ=255)       cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATc                                                                                                                                        >  1:774635/1‑151 (MQ=255)           gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAAcc                                                                                                                                    <  2:31368/151‑1 (MQ=255)            ggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAAccc                                                                                                                                   <  2:20873/151‑1 (MQ=255)                ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGcc                                                                                                                                >  2:811239/1‑150 (MQ=255)                 taAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCat                                                                                                                              >  1:298005/1‑151 (MQ=255)                            aCCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcgg                                                                                                                   >  2:267409/1‑151 (MQ=255)                                 gAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGtt                                                                                                                >  2:399950/1‑149 (MQ=255)                                 gAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTa                                                                                                               >  1:592093/1‑150 (MQ=255)                                       tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATAtt                                                                                                        >  2:295480/1‑151 (MQ=255)                                       tttAATTAAAAGTTTTGGGTTGCTAATTTTTTCGATTTTTATATTTTTTTTAAATACCTTAATGTTTTTGATATTGATTTTTTTTTTAACAATTTTGGG‑TTTTAAAAAATAAAAAAGCCAACCCGCCATATGAACGGCGGGTTAAAata                                                                                                          <  2:573829/149‑1 (MQ=255)                                       tttAATTAAAAGTATGGGGTTGTTGATATTTTGGGTTGTTTTTTTTTTTTTAAATATTAATATAGGTTTGTTAAGGGATTTTTTTTTTAAAATTTTTGT‑TTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAatat                                                                                                         <  2:886040/150‑1 (MQ=255)                                          aaTTGAAAGTATTGGGTTGCTTATAATTTGAGCTTTTCTATTTTTTTAAAATATCTATATAGGTTTTTAAATGGTTTTTTTTTTTACAATTTTTTG‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                       <  1:779714/149‑1 (MQ=255)                                          aaTTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                       >  2:684201/1‑149 (MQ=255)                                                 aGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                                >  1:790075/1‑149 (MQ=255)                                                      tGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCaa                                                                                           >  1:361560/1‑149 (MQ=255)                                                       gggTTGCTAATAATTTGAGTTGTTCTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTATTTTTACATTTTTTTG‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATc                                                                                         <  1:227968/150‑1 (MQ=255)                                                             cgctcttccgattttGTTTTTTTTTTTTAAAATATTTATATAGGTTTGTTAATGGATTTTATTTTTAAAATTTTTTG‑TGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                                   <  2:828171/136‑1 (MQ=255)                                                               gATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                                 >  1:459065/1‑150 (MQ=255)                                                                  aaTTTGAGCTTTTCTTTTCTTTTTAAATATCTATATAGGTCTGTTAAGGAATTTTATTTTTACATTTTTTTG‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                              <  1:84738/150‑1 (MQ=255)                                                                      tgattttttCTTTTTTTTTAAAATATTTATATAGGTTTGTTAATGGATTTTATTTTTAAATTTTTTTG‑TGTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTg                                                                          <  1:128857/143‑1 (MQ=255)                                                                       gAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGa                                                                         >  2:359305/1‑150 (MQ=255)                                                                           tttttttttttttttAAATTTTTATATAGGTTTTTAAAGGGTTTTTATTTTTAAAATTTTTTT‑TTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                     <  2:193953/140‑1 (MQ=255)                                                                           tGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAcc                                                                                     >  1:828238/1‑134 (MQ=255)                                                                            gTTCTTTTCTTTTTAAATATCTATATAGGTCTGTTAATGAATTTTATTTTTACATTTTTTTG‑TGTTTAGGCATATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCg                                                                    <  1:584429/150‑1 (MQ=255)                                                                             ttttttttttttaAAATATCTATATAGGTCTTTTAAGGTATTTTTTTTTTACATTTTTTTG‑TGTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                     <  1:658645/140‑1 (MQ=255)                                                                             ttCTATTTTTTTAAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTACAATTTTTTG‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                                     <  1:274352/148‑1 (MQ=255)                                                                                tATTTTTTTTAAATATCTATATAGGTTTGTTAATGGATTTTTTTTTTAAAATTTTTTG‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtt                                                                 <  1:404241/149‑1 (MQ=255)                                                                                tATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTATAATTTTTATAACTTAGCAATCAACCATTAACGCTTGCTATCGCtttt                                                               >  2:746437/1‑151 (MQ=255)                                                                                 aTTCTTTTTAAATATCTATATAGGTCTGTTAATGGCTTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAATAATTTTTTAACTTTGCAATCAACCCATAATGCTTGATATCGCcttt                                                               >  2:14380/1‑150 (MQ=255)                                                                                   tttttttAAATATCTAAATAGGTTGGTTAATGGTTTTTATTTTTAAAATTTTTTG‑TTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaa                                                             <  1:253149/148‑1 (MQ=255)                                                                                        ttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcg                                                        >  1:335786/1‑150 (MQ=255)                                                                                        ttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgc                                                       >  1:280117/1‑151 (MQ=255)                                                                                                 tatatgGGTCTGTTAATGGTTTTTTTTTTTACATTTTTTTG‑TGTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                               <  1:839599/150‑1 (MQ=255)                                                                                                 tatataGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAAATTAACAAACAACCATTAAAGCGTGGTATTCCCTTTTAAGTACCGttttt                                                 >  2:137095/1‑148 (MQ=255)                                                                                                      aGGTCGTTTAAGGAATTTTTTTTTTACATTTTTTTG‑TGTTGAGGAATAAAAAACTCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                                <  1:55124/144‑1 (MQ=255)                                                                                                       ggTCTTTTAATGGTTTTTTTTTTTAAAATTTTTTG‑TGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                              <  1:780288/145‑1 (MQ=255)                                                                                                         tttttaaaTGGTTTTTTTTTTTAAATTTTTTTT‑TGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                       <  1:810865/144‑1 (MQ=255)                                                                                                         tCTGTAAATGGTTTTTATTTTTAAAATTTTTTT‑TGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCt                                         <  2:863307/148‑1 (MQ=255)                                                                                                              tAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATTTTTACAACTTAACAAACCAACATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGTATACa                                  >  2:605092/1‑150 (MQ=255)                                                                                                               aaTGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGTATACAGc                                >  2:111412/1‑151 (MQ=255)                                                                                                                      tttattttTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATACTGTATACAGCTGAcg                           >  1:821068/1‑149 (MQ=255)                                                                                                                       ttatttatACATTTTTTTG‑TTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg                         <  1:679767/150‑1 (MQ=255)                                                                                                                        tattttTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                        >  1:704603/1‑150 (MQ=255)                                                                                                                          tttttAAAATTTTTTG‑TTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                      <  2:275610/150‑1 (MQ=255)                                                                                                                          tttttA‑AATTTTTTGTTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                      <  1:612506/150‑1 (MQ=255)                                                                                                                          ttttaAAAATTTTTTT‑TTTTTAGTAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                        <  2:372363/148‑1 (MQ=255)                                                                                                                           ttttACAATTTTTTG‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACg                     <  1:418196/150‑1 (MQ=255)                                                                                                                               aCATTTTTTTT‑TTTTAAGCCATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                  <  2:875373/149‑1 (MQ=255)                                                                                                                                cAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAAATTTTATAAATTTGCTAACCACCAATAACTCCTGATAGTGCTTTTAAAGTCGCGTTTTTTATATACCCTATACAGTCGACGTGGagaggcca                >  2:361570/1‑143 (MQ=255)                                                                                                                                 atttttttG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa                 <  1:61575/146‑1 (MQ=255)                                                                                                                                  aTTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCCTTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGGCGGGGAATCt             >  1:699098/1‑151 (MQ=255)                                                                                                                                   ttttttG‑TTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCaa                <  1:791647/147‑1 (MQ=255)                                                                                                                                    tttttt‑ttttTAGAAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt            <  1:184753/150‑1 (MQ=255)                                                                                                                                    tttttG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTTGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAAGCGGACGGGCAATCtt            >  1:820796/1‑150 (MQ=255)                                                                                                                                     ttttG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTc           >  1:836808/1‑150 (MQ=255)                                                                                                                                      tttttttttaagcaatatAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTc           <  1:337971/136‑1 (MQ=255)                                                                                                                                       ttG‑TTTTAAGGCATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCAt         <  2:527755/150‑1 (MQ=255)                                                                                                                                        tg‑tgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCAt         >  1:527805/1‑149 (MQ=255)                                                                                                                                              ttAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTc   >  2:367492/1‑150 (MQ=255)                                                                                                                                                aGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCa  <  1:307605/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                              CATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCA  >  NC_000913/752976‑753266

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