Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,001,918	0	A	C	44.8%	90.9 / 14.5	58	intergenic (‑129/+192)	fliA/fliC	RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (25/7);  new base C (0/26);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 2.67e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001768‑2002066                                                                                                                                                        |                                                                                                                                                     cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTGATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaaa                                                                                                                                                      <  2:499923/151‑1 (MQ=255)    cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAccc                                                                                                                                                   <  2:374929/151‑1 (MQ=255)     ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                                  >  2:785786/1‑151 (MQ=255)        ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGCTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc                                                                                                                                               >  2:689160/1‑151 (MQ=255)          agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                             <  2:593477/151‑1 (MQ=255)           gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                             >  1:750814/1‑150 (MQ=255)                 aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                                      >  1:469308/1‑151 (MQ=255)                    cACGATAAACGGCCCTGCGTTATATGAGTTATCGGCATGATTTCCCGTTTCTGCAGGGTTTTT‑AATCGGACGTTTAGTGGGTGAAATGAGGGTTTATTGGGGGGTTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                   <  1:497153/151‑1 (MQ=255)                    cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCGGCAGGTTTTTT‑AATCGGACGATTGGTGGGTAAAAGGAGGGGTTATTTGGGGTTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCgggg                                                                                                                                     <  1:166267/149‑1 (MQ=255)                        aaaaaCAGCCCTGCTTTATATGAGTTATCGGAATGATTACCCTTTTCTGCAGGTTTTTT‑AATCGGACGTTTAGTGGGTGAAAGGAGGGTTTTTTGGGGGTTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCa                                                                                                                                <  1:308118/148‑1 (MQ=255)                          aaaCAGCCTTGTTTTAAAGGGTTTTTCGGATTGTTTTTCCTTTTCTGCGGGTTTTTT‑AATCGGACGTTTAGTGGGTGAAAGGGGGGTTTTTTGGGGGGTTACAGTAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                                 <  1:730720/147‑1 (MQ=255)                           aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTACCCGTTTCTGCGGGTTTTTT‑AATCGGACGATTAGGGGGTGAAAGGAGGGTTTTTTTGGGGGTTAAAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                            <  2:660591/151‑1 (MQ=255)                                cccTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTT‑AATCGGACGATTAGGGGGTGAAATGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                       <  1:90914/151‑1 (MQ=255)                                cccTGCGTTATATGATTTATCGGAATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGTTTACGGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                                          <  1:92993/148‑1 (MQ=255)                                   tGCGTTATATGAGTTATCGGCATGATTATCCTTTTCTGAGGGGTTTTG‑AATCGGAAGAGTTGTGGGTAAAATGAGGGTTTATGTGGGGGTCACGGGAAATTCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                        <  1:25520/147‑1 (MQ=255)                                       ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                >  1:444682/1‑151 (MQ=255)                                          tatGAGTTTGCGGAATGATTATCCGTTTCGGCGGGTTTTTT‑AATCGGACGTTTAGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                               <  1:903672/149‑1 (MQ=255)                                             gatttAGCGGCTTGTTTTTCCGTTTCTGCAGGTTTTTT‑AACCGGACGATTGGTGGGTAAAATGGGGGTTTTTTTGGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                <  1:896509/142‑1 (MQ=255)                                             gAGTTCTCGGCATGTTTACCCTTTTCTGCAGGTTTTTT‑AATCGGACGATTAGGGGGTAAAAGGGGGGTTTTTTTGGGGGTAACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                           <  2:543271/150‑1 (MQ=255)                                              aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                          >  2:482555/1‑150 (MQ=255)                                                ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                        >  2:51728/1‑150 (MQ=255)                                                ttAGCGGCATGATTATCCTTTTCTGCAGGTTTTTT‑AATGGGACTTTTAGGGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGATAATTCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                          <  2:895945/148‑1 (MQ=255)                                                 tttCGGCTTGATTTTCCTTTTCTGAGGGTTTTTT‑AATCGGAGGTTTTGGGGGGAAAAATGGGGGTTTTTTGGGGGTTAAAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                            <  1:851595/143‑1 (MQ=255)                                                             tccctgTTTCTGCAGGGTTTTA‑AATCGGACGATTAGGGGGGGAAATGAGGGTTTATTTGGGGGTTACAGGTAATTTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                           <  2:789164/145‑1 (MQ=255)                                                                   tttCTGCGGGTTTTTT‑AATCGGACGTTTGGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                        <  2:393508/147‑1 (MQ=255)                                                                    ttCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                   >  1:254808/1‑151 (MQ=255)                                                                       tGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                                >  2:903510/1‑151 (MQ=255)                                                                         gaGGGTTTTTAAAGCGAACGATAGGTGGAAAAAAGGAGGGGTTGTTTGGGGGTACTAGAAAAATCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                   <  1:689256/146‑1 (MQ=255)                                                                         cAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              >  2:897442/1‑151 (MQ=255)                                                                          aGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              >  2:330030/1‑150 (MQ=255)                                                                           gggTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                             >  2:411025/1‑150 (MQ=255)                                                                           gggTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                            >  1:533905/1‑151 (MQ=255)                                                                           gggTTTTT‑AATCGGACGATGAGGGGGAAAAAGGGGGGGTTTTTTGGGGGTAACAGGAAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                              <  2:828434/149‑1 (MQ=255)                                                                              ttttt‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                         >  2:94227/1‑151 (MQ=255)                                                                              ttttt‑AACGGGACGATTGGGGGGTAAAAGGGGGGGTTTTTTGGGGGTTACAGGTAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                           <  2:351823/149‑1 (MQ=255)                                                                                ttt‑AAGCGGACGTTTGGGGGGTAAAAGGGGGGTTTTTTTGGGGGTTAAAGTTAAATCCCGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                        <  2:403405/150‑1 (MQ=25)                                                                                    aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                     >  1:926349/1‑150 (MQ=255)                                                                                      tCGGAGGATTGGTGGGTAAAAGGAGGGGTTATTGGGGGTTTAAGGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                   <  1:593525/150‑1 (MQ=37)                                                                                      ccGGACGATTAGGGGGTAAAATGGGGGTTTTTTTGGGGTTTAAAGGAAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                  <  2:539486/150‑1 (MQ=25)                                                                                            gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                            >  1:105149/1‑151 (MQ=255)                                                                                             aTTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                            >  2:811923/1‑150 (MQ=255)                                                                                              ttggggggTAAAAGGGGGGTTTTTTTGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                             <  2:122933/143‑1 (MQ=25)                                                                                               tAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                         >  2:315484/1‑151 (MQ=255)                                                                                                   gggTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                     <  1:643897/151‑1 (MQ=255)                                                                                                   gggTGAAAGGAGGGTTTTTTTGGGGTTTAAAGTTAAATCCCAGACAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                     <  1:240158/151‑1 (MQ=37)                                                                                                   gggTAAAATGAGGTGTTATTGGGGGTTTACAGGTAATTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                        <  1:868346/148‑1 (MQ=38)                                                                                                      tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                   >  2:466912/1‑150 (MQ=255)                                                                                                            ggggggTTTTTTGGGGTTTAAAGTAAAATTCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                             <  2:837228/148‑1 (MQ=25)                                                                                                             aGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                           >  1:490507/1‑151 (MQ=255)                                                                                                              gggtTTTTTTGGGGTTTACAGGAAAATTCCAGAAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                           <  2:585858/150‑1 (MQ=37)                                                                                                                     tttGGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                                      <  1:118141/148‑1 (MQ=255)                                                                                                                     tttGGGGGTTACAGTTAAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                    <  1:660764/150‑1 (MQ=255)                                                                                                                       tGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                  >  2:428689/1‑150 (MQ=255)                                                                                                                           ggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                              >  1:686513/1‑150 (MQ=255)                                                                                                                              tACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                          >  2:31206/1‑151 (MQ=255)                                                                                                                               aCAGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                           <  1:409450/149‑1 (MQ=255)                                                                                                                                            aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc              >  2:740113/1‑149 (MQ=255)                                                                                                                                             ggaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag            <  2:76026/145‑1 (MQ=255)                                                                                                                                             ggaaGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca             <  2:533819/146‑1 (MQ=255)                                                                                                                                                 gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt       >  1:22915/1‑151 (MQ=255)                                                                                                                                                   aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa     <  1:636280/151‑1 (MQ=255)                                                                                                                                                       aaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc  >  1:389919/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                     CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGC  >  NC_000913/2001768‑2002066

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