Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,001,911	0	G	A	40.7%	77.0 / 28.4	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 (32/3);  new base A (0/24);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.35e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001769‑2002059                                                                                                                                                    |                                                                                                                                                      aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑c                                                                                                                                              <  2:173861/151‑1 (MQ=255)    ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGATATCGGCATGATTATCCGTTTCAGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGCGAGATGCGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATGCCAGGCAGAAGAAAA‑accc                                                                                                                                           >  2:60935/1‑151 (MQ=255)     gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCccg                                                                                                                                          >  1:31502/1‑151 (MQ=255)            gTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTgg                                                                                                                                   >  1:918055/1‑151 (MQ=255)                aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCggg                                                                                                                               >  1:702438/1‑151 (MQ=255)                  tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCgggg                                                                                                                              >  1:150213/1‑150 (MQ=255)                       aaaaaCGGCCCTGCTTTATATGGTTTATCGGCTTGTTTACCTTTTTCTGGGGGTTTTTT‑AGTCGAAC‑TATTGGTGGGGAAAAGGGGGGTTTATTGGGGGG‑TT‑ACAGG‑TAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAg                                                                                                                           <  2:873941/146‑1 (MQ=255)                       aaaaaCAGCCATGCGTTATATATGTTAGGGTCATTATTATCCTTTTCTGCAGGTTTTTT‑AATCGGAT‑GTTTAGGGGGTAAAAAGAGGGGTTTTTGGGGGT‑TT‑ACAGA‑TAATTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGc                                                                                                                          <  1:293510/147‑1 (MQ=255)                       aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCAc                                                                                                                        >  1:45075/1‑151 (MQ=255)                         aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                      >  1:51984/1‑151 (MQ=255)                         aaaCAGCCCGTCGTTATATGGTTTAGGGGATTTATTCCCCGTTTCCGCGGGTTTTTT‑AAGCGAAG‑GATTGGTGGGTAAAATGAGGGGTTATTGGGGGG‑TT‑ACAGG‑TAAATCCAAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACg                                                                                                                       <  1:945504/150‑1 (MQ=255)                              gCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                 >  1:204382/1‑151 (MQ=255)                              gCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                 >  1:615747/1‑151 (MQ=255)                                   gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                             >  1:903000/1‑150 (MQ=255)                                    cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                           >  1:918266/1‑151 (MQ=255)                                     gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                          >  2:716925/1‑151 (MQ=255)                                      ttATTTGAGTTATGGGTATGTTTATCCGTTTCTG‑AGGGTTTTT‑TAACCGACAGATAAGGGGGTGAAAGGGGGGTTTTTTTGGGTT‑TC‑ACGGG‑AAATTTCAAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                          <  1:534395/150‑1 (MQ=255)                                       tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                        >  1:519503/1‑151 (MQ=255)                                       tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                         >  2:242767/1‑150 (MQ=255)                                         tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                       >  1:143603/1‑150 (MQ=255)                                                 atgggCTTGATTTTCCGTTTCGGCAGGTTTTTT‑AATGGGAG‑GATTGGGGGGTAAAAGGGGGGTTTTTTGGGGGG‑TA‑ACAGA‑TAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                              <  2:902876/148‑1 (MQ=255)                                                 aTCGGCATGATTATCCTTTTCTGCGGGTTTTTT‑AATCGGAC‑GATTGGGGGGTAAAATGGGGGGTTTTTTGGGGT‑TT‑ACAGG‑TAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                               <  1:771322/150‑1 (MQ=255)                                                      cATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         >  1:13587/1‑151 (MQ=255)                                                        tGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                        >  2:762216/1‑150 (MQ=255)                                                           ttATCGTTTTCTGGAGGTTTTTT‑AAGGGGAC‑AATTAGGGGGAAAAATGGGGGTTTTTTTGGGGT‑TT‑AAAGT‑AAAATCCGAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         <  1:944437/146‑1 (MQ=255)                                                               ccttttCGGCTGGTTTTTT‑AATCGGAC‑ATTTGGTGGGTAAAAGGGGGGTTTATTGGGGGG‑TT‑‑CGGG‑TAAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                  <  2:496822/146‑1 (MQ=255)                                                                  tttCTGCAGGGTTTTT‑AATGGGAC‑GATTGGGGGGTAAAATGGGGGTTTTTTTGGGGT‑TT‑ACAGG‑AAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                               <  2:913999/149‑1 (MQ=255)                                                                   ttCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                            >  1:718772/1‑151 (MQ=255)                                                                        cAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       >  1:187050/1‑151 (MQ=255)                                                                         cgggTTTTTAAATCGGAG‑ATTTAGGGGGT‑AAAAGAAGGGTTTTTTGGGGG‑TTAACAGT‑TAATTCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       <  2:34558/150‑1 (MQ=255)                                                                          gggTTTTA‑AATGGGAC‑GATTAGGGGGTAAAATGAGGGTTTATTGGGGGG‑TT‑ACAGG‑TAAATACAGGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                      <  1:580908/150‑1 (MQ=255)                                                                             ttttt‑AATGGGAC‑GTTTGGGGGGTGAAATGAGGGTTTTTTGGGGGG‑TT‑ACAGT‑TAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                   <  1:877077/150‑1 (MQ=255)                                                                              tttt‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                 >  1:599052/1‑151 (MQ=255)                                                                              tttt‑AATCGGAC‑GATTAGTGGGAGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGAGGCGGGGAAGCCCGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                  >  2:869262/1‑150 (MQ=255)                                                                                 t‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                              >  1:471257/1‑151 (MQ=255)                                                                                   aaTCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                              >  1:743644/1‑150 (MQ=255)                                                                                       ggAC‑GATTAGTGGGGAAAATGGGGGTTTTTTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                         <  1:570757/151‑1 (MQ=38)                                                                                         ggaagtttagggggGAAAAGGGGGGTTTTTTGGGGG‑GT‑TGAGGATAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                         <  2:1178/141‑1 (MQ=25)                                                                                         aC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        >  1:814774/1‑150 (MQ=255)                                                                                               tAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                   >  1:391711/1‑150 (MQ=255)                                                                                                 gTGGGTAAAATGAGGGGTTATTTGGGGT‑TT‑ACAGT‑AAAATCCCAGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                 <  1:100599/150‑1 (MQ=38)                                                                                                   gggTAAAATG‑GGGGTTTTTTGGGGGTTT‑ACAGG‑TAATTCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  <  2:79691/147‑1 (MQ=37)                                                                                                   gggTAAAAGGGGGGGTTATTGGGGGG‑TT‑ACAGG‑TAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  <  1:136737/147‑1 (MQ=37)                                                                                                    ggTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAACGTGCAAAtt                                             >  2:40161/1‑151 (MQ=255)                                                                                                     gTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                            >  2:746115/1‑151 (MQ=255)                                                                                                       gAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                           >  2:405069/1‑150 (MQ=255)                                                                                                        aaaGGAGGGTTTTTTTGGGGT‑TT‑ACAGG‑AAAATTCAAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                          <  1:354174/150‑1 (MQ=37)                                                                                                         aaggAGGGTTTATTTGGGGG‑TT‑AAAGG‑AAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                        <  1:605224/148‑1 (MQ=37)                                                                                                            gggggtTTTTTTGGGGG‑TT‑ACAGG‑TAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                     <  1:23727/149‑1 (MQ=37)                                                                                                              ggggTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                   >  1:684376/1‑151 (MQ=255)                                                                                                                   tATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                              >  2:501896/1‑151 (MQ=255)                                                                                                                      ttGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             >  2:528188/1‑149 (MQ=255)                                                                                                                      tgggggg‑TT‑ACAGG‑AAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                            <  1:905279/148‑1 (MQ=255)                                                                                                                        ggggg‑TT‑AAAGG‑AAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                         <  1:319941/151‑1 (MQ=37)                                                                                                                        ggggg‑TA‑ACATG‑AAAATACCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                         <  1:372431/151‑1 (MQ=37)                                                                                                                         gggg‑TT‑AAAGG‑TAAATTCCAGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                            <  1:315207/148‑1 (MQ=255)                                                                                                                                        aaaTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt             <  1:696422/150‑1 (MQ=255)                                                                                                                                             ccAGGGAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa           <  2:636890/147‑1 (MQ=255)                                                                                                                                                  cAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  1:1009074/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGAC‑GATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑ACAGG‑TAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001769‑2002059

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