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%	117.5 / 21.5	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 (41/1);  new base C (0/15);  total (41/16)
Fisher's exact test for biased strand distribution p-value = 7.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGC  >  NC_000913/2001779‑2002054                                                                                                                                            |                                                                                                                                         agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaaa                                                                                                                                          >  1:151630/1‑140 (MQ=255)  gtaaaTTCACAAAAAACACCCCTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAACCGGACGATTGGGGGGTGAAATGGGGGGTTTTTGGGGGAATACGGGTAAATCCCAGGCA‑AAAAAAAcc                                                                                                                                        <  2:92386/137‑1 (MQ=255)   tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                        >  1:1121092/1‑140 (MQ=255)      aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                     >  1:641314/1‑140 (MQ=255)      aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                     >  2:428813/1‑140 (MQ=255)      aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                     >  1:223020/1‑140 (MQ=255)           cGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGt                                                                                                                                >  1:459371/1‑140 (MQ=255)             aaaaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGGTTTTAAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTgg                                                                                                                              <  2:735612/138‑1 (MQ=255)                      ccTGCTTTATATGAGTTATGGGCATAATTTTCCTTTTCTGCAGGGTTTTAAATCGGACGTTTAGTGGGTAAAATGGGGGTTTTTTGGGGGGTTACAGAAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                     <  2:800421/140‑1 (MQ=255)                       cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCa                                                                                                                    >  1:859481/1‑140 (MQ=255)                          cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGGAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                 >  2:1154233/1‑140 (MQ=255)                             tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                              >  1:58874/1‑140 (MQ=255)                               tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                            >  2:761604/1‑140 (MQ=255)                               tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                            >  1:434061/1‑140 (MQ=255)                                    gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                       >  2:64056/1‑140 (MQ=255)                                    gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                       >  1:498013/1‑140 (MQ=255)                                    gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                       >  1:953143/1‑140 (MQ=255)                                       aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                    >  1:733278/1‑140 (MQ=255)                                       aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                    >  1:40401/1‑140 (MQ=255)                                       aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                    >  1:597245/1‑140 (MQ=255)                                         cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                  >  1:714596/1‑140 (MQ=255)                                         cGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                  >  2:610301/1‑140 (MQ=255)                                                  tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                         >  1:213456/1‑140 (MQ=255)                                                         ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                  >  1:655332/1‑140 (MQ=255)                                                          tcggCAGGTTTTTTAATCGGAGGATTAGGGGGTAAAATGAGGGTTTATTTGGGGGTTACAGGAAAATCCAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                 <  2:223020/137‑1 (MQ=255)                                                          tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                 >  2:1043937/1‑140 (MQ=255)                                                          tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                 >  2:623654/1‑140 (MQ=255)                                                                  gTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                         >  1:614442/1‑140 (MQ=255)                                                                             gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                              >  2:973613/1‑140 (MQ=255)                                                                             gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                              >  1:700303/1‑140 (MQ=255)                                                                             gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                              >  1:521018/1‑140 (MQ=255)                                                                              aCGATTGGGGGGTAAAAGGGGGGGTTTTTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                             <  2:641314/140‑1 (MQ=37)                                                                              aCGATTAGGGGGTGAAATGAGGGGTTATTGGGGGGTTACAGAAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                             <  2:655332/140‑1 (MQ=39)                                                                                        ggTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGTGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                   >  2:309510/1‑140 (MQ=255)                                                                                          tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                 >  1:317144/1‑140 (MQ=255)                                                                                          tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                 >  1:191698/1‑140 (MQ=255)                                                                                          tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                 >  1:408989/1‑140 (MQ=255)                                                                                           gAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                >  1:426818/1‑140 (MQ=255)                                                                                            aaaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                               >  2:463678/1‑140 (MQ=255)                                                                                             aaTGAGGGTTTTTTTGGGGGTTACAGGAAAATTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                              <  2:191698/140‑1 (MQ=37)                                                                                                gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTCCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGaaa                                           >  2:598272/1‑140 (MQ=255)                                                                                                       tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                    >  1:336873/1‑140 (MQ=255)                                                                                                           ggggggTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCCACAAAAATGTGCAAATTCAATAAAt                                <  2:1121092/139‑1 (MQ=37)                                                                                                              gggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                             >  2:129779/1‑140 (MQ=255)                                                                                                               gtttACAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGTAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                            <  2:953143/138‑1 (MQ=25)                                                                                                               ggTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                            <  1:623654/140‑1 (MQ=255)                                                                                                               ggTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                            <  2:58874/140‑1 (MQ=255)                                                                                                               ggTTACAGGAAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                            <  1:1155316/140‑1 (MQ=255)                                                                                                                   aCAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                        >  2:1197780/1‑140 (MQ=255)                                                                                                                         aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                  >  2:1095505/1‑140 (MQ=255)                                                                                                                         aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                  >  1:205879/1‑140 (MQ=255)                                                                                                                           aTTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                <  1:203272/140‑1 (MQ=255)                                                                                                                           aTTCCAGGCAGAAAAAAACCCCGCCGGTGGTGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                <  1:309510/140‑1 (MQ=255)                                                                                                                              ccAGTCATAAAAACCCACCGCCGGTGGCGGGGAAGCGCGGTACTGGCAAACTGCGCTTTATGGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc             <  2:436675/140‑1 (MQ=21)                                                                                                                              ccAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc             <  1:761604/140‑1 (MQ=255)                                                                                                                              ccAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc             <  1:428813/140‑1 (MQ=255)                                                                                                                              ccAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc             <  2:841776/140‑1 (MQ=255)                                                                                                                                 ggCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGGCTTATCCGGCCTACAAAAATGTGGAAATGCAATAAATTGCAATTCAACTTGTAGGCCTg          >  1:82946/1‑140 (MQ=255)                                                                                                                                 ggCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg          >  2:1149709/1‑140 (MQ=255)                                                                                                                                 ggCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg          >  1:1148726/1‑140 (MQ=255)                                                                                                                                      aaaaaaaCCCCGCCGGTTGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag     >  1:1089304/1‑140 (MQ=255)                                                                                                                                         aaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc  <  2:426818/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                         AGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGC  >  NC_000913/2001779‑2002054

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