Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,845,823	0	G	T	44.3%	93.5 / 27.3	61	R449S (CGC→AGC)	uhpT	hexose phosphate transporter
Reads supporting (aligned to +/- strand):  ref base G (27/7);  new base T (0/27);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 3.53e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTG  >  NC_000913/3845674‑3845968                                                                                                                                                          |                                                                                                                                                   cAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGt                                                                                                                                                   <  2:408078/151‑1 (MQ=255)  aGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGtt                                                                                                                                                  <  2:191995/151‑1 (MQ=255)     ggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGttctt                                                                                                                                               >  2:724472/1‑151 (MQ=255)         aTTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCAt                                                                                                                                           >  1:86195/1‑151 (MQ=255)            tGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTa                                                                                                                                         <  1:935076/150‑1 (MQ=255)                    agaAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAGTTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACt                                                                                                                                >  1:195503/1‑150 (MQ=255)                     gaAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTAt                                                                                                                              >  1:928233/1‑151 (MQ=255)                        gTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGc                                                                                                                           >  1:400320/1‑151 (MQ=255)                          ttGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGcc                                                                                                                          >  2:317049/1‑150 (MQ=255)                           tGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCACCACTATCGCCAt                                                                                                                        >  1:563901/1‑151 (MQ=255)                                    gacgtgggAAGTCAGGAAAAAAAAAAAAAGGCCGGGCGTCATGGGGGGAAAAAGTCACCATTTA‑GTTTTTGCCCCTTTTCAACGGCT‑AATTTTTTT‑CCCGGCGGCGTTTTTT‑GCTTTTTCCAATTA‑CCGCCACTATCGCCATCAGGCa                                                                                                                  <  1:878370/142‑1 (MQ=255)                                     aCGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGAt                                                                                                               >  2:561463/1‑150 (MQ=255)                                      cGCTGGAAGCCACGAAAAAAAAAAAAAACCCCGGGGTCGGGCGGGAAAAAAGCCCCCAGTTACTTTTAGCCCCC‑TGCCACCTGCTAATTTTTTTC‑CCC‑GGGGGGTTTTTT‑GTTTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATa                                                                                                              <  1:494499/150‑1 (MQ=255)                                            aggTCAAGAAAAAAAAAAAAACCCCCGGCGTCAGGCGGGGAAAAATCCACCAGTTCGTTTTATCCCAC‑TACCACCTGC‑GAATTTTTTTCCCC‑GCGGCGGTTTTT‑GCTTTCTCCCATTA‑CCGCCACTATCGCCATCAGGCAGATAc                                                                                                             <  2:659561/143‑1 (MQ=255)                                            aaGTCACGAAAAAAAAAAAAAACCCCGGCCTCATCCCGGTAAAAACTCCCCAGTTACTTTTGTGCCCC‑TTCCAGTCGCTAATTTTTTTT‑CCC‑GCGGGGGTTTTT‑GGTTTTTTCCATTC‑CCGCCACTATCGCCATCAGGCAGATACaa                                                                                                           <  2:473712/147‑1 (MQ=255)                                                    gAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggcg                                                                                               >  1:170056/1‑151 (MQ=255)                                                     aaaaaaaaaaaacGCCCGGCGTCAGGCCGGGAAAAAGCCCCCAGTTACGTTTATCCCCC‑TTTCAACTGCTAATTTTTTTT‑CCC‑GCGGGGTATTTG‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggc                                                                                                <  2:722367/143‑1 (MQ=255)                                                     aaaaaaaaaaaaaCCCGGGCGTCTTGCGGGAAAAATGTCCCAAGTTCTTTTTATCCCCC‑TGCCACCTCGTAATTTTTTTT‑CTC‑GGGGTGATTTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGa                                                                                              <  2:216382/145‑1 (MQ=255)                                                     aaaaaaaaaaaaaCCCCGGCGTACTGGCGGGCAAAAAGCCCCCGTTATTTTTTTTCCCC‑CTGCAAGCGCCTAATTTTTTTCCGC‑GCGTCGATTTTG‑GCTTTCTCCCATTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggcg                                                                                               <  1:303821/145‑1 (MQ=255)                                                      ataataAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGat                                                                                             >  2:315543/1‑151 (MQ=255)                                                         ataaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATc                                                                                          >  2:414906/1‑151 (MQ=255)                                                          aaaaaaaaCCCCGGCTTCAGGCGGAAAAAAAGCCCCCGTTTACGTTTACCCCAC‑TTTCAACTGTTAATTTTTTTC‑CCC‑GCGGCGGTTTTT‑GCTTTCTCCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                                                                                         <  1:210985/150‑1 (MQ=255)                                                              aaaaCCCGGGCATCGTGCCGGGAAAAAGCCCCCTGTTATGTTTACCCCCC‑TGAAAGCTGTTTATTTTTTTC‑CCC‑GGGGCGATTTTT‑GCTTTCTCCCTTTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAgcg                                                                                     <  1:745770/151‑1 (MQ=255)                                                                aaccccGGCGTCTGGCGGGAAAAAAGCCACCTTTTGCTTTTATCCCAC‑TGTCAACTGCTAATTTTTTTC‑CGG‑GGGGTGGTTTTT‑TCGTTTTCCCTTTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGcgcg                                                                                   <  2:719752/148‑1 (MQ=255)                                                                   cccGGCGTCAGGCGGGGAAAAGGACCCCATTTACTTTTATCCCAC‑TGAAAACTACTAATTTTTTTT‑CCC‑GCGGCGTTTTTT‑GCTTTCTCCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcggcg                                                                                <  2:387041/151‑1 (MQ=255)                                                                        cGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAgg                                                                            >  2:465862/1‑150 (MQ=255)                                                                         gTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCCATCGCGGCGATATCAAACGCGGCGAAGGtg                                                                          >  1:555846/1‑151 (MQ=255)                                                                          tCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGtgc                                                                         >  1:933303/1‑151 (MQ=255)                                                                          tCACGCCGGGAAAAACCCCCCATTTCTTTTTATCCCTC‑TGCCACCTG‑TTAATTTTTTT‑CCCGGCGGCGGTTTTT‑GCTTTCTTCCTTTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGt                                                                           <  1:875534/149‑1 (MQ=255)                                                                                cGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGccc                                                                   >  2:632230/1‑151 (MQ=255)                                                                                   gCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGc                                                                >  2:19384/1‑151 (MQ=255)                                                                                           cACCAGTTACGTTTATGCCAC‑TGTCAACCGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAgg                                                        >  2:506830/1‑151 (MQ=255)                                                                                               aGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCg                                                     >  1:158136/1‑150 (MQ=255)                                                                                                   atttttATCCCAC‑TGACACCTGTTATTTTTTTTT‑CTC‑GGGGGGTTTTTT‑GCTTTCTCCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAAt                                                  <  2:823114/146‑1 (MQ=255)                                                                                                      tttagcccAC‑TGTCACCTGCTAATTTTTTTT‑CCC‑G‑GGGGGTTTTT‑GGTTTTCTCCCTTACCCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATAcc                                               <  1:535785/143‑1 (MQ=255)                                                                                                      tttATCCCTC‑TGTCATCTGTTAATTTTTTTT‑CCC‑GCGGGGTTTTTT‑GGTTTCTTCCTTTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATAc                                                <  2:559895/148‑1 (MQ=255)                                                                                                        tATCCCTC‑TGTCAACTGCTGTTTTTTTTC‑CTC‑GCGGCGATTTTT‑CCTTTCTTCCATTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGc                                            <  2:367884/150‑1 (MQ=255)                                                                                                            cccc‑TGTCACCTTCTAATTTTTTTT‑CCC‑GCGGGGGTTTTT‑GCTTTCTCCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTc                                         <  2:168505/146‑1 (MQ=255)                                                                                                            ccAC‑TGTCAACTGCTGATTTTTTTT‑CGG‑GGGGGGTTTTTT‑GCTTTTTTCCTTTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGc                                            <  2:895515/146‑1 (MQ=255)                                                                                                                  gTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGCTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCg                                   >  1:43088/1‑150 (MQ=255)                                                                                                                     aaCTGTTAATTTTTTTC‑CCC‑GGGGGGTTTTTT‑GCTTTCTCCCTTTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCaa                               <  2:901434/151‑1 (MQ=255)                                                                                                                       acggagtaTTTTTTT‑CCC‑CGGGGGGGTTTT‑TGTTTCCTCCTTTACCCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCa                                <  1:785115/142‑1 (MQ=255)                                                                                                                        tgtttatttttttt‑CCC‑GGGGGGTTTTTT‑GCTTTCTCCCATTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAAt                              <  2:523207/142‑1 (MQ=255)                                                                                                                        tgctgAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATc                             >  1:254089/1‑150 (MQ=255)                                                                                                                         gctgAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCAt                           >  1:757796/1‑151 (MQ=255)                                                                                                                           acttttttttt‑CCC‑GGGGGGGTTTTTCGTGTTCTTCCCTTC‑CCGCACCTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCa                            <  2:330546/145‑1 (MQ=255)                                                                                                                              aTTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCATCGCGGCGAAGGAGCCTGCACAGTCGGTAAGGCCGAATATCGGCGTCCCACCGGCAATCATTccc                       >  1:722359/1‑150 (MQ=255)                                                                                                                               ttttttt‑CCC‑GGGGGGGTTTTT‑GCTTTCTCCCATTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTcc                        <  2:426460/148‑1 (MQ=255)                                                                                                                                 ttttt‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGAc                   >  2:732020/1‑151 (MQ=255)                                                                                                                                     t‑ctc‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTa                >  2:317177/1‑150 (MQ=255)                                                                                                                                        tc‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACt             >  2:88764/1‑151 (MQ=255)                                                                                                                                         c‑gcgGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACtt            >  1:715608/1‑151 (MQ=255)                                                                                                                                             ggggattttt‑CCTTTCTTCCATAC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTg           <  2:86203/143‑1 (MQ=255)                                                                                                                                              gggtttttt‑GCTTTTTCCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCa        <  2:973669/146‑1 (MQ=255)                                                                                                                                                ggATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaa       >  1:47111/1‑150 (MQ=255)                                                                                                                                                  ttttg‑GCTTTCTCCCTTTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAg     <  2:815888/149‑1 (MQ=255)                                                                                                                                                   tttt‑GCTTTCTCCCTTTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaaa      <  1:309802/148‑1 (MQ=255)                                                                                                                                                   tttt‑GCTTTCTCCCTTCC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCaaa      <  2:136331/148‑1 (MQ=255)                                                                                                                                                   tttt‑GCTTTCTCCAATGA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCa        <  1:620022/146‑1 (MQ=255)                                                                                                                                                    ttt‑GCTTTCTCCCATTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCt   <  2:503772/150‑1 (MQ=255)                                                                                                                                                     tt‑GCTTTTTCCCATTC‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTg  <  1:189310/150‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                   CAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAC‑TGTCAACTGCTGAATTTTTTT‑CTC‑GCGGCGGATTTT‑GCGTTCTTCCATTA‑CCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTG  >  NC_000913/3845674‑3845968

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