Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,001,911	0	G	A	44.1%	56.5 / 39.8	59	intergenic (‑122/+199)	fliA/fliC	RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (31/2);  new base A (0/26);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 9.54e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGC  >  NC_000913/2001764‑2002054                                                                                                                                                        |                                                                                                                                                 ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGaaa                                                                                                                                               <  1:446645/151‑1 (MQ=255)           tataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCccgc                                                                                                                                    >  1:740/1‑151 (MQ=255)              agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCgg                                                                                                                                 >  2:660824/1‑151 (MQ=255)              agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCgg                                                                                                                                 >  2:566024/1‑151 (MQ=255)                  tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGc                                                                                                                             >  1:937253/1‑151 (MQ=255)                         aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGa                                                                                                                        >  1:617248/1‑149 (MQ=255)                            aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCa                                                                                                                    >  1:300700/1‑150 (MQ=255)                              aaaCGGCCCTGCTTTATATGAGTTATCGGCATTATTATCCTTTTCGGCGGGGTTTTTAATCGGACGATTAGTGGGTGAAAGGGGGGGTTATTG‑GGGGGTT‑AAAGAAAAATTCC‑GGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAg                                                                                                                      <  1:680469/146‑1 (MQ=255)                                   ccccTGGGTTATTGGAGTTATCGGCATGTTTACCCTTTTCGGCAGGTTTTTTAATCGGACTTTTAGTGGGTAAAAGGGGGGTTTTTTG‑GGGGGTT‑ACAGGAAAATTCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                            <  2:623147/150‑1 (MQ=255)                                    cccTGTTTTATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTATTG‑GGGGGTT‑ACAGGTAAATCCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                               <  2:423548/147‑1 (MQ=255)                                           ttATATGATTTATCGGAAGGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGGTTTTTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                          <  2:166071/145‑1 (MQ=255)                                            tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTG‑GGGGGTT‑ACAGGTAAATTCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                    <  1:954915/150‑1 (MQ=255)                                            tataGGATTTCTCGGCTTGTTTATCTGTTTCTGGAGGGTTTTAAAGCGGACAATTGGGGGGTGAAAGGGGGGGTTTTTG‑GGGGGTT‑ACAGG‑AAATTCC‑CGGCAAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                    <  1:13967/150‑1 (MQ=255)                                             atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                  >  2:708365/1‑151 (MQ=255)                                               atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                >  1:789514/1‑151 (MQ=255)                                                 gAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                              >  1:905459/1‑151 (MQ=255)                                                     tATCGGCATGATTATCCTTTTCTGAGGGTTTTTTAATGGGACGATTAGGGGGTAAAAGGGGGGGTT‑TTT‑GGGGGTTAAAAGGAAAATCCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                           <  1:905412/150‑1 (MQ=255)                                                      acggccatttttATCCTTTTCTGGAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTTTTG‑GGGGGTT‑AAAGGTAAATCCA‑AGGC‑AAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                          <  1:563474/140‑1 (MQ=255)                                                      aTCGGCATGTTTATCCTTTTCTGCGGGGTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGGTTTTTT‑GGGGGTT‑ACAGGTAATTTCC‑AGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                         <  2:588551/151‑1 (MQ=255)                                                      aTCGGCATGATTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGTTTATTT‑GGGGTTT‑ACAGGTAAATTCC‑AGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                          <  2:513197/150‑1 (MQ=255)                                                      aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                         >  2:704646/1‑151 (MQ=255)                                                      aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGGTATTT‑GGGGGTT‑ACAGGTAAATTCC‑ACGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                          >  2:448157/1‑150 (MQ=255)                                                        cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                       >  2:494405/1‑151 (MQ=255)                                                        cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                       >  1:17022/1‑151 (MQ=255)                                                             tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                  >  2:615809/1‑151 (MQ=255)                                                              gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                 >  2:531028/1‑151 (MQ=255)                                                                ttACCCTTTTTTGCAGGTTTTTTAATCGAACGATAAGGGGGTGAAAGGGGGGTTTTTTG‑GGGGGTT‑AAAGGAAAATCCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                   <  1:960059/147‑1 (MQ=255)                                                                 tttCCTTTTCGGCGGGTTTTTAAAGGGGAAGATTGGTGGGTAAAAGGAGGGTTTATTG‑GGGGGTT‑ACA‑GTAAATTCCAAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                  <  1:672866/145‑1 (MQ=255)                                                                  aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                             >  2:160031/1‑151 (MQ=255)                                                                   tCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                            >  1:456217/1‑151 (MQ=255)                                                                    ccGTTTCTGCGGGGTTTTTAATGGGACGTTTGGGGGGTAAAAGGGGGGGTTTTTG‑GGGGGTT‑AAAGTAAAATCCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                           <  1:924950/151‑1 (MQ=255)                                                                      gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                          >  1:702006/1‑150 (MQ=255)                                                                       ttttTGCGGGTTTTTTAACCGGACGATTAGTGGGTAAAATGGGGGTTTATTG‑GGGGGTT‑ACAGGTAAATCCC‑AGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                           <  2:473642/148‑1 (MQ=255)                                                                               ggttttttAAGGGGACGATTGGGGGGTAAAATGAGGTTTTATTG‑GGGGGTT‑ACAGGTAAATTCC‑AGAA‑AAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                       <  2:34422/141‑1 (MQ=255)                                                                               gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                >  2:224740/1‑151 (MQ=255)                                                                                  tttttAAGCGGGAGATTGGGGGGTAAAATGAGGGGTTTTTG‑GGGGGTG‑AAAGTTAAATCCC‑GGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                               <  1:788618/149‑1 (MQ=255)                                                                                  ttgttAATCGGACGGTTAGGGGGTGAATTGAGGGTTTATTT‑GGGGGGT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAGCGCCTTATCCg                                                             >  1:137967/4‑151 (MQ=255)                                                                                         tCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                       >  1:145217/1‑150 (MQ=255)                                                                                            gACGATTGGGGGGTGAAATGGGGGGTTATTG‑GGGGGTT‑ACAGATAAATCCA‑AGGA‑AGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                   <  1:322042/151‑1 (MQ=37)                                                                                               gATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                >  1:580569/1‑151 (MQ=255)                                                                                                 ttAGTGGGTGAAAAGGGAGTTTCTTG‑GAGGTTT‑GGTGATAAATCCC‑GGGC‑AAAAAAACA‑CGCCGCCGGTGGCGGGGAAGCACGTGGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                               <  1:537876/150‑1 (MQ=21)                                                                                                 ttAGTGGGGGAAATGAGGGGTTTTTT‑GGGGGTT‑AAGGAAAAATTCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                 <  1:217348/148‑1 (MQ=37)                                                                                                    gggggTAAAATGAGGGTTTTTTG‑GGGGGTT‑ACAGGTAAATCCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                            <  1:780997/148‑1 (MQ=37)                                                                                                    gggggTAAAAGGAGGGTTTTTTG‑GGGGTTT‑ACAGGAAATTCCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                              <  1:916954/146‑1 (MQ=25)                                                                                                         tGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                       >  1:777615/1‑150 (MQ=255)                                                                                                            aaaggggggTTTTTTGGGGGGTT‑ACAGTAAAATTCC‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                    <  1:240925/146‑1 (MQ=25)                                                                                                               gAGGGGTTTTTG‑GGGGGTT‑ACAGGTAAATCCC‑AGAA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                 <  1:837275/150‑1 (MQ=37)                                                                                                                 ggggTTATTG‑GGGGGTT‑ACAGGAAAATTCC‑AGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                              <  1:670023/151‑1 (MQ=255)                                                                                                                     ttATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                          >  2:182691/1‑151 (MQ=255)                                                                                                                              gggTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                  >  1:378904/1‑151 (MQ=255)                                                                                                                                      aGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt           >  1:739915/1‑151 (MQ=255)                                                                                                                                        gTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg          >  1:704172/1‑150 (MQ=255)                                                                                                                                          aaaTTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt           >  2:451997/1‑147 (MQ=255)                                                                                                                                           aaatCC‑AGGG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa       <  1:428978/148‑1 (MQ=255)                                                                                                                                               cc‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc  >  1:924608/1‑151 (MQ=255)                                                                                                                                               cc‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc  >  1:359578/1‑151 (MQ=255)                                                                                                                                               cc‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc  <  2:21379/151‑1 (MQ=255)                                                                                                                                               cc‑AGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag     <  2:533973/148‑1 (MQ=255)                                                                                                                                                c‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc  >  1:69443/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                 CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTT‑ACAGGTAAATTCC‑AGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGC  >  NC_000913/2001764‑2002054

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