Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,271	0	T	G	42.6%	97.4 / 10.4	62	G165G (GGT→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base T (4/31);  new base G (26/0);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 2.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGC  >  NC_000913/1157122‑1157416                                                                                                                                                      |                                                                                                                                                    catATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGt                                                                                                                                                     <  1:754028/150‑1 (MQ=255) catATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGt                                                                                                                                                     >  2:36351/1‑150 (MQ=255)             ggTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATaat                                                                                                                                             >  1:643659/1‑144 (MQ=255)             ggTCGTGAACTGAACAAGCCGGTTAACGTCCAGACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGAGTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGAtc                                                                                                                                        <  2:153707/151‑1 (MQ=255)                   gaactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGAt                                                                                                                                  <  2:926038/151‑1 (MQ=255)                     actgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGGTCTCGGGttt                                                                                                                                >  1:738877/1‑151 (MQ=255)                     actgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGAGATCGGGttt                                                                                                                                >  2:344321/1‑151 (MQ=255)                        gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTa                                                                                                                              <  2:177542/150‑1 (MQ=255)                          acAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACAt                                                                                                                           <  1:934170/151‑1 (MQ=255)                          acAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGTCTCGGGTTTTaaa                                                                                                                            >  2:71095/1‑148 (MQ=255)                              gCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑cc                                                                                                                      <  2:348673/151‑1 (MQ=255)                                 ggTTATCGTCCATGCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑ttt                                                                                                                  <  2:131062/151‑1 (MQ=255)                                 ggTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑ttt                                                                                                                  <  2:1109581/151‑1 (MQ=255)                                  gTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGGACTCGGATTTTatcattactttt                                                                                                                   >  2:343655/1‑140 (MQ=255)                                   ttATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTTGATCTCGGGGTTTATATCT‑CC‑tttt                                                                                                                 >  2:1090517/1‑150 (MQ=255)                                     aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCg                                                                                                              <  2:972551/151‑1 (MQ=255)                                          ccATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCAtt                                                                                                         <  2:911364/151‑1 (MQ=255)                                            aTACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATtgt                                                                                                       <  1:305858/151‑1 (MQ=255)                                             tACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTTCCGCAtttt                                                                                                       >  1:106548/1‑148 (MQ=255)                                             tACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTTGGTCTCGGATTTTTCATCT‑CC‑CTTTTCCGCAtttt                                                                                                       >  1:385825/1‑148 (MQ=255)                                                   tGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCtt                                                                                                 <  1:869519/150‑1 (MQ=255)                                                   tGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATAATGGAACTCGGGTTTTTCATAT‑CC‑TTTTTCCGCaatttggc                                                                                                    >  2:701381/1‑140 (MQ=255)                                                   tGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATTATGGTTACCGGGCTTTTAATTT‑CC‑TTTTTCGGCACtgttt                                                                                                     >  2:269268/1‑144 (MQ=255)                                                    gACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTcc                                                                                               <  2:985997/151‑1 (MQ=255)                                                     aCGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCg                                                                                              <  2:851698/151‑1 (MQ=255)                                                     aCGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTATTGGGGCTTCGATTTTATATCT‑CC‑TTTTTTCGCATTGTGGCCtttc                                                                                               >  2:952028/1‑148 (MQ=255)                                                       gCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGCGAATTAATGGGTCTCGGCTTTTACATCT‑CC‑TTTTTCTGTATtgtggacctcc                                                                                               >  1:511558/1‑141 (MQ=255)                                                        cccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTTTATCT‑CC‑TTTTTCGGCATTGTGACCTTCCGTa                                                                                            >  1:195496/1‑150 (MQ=255)                                                         ccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGGTCTCGGATTTTACATCT‑CC‑TTTTTCCGCATTTTGGCCTTCCGTaa                                                                                           >  2:778253/1‑150 (MQ=255)                                                           gTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATAATGGATCTCGGGTTTTTCAACT‑CC‑CTTTTCTGCCTTTTGACCtttcctta                                                                                           >  2:187159/1‑142 (MQ=255)                                                               cGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGATCTCGGGTTTTACATCT‑CC‑CTTTTCGGCATTGTGACCTTCCGTAATGCgg                                                                                      >  1:485376/1‑149 (MQ=255)                                                               cGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATCTCGGGTTTTACATCT‑CC‑TTTTTCCGCATTGTGACCTTCCCTaaatcgggg                                                                                    >  1:721770/1‑144 (MQ=255)                                                                 aTACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCa                                                                                  <  2:206608/151‑1 (MQ=255)                                                                 aTACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGTCTCGGGTTTTTCATCT‑CC‑TTTTTCCGCATTTTGAACTTCCGTAATGCgggg                                                                                    >  1:369721/1‑147 (MQ=255)                                                                     gCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAATTAATTGATCTCGGGGTTTACAACT‑CC‑TTTTTCCGCATTGTGACCTTCCGTAATGCGGGGCAACTg                                                                              >  2:272477/1‑151 (MQ=255)                                                                         gCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGGGCTCGGGTTTTATATCT‑TC‑TTTTTCCGCGTTGTGGCCCTCCCTTAAGGGGGGGAACcgcgcg                                                                          >  2:426525/1‑151 (MQ=255)                                                                         gCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGGTCTCGGATTTTACATCT‑CC‑TTTTTCGGCATTGTGACCTTCCGTAATGCGGAGCCAATGcgcg                                                                          >  1:830817/1‑151 (MQ=255)                                                                          cAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGa                                                                         <  1:179187/151‑1 (MQ=255)                                                                           aaTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGa                                                                         <  2:537822/150‑1 (MQ=255)                                                                           aaTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATCTCGGGTTTTACATCT‑CC‑TTTTTCTGCATTGTGACCTTCCGCAATGCGGGGCAAATgcgc                                                                           >  1:667833/1‑148 (MQ=255)                                                                            aTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGATCTCCGGTTTTACATCT‑CC‑TTTTTCCGCATTGTGAACTTCCGTAATGCGGAGAAACCGCGCGGtg                                                                       >  1:704588/1‑151 (MQ=255)                                                                                 gcgcGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTgcg                                                                  <  2:389327/151‑1 (MQ=255)                                                                                  cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATACTGGATCTCGGGGTTTACATCT‑CC‑TTTTTCGGCATTTTTACCCTCCGTAATGTGGAGCAACAGCGCGATGCTgcg                                                                  >  1:1020046/1‑150 (MQ=255)                                                                                   gcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTgcgc                                                                 <  2:3652/150‑1 (MQ=255)                                                                                    cgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTTCTGGATCTCGGGTTTTAAATCT‑CC‑TTTTTCGGCATTGTGAGCTTTCCTAATGCGGGGGAACTGCGCGATGCTGcgcg                                                                >  2:491450/1‑150 (MQ=255)                                                                                    cgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGATCTCGGGTTTTAAATCT‑CC‑CTTTTCGGCATTGTGAACTTCCGGAAAGCGGGGGAACTGcgcggt                                                                        >  2:439808/1‑140 (MQ=255)                                                                                        gaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATg                                                           <  1:541177/151‑1 (MQ=255)                                                                                            aaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTccc                                                       <  1:22754/151‑1 (MQ=255)                                                                                               tGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTTCATCT‑CC‑TTTTTCGGCCATGTGACCTTCCGTTATGCGGGGGAACCGCGCGGTGGTGCGCGCTTTGTccccc                                                     >  1:1031846/1‑150 (MQ=255)                                                                                                  cGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATTTTGGGTCTCGGGTTTTTTATTT‑CC‑TTTTTTGGCCTTGTGACCTTCCCTAATGCGGGGCAACTGCGCGGTGCTGCGCGGTTTGTTCCCCCgg                                                  >  2:461086/1‑150 (MQ=255)                                                                                                   ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGAt                                                <  2:451960/151‑1 (MQ=255)                                                                                                     aTTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCg                                              <  2:690465/151‑1 (MQ=255)                                                                                                       tGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCgg                                             <  1:550261/150‑1 (MQ=255)                                                                                                        gCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGtt                                           <  2:1073790/151‑1 (MQ=255)                                                                                                        gCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGAAATCGGGGTTTAAATCTCCC‑TTTTCCGGC‑TTGTGAGCTTTCGTAAAGCGGAGCAACAACGCGGAGCTGCGCGGTATTTcccccccgatagggt                                           >  1:36832/1‑141 (MQ=255)                                                                                                             ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACtg                                       <  2:619584/150‑1 (MQ=255)                                                                                                             ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACtg                                       <  1:694527/150‑1 (MQ=255)                                                                                                              gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACtggt                                     <  2:738851/151‑1 (MQ=255)                                                                                                               cGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATAATGGATATCCGGGTTTACATCT‑CCTTTTTTCGGGCTT‑TGAGCTTCCCTTAAGCGGGGGAAACAGGCGATGCCTGGCGGTTTTTTCCCCTGGATCGGTTACttg                                      >  1:615695/1‑147 (MQ=255)                                                                                                                     acacTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGaaa                                <  2:195487/149‑1 (MQ=255)                                                                                                                      cacTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTg                             <  2:221158/151‑1 (MQ=255)                                                                                                                       acTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATTTCGGATTTTAAATCC‑CC‑TTTTTCGGCATTTTGGGCTCCCGTCATGCGGGGGAATATCGCGCTGCTGCGCGGTTTTTCCCCCCGGATCGGGTTCTGGTGGgaa                                >  2:877883/1‑147 (MQ=255)                                                                                                                                            aaaCGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCACCTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGt       <  2:574269/151‑1 (MQ=255)                                                                                                                                                 gcggGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGc  <  2:125587/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CC‑TTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGC  >  NC_000913/1157122‑1157416

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