Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,001,918	0	A	C	26.3%	125.7 / 10.4	57	intergenic (‑129/+192)	fliA/fliC	RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein
Reads supporting (aligned to +/- strand):  ref base A (37/5);  new base C (0/15);  total (37/20)
Fisher's exact test for biased strand distribution p-value = 7.03e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAG  >  NC_000913/2001779‑2002051                                                                                                                                            |                                                                                                                                      agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaaa                                                                                                                                       >  2:875690/1‑140 (MQ=255) agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaaa                                                                                                                                       >  1:74517/1‑140 (MQ=255)   tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                     >  2:748061/1‑140 (MQ=255)     aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                   >  2:933742/1‑140 (MQ=255)      aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                  >  1:609487/1‑140 (MQ=255)      aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                  >  2:1125398/1‑140 (MQ=255)      aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                  >  1:206714/1‑140 (MQ=255)       ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                 >  1:1027672/1‑140 (MQ=255)              aaaaCAGCCCTGCGTTATATGAGTTCTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGAAGAAAAAAACCCCGCCGGTGGc                                                                                                                          <  1:515426/139‑1 (MQ=255)                  cAGCCCTGCTTTTTATGGGTTTTCGGAATTTTTATCCTTTTCTGGGGGGTTTTTAATGGGTCAATTAGGGGGAAAAGTGAGGGGTTGTTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCgggg                                                                                                                      <  2:442233/140‑1 (MQ=255)                      cctggtttATATGAGTTATCGGCATGATTATCCTTTTCGGCGGGTTTTTTAATCGGCCGTTTGGGGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                  <  2:206714/134‑1 (MQ=255)                        tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc                                                                                                                >  1:442233/1‑140 (MQ=255)                            ttATATGATTTACGGGCATGTTTCCCCTTTTCGGCGGGTTTTTAAATGGGACGTTTGGGGGGTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                            <  2:451765/140‑1 (MQ=255)                               tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                         >  1:451765/1‑140 (MQ=255)                                    gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                    >  1:701748/1‑140 (MQ=255)                                    gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                    >  2:711409/1‑140 (MQ=255)                                      tATGGGCATGATTATCCGTTTCTGGGGGGTTTTTAATGGGACGATTAGGGGGTAAAATGGGGGGTTTTTGGGGGGTTAAAGGTAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                  <  1:117763/140‑1 (MQ=255)                                      tATCGGCATGTTTTTCCTTTTTCGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                  <  2:1002086/140‑1 (MQ=255)                                       aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                 >  1:1002086/1‑140 (MQ=255)                                          ggCATGATTATCCTTTTCTGCGGGTTTTTTAATGGGACGATTGGGGGGTAAAATGGGGGGTTTTTGGGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                              <  2:77904/140‑1 (MQ=255)                                           gCATGATTAGCCTTTTCGGCGGGGTTTTAAATCGGATGATGGGGGGGTGAAGTGGGGGTTTATTGGGGGGTTACAGGTAAATCCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                             <  2:612590/140‑1 (MQ=255)                                             aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAGTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                           >  1:612590/1‑140 (MQ=255)                                             aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                           >  2:1539338/1‑140 (MQ=255)                                             aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                           >  1:1518834/1‑140 (MQ=255)                                                      cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGccc                                                                                  >  1:77904/1‑139 (MQ=255)                                                         ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                               >  2:794774/1‑140 (MQ=255)                                                          tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                              >  2:331482/1‑140 (MQ=255)                                                           cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                             >  1:1243111/1‑140 (MQ=255)                                                                      ttAATCGGAGGATTAGTGGGTAAAATGAGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                  <  1:933742/140‑1 (MQ=255)                                                                            ggACGATTGGTGGGAAAAAGGGGGGTTTTTTGGGGGGTTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                            <  2:1509453/140‑1 (MQ=37)                                                                             gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                           >  1:1003476/1‑140 (MQ=255)                                                                             gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                           >  1:267424/1‑140 (MQ=255)                                                                              aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                          <  1:1169355/140‑1 (MQ=255)                                                                                            aaaTGGGGGTTTATTGGGGGTTTACAGGTAATTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                            <  2:1768100/140‑1 (MQ=37)                                                                                            aaaTGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                            <  2:689962/140‑1 (MQ=38)                                                                                            aaaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                        >  2:466469/1‑128 (MQ=255)                                                                                            aaaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                        <  1:466469/128‑1 (MQ=255)                                                                                            aaaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                            >  2:1457008/1‑140 (MQ=255)                                                                                              aggAGGGGTTATTTGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                          <  2:701748/138‑1 (MQ=38)                                                                                                gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                        >  1:1768100/1‑140 (MQ=255)                                                                                                gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                        >  1:689962/1‑140 (MQ=255)                                                                                                       tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                 >  1:789250/1‑140 (MQ=255)                                                                                                        aTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                >  2:1531460/1‑140 (MQ=255)                                                                                                        aTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                >  2:1284884/1‑140 (MQ=255)                                                                                                           tGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                             >  2:879181/1‑140 (MQ=255)                                                                                                           gggggcTTACAGTTAAATTCCGGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                             <  1:748061/139‑1 (MQ=255)                                                                                                                gTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                        >  2:1338148/1‑140 (MQ=255)                                                                                                                gTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                        >  2:869220/1‑140 (MQ=255)                                                                                                                  tACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                      <  1:1457008/140‑1 (MQ=255)                                                                                                                   aCAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                     >  2:371919/1‑140 (MQ=255)                                                                                                                       gTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                 >  2:962578/1‑140 (MQ=255)                                                                                                                       gTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                 >  1:1074935/1‑140 (MQ=255)                                                                                                                         aattgccAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg               <  2:762329/135‑1 (MQ=37)                                                                                                                         aaaTTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg               <  1:794774/140‑1 (MQ=255)                                                                                                                         aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg               >  2:1071737/1‑140 (MQ=255)                                                                                                                          aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt              >  2:706651/1‑140 (MQ=255)                                                                                                                          aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt              >  1:481327/1‑140 (MQ=255)                                                                                                                           aTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa             <  1:1539338/140‑1 (MQ=255)                                                                                                                              ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc          >  2:795899/1‑140 (MQ=255)                                                                                                                                   cAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt     <  1:179800/140‑1 (MQ=255)                                                                                                                                      aaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag  <  2:1003476/140‑1 (MQ=255)                                                                                                                                      aaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag  <  1:1284884/140‑1 (MQ=255)                                                                                                                                      aaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag  >  2:1577660/1‑140 (MQ=255)                                                                                                                                      aaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag  >  1:742423/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                      AGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAG  >  NC_000913/2001779‑2002051

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