Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	835,937	0	T	G	42.4%	106.1 / 11.8	61	G230G (GGT→GGG)	ybiB	putative family 3 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (12/22);  new base G (25/0);  total (39/22)
Fisher's exact test for biased strand distribution p-value = 6.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGATTGAGCG  >  NC_000913/835788‑836086                                                                                                                                                        |                                                                                                                                                      gCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGt                                                                                                                                                       <  2:722700/150‑1 (MQ=255)            ggCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATg                                                                                                                                             <  1:218989/149‑1 (MQ=255)              cGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACCAATGC‑ATGGT‑ACTGAAGGGGAAGGGTATggt                                                                                                                                           >  1:475238/1‑147 (MQ=255)              cGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTaa                                                                                                                                         <  2:361945/151‑1 (MQ=255)               gAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAATGC‑ATGGT‑GCTGAAGGGGGAGAGGATTCTaa                                                                                                                                         >  1:519575/1‑150 (MQ=255)                aaaCTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑CTGGT‑ACTGAAGGGGAAGGGTATGCTaa                                                                                                                                         >  2:334009/1‑149 (MQ=255)                 aaCTGGCGACTCCATTTGCCGAAATTGAGGCTCTCAGTCTTTCCCACCTTTCGCATCCGGAATACATTCTTCCCGTCGCGAAGTTCTTTAGCGATATCGTCAGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAACCTTAAGTGTATGCTAATcc                                                                                                                                      <  2:1103348/151‑1 (MQ=255)                      gCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGCGCCCACAATAA‑ATTGT‑AATTAAGGTGGAGGGGATTCTAATCCGCa                                                                                                                                   >  2:427217/1‑149 (MQ=255)                      gCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAgc                                                                                                                                 <  2:87623/151‑1 (MQ=255)                        gACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCTCCAATGC‑AAGGG‑ACTGGAGGGGAAGTGTGTGCTAATACGCCgcgc                                                                                                                               >  1:802227/1‑151 (MQ=255)                        gACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAgcgc                                                                                                                               <  1:185576/151‑1 (MQ=255)                            ccATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCCCCTAATGA‑AAGGTAATTGAAGG‑GAAGGGTATGCTAATCCCCAgcgccccc                                                                                                                           >  1:480841/1‑147 (MQ=255)                             cATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAATGC‑ATGGG‑ACTGGAGGGGGAGGGGATGCTAATTCGCAGCGCCGccc                                                                                                                          >  1:527581/1‑151 (MQ=255)                             cATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCAGGACTACATTGGACGCGTCGCGAAGTTCTTTAGAGATATCGGCGGGGGGGGGGCCCTCATCA‑AAAGT‑ACCGTAAGGGAAGTGACTGCCAGTCCGCAGCGCTGccc                                                                                                                          >  2:706758/1‑151 (MQ=255)                              aTTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAATGC‑ATGGT‑ACTGAAGGGGAAGAGTATGCTAATCCGCAGCGCTGccc                                                                                                                          >  2:408560/1‑150 (MQ=255)                                    cGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATc                                                                                                                   <  2:425594/151‑1 (MQ=255)                                          tGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACCAATGC‑ATGGT‑ACTGGAAGGGGAGTGTATGACAATTCGCCGCGCAGACCGCCGCTCaa                                                                                                                 >  1:86068/1‑147 (MQ=255)                                                gcTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAATGC‑ATGGT‑ACTGAAGGGGGAGGGTATGCTAATCCGCAGCGCCGCCCGCCGATCAATATCAtt                                                                                                          >  2:993667/1‑148 (MQ=255)                                                gcTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGGGGAGGGTATGCTAATCCCCAGCGCTGCCCGCCGATCAATCTCAtt                                                                                                          >  2:55222/1‑148 (MQ=255)                                                gcTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCCACTAAAAC‑ATGGG‑ACTGAAAGGGAAGGGGATGCCAATTCCCCGCGCTGCCCCCCGAACAATCTCAttt                                                                                                         >  2:950048/1‑148 (MQ=255)                                                  tGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACTAATGA‑ATGGG‑ACTTCAGGGGAAGTGTTTGCTAATTCGCAGCGCTGCCCGCCGGTCAAACTCAttt                                                                                                         >  1:199326/1‑146 (MQ=255)                                                   gCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCg                                                                                                     <  1:265263/150‑1 (MQ=255)                                                     gTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCCACTAAAGC‑ATGGT‑ACTGAAGGGGAAGTGTATGCTAATCCGCAGCGCTGCCCCCAGATCAAACTCAtt                                                                                                          >  2:585176/1‑143 (MQ=255)                                                            cAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCCACCAATGC‑AAGGG‑GATGCAGGGGAAGGGGGTGCTCATACCCCGCGCTGCCCGCCGATCAATCTCACTTGCCGTGGAgg                                                                                               >  1:24224/1‑147 (MQ=255)                                                              gCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCg                                                                                         <  2:151278/151‑1 (MQ=255)                                                              gCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCg                                                                                         <  2:968108/151‑1 (MQ=255)                                                               cGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGGTGCTAATGG‑ACGGG‑ATTTAAGGGGGAGTGGGGGTTAAACCGCCGGCGTTCCCCCAGAGCAAACTCAtttttccgggaggg                                                                                              >  1:93445/1‑133 (MQ=255)                                                                gTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGGTACTAATGC‑ATGGG‑ACTGAAGGGGAAAGGGATGATAAAACGCCGCCGCGCCCGCAGATCAATCTCATTGAGCGTGGAgggggggggg                                                                                       >  1:2456/1‑144 (MQ=255)                                                                   tCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTg                                                                                     <  1:624065/150‑1 (MQ=255)                                                                     gCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGGGCTACTAATTC‑ATGGT‑AATGAAGGTGGAGTGTATGCTAAACCCCAGCGCTGCCCCCCGATCAAACTCATTTAACGTGAAGGGGTGCggggg                                                                                     >  2:475336/1‑146 (MQ=255)                                                                       aTCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTa                                                                                <  1:307060/151‑1 (MQ=255)                                                                       aTCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCAGGCGGGACTCAAGC‑CTTGT‑ACTATAGGGGAAAAGAATGCTAATCCGCCGCGCTGCCCGCCGATCAACGTCATTTACCGGGGAGGGCGGCGGGt                                                                                      >  2:1092955/1‑145 (MQ=255)                                                                           ggAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGaa                                                                            <  1:534655/151‑1 (MQ=255)                                                                              aTACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGGGGGCGCCAATCATAC‑ATGGC‑ACTAAAAGGGGAGGGGTTTCTCAGTCGCAGCGCTCCCCCCCAAACCATTTCATTGGCCGCGTAAGGATGCGGGTGCGGGATGaaaa                                                                          >  2:1055207/1‑150 (MQ=255)                                                                                aCATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGGGACTAATGC‑ATGGT‑ACTGAAGGGGGAGAGGATGGTAATCCGCCGCGCCGCCCGCAGATCAATCTCCTTTATCGGGAAGGGGTGCGGGGGGTGGATGaaaaa                                                                         >  1:254093/1‑149 (MQ=255)                                                                                         gcgTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCCACTAATGC‑ATGGT‑ACTGGAGGGGGAGGGGATGCTAATCCGCAGCGCCGCCCCCAGATCAAACTCGTTGTCCGTGGAAGGAGGCGGGGGGTGGATGAAAAAAAGGAAAatgc                                                              >  2:516086/1‑151 (MQ=255)                                                                                           gTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGGGCCACTAATGC‑CTGGG‑ACTGAAGGGGAAGGGGATGCTAATCCGCAGCGCCGCCCCCAGAGCAATCTCATTGACCCCGAAGGGAGGCGGGGGGTGGATGAAAAACAGGAAAcagctg                                                            >  2:261483/1‑151 (MQ=255)                                                                                                gAAGTTCTTTAGCGATATCGGCGGGGGGGGGCCACTAATGC‑ATGGT‑ACTGGAGGGGAAGGGTATGCTAATCCGCAGCGCCGCCCGCAGAACAATCTCATTGACCGGGAAGGGGGGGGGGGGGTGGATGAAAAACAGGACAATGATGGGAg                                                        >  2:201909/1‑150 (MQ=255)                                                                                                 aaGTTCTTTAGCGATATCGGCGGGCGGGGGCTACTAAATC‑ATGGT‑ACTGAAGGGGAAGGGTATGCTAAACCCCAGCGCTGCCCGCAGAACAAACTCATTGACCGTGGAGGGGTGCGGGGGCTGTATGAAAAAAAGGACACTGCTGGTAGCg                                                      >  2:579046/1‑151 (MQ=255)                                                                                                  aGTTCTTTAGCGATATCGGCGGGCGGGGGGTACTTATGCGAGGGA‑ATGGAAGGGGAGGGGGGGGTTAA‑CCGCCGCGCCGGCCGCAGATCCAACTCATTTTCCGTGCAGGGGTGCGGGGGCTGTATGAAAAACAGGGCAccgct                                                             >  1:454053/1‑143 (MQ=255)                                                                                                          aGCGATATCGGCGGGCGGGCGCTTCTAATGC‑CTGGG‑AATGAAGGGGGAGGGGGTGCTAATCCGCAGCGCCGCCCGCCGACCAATCTCATTGAACGGGAAGGGGGGCGGGGGGGGGATGAAAAAAAAGACACTGCTGGGAGCggggttct                                               >  1:506071/1‑142 (MQ=255)                                                                                                                    gcgggcgggcgCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCaa                                    >  1:1043875/1‑150 (MQ=255)                                                                                                                      gggcgggcgCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCaaa                                   >  2:665097/1‑149 (MQ=255)                                                                                                                           ggcgCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCg                            >  1:440115/1‑151 (MQ=255)                                                                                                                             cgcTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCgg                           >  1:507282/1‑150 (MQ=255)                                                                                                                              gcTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGa                          >  2:841965/1‑150 (MQ=255)                                                                                                                               ctactaATGC‑ATGGT‑ACTGAAGGCGAAGTGTAGGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGaa                         >  2:382903/1‑150 (MQ=255)                                                                                                                                 actaATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGaaa                        <  1:261426/149‑1 (MQ=255)                                                                                                                                 actaATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAAcc                      <  1:472118/151‑1 (MQ=255)                                                                                                                                    aaTGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCa                     <  2:151044/149‑1 (MQ=255)                                                                                                                                    aaTGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACg                   <  2:784682/151‑1 (MQ=255)                                                                                                                                     aTGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACgg                  <  1:377711/151‑1 (MQ=255)                                                                                                                                      tGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGgc                 >  2:342357/1‑151 (MQ=255)                                                                                                                                        c‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGgcgc               >  2:180280/1‑151 (MQ=255)                                                                                                                                              t‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTgg          >  2:1083810/1‑151 (MQ=255)                                                                                                                                                   gAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGAtt       <  2:565924/150‑1 (MQ=255)                                                                                                                                                   gAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGAtt       <  1:595036/150‑1 (MQ=255)                                                                                                                                                   gAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGAt        <  1:8733/149‑1 (MQ=255)                                                                                                                                                   gAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGATTg      <  2:1051209/151‑1 (MQ=255)                                                                                                                                                     aGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGATTg      >  1:134142/1‑149 (MQ=255)                                                                                                                                                     aGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGATTGAg    <  2:527467/151‑1 (MQ=255)                                                                                                                                                       gTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGATTGAGCg  >  1:843729/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GCACATACCCTGGCGAAACTGGCGACGCCATTTGCCGAAGGTGAGGCGCTGCGTCTTTCCAGCGTTTCGCATCCGGAATACATTGGACGCGTCGCGAAGTTCTTTAGCGATATCGGCGGGCGGGCGCTACTAATGC‑ATGGT‑ACTGAAGGTGAAGTGTATGCTAATCCGCAGCGCTGCCCGCAGATCAATCTCATTGACCGTGAAGGGATGCGGGTGCTGTATGAAAAACAGGACACTGCTGGTAGCGAGTTACTGCCACAAGCAAAAGATCCGGAAACCACGGCGCAGTGGATTGAGCG  >  NC_000913/835788‑836086

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