Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	445,558	0	A	C	30.6%	71.7 / 26.7	49	V370G (GTA→GGA)	yajR	putative major facilitator superfamily transport protein YajR
Reads supporting (aligned to +/- strand):  ref base A (24/10);  new base C (0/15);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 2.90e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTG  >  NC_000913/445425‑445693                                                                                                                                       |                                                                                                                                        aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCGTGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCt                                                                                                                                   <  1:1162804/140‑1 (MQ=255) aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCt                                                                                                                                   <  1:1073194/140‑1 (MQ=255)   gCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGa                                                                                                                                 <  2:1003842/140‑1 (MQ=255)   gCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGa                                                                                                                                 <  2:1191169/140‑1 (MQ=255)            aCGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATAccccc                                                                                                                                    <  1:1178234/128‑1 (MQ=255)            aCGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATAccccc                                                                                                                                    >  2:1178234/1‑128 (MQ=255)            aCGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATAccccc                                                                                                                                    <  1:408667/128‑1 (MQ=255)            aCGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATAccccc                                                                                                                                    >  2:408667/1‑128 (MQ=255)             cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACa                                                                                                                       <  2:746258/140‑1 (MQ=255)              gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                      <  2:836616/140‑1 (MQ=255)                           cAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                         <  1:892055/140‑1 (MQ=255)                              tGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGGACTGGCGGCTGTCAGCCACACTGCGGCCCGCCTT‑TCGCCCGCGAGAAAAACCCCCCCACCGTCAAACATGCCGTTAATCCAgccc                                                                                                      >  2:1491608/1‑139 (MQ=255)                                ttCGCCGGGATTTCAATGCGCAACCTGCTGACATCCGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAgccgcc                                                                                                    <  1:1491608/140‑1 (MQ=255)                                          tttCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCCGCGAACCg                                                                                          >  1:601974/1‑140 (MQ=255)                                          tttCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAAccc                                                                                          >  1:1263664/1‑139 (MQ=255)                                                  cgcAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATT‑GCGCCCGCGAGAAAAACCCCCCGACCGTCAAACATGCCGTTAATCCAGCCGCCCCGCGAACCGCCCAAcgc                                                                                  >  2:391642/1‑140 (MQ=255)                                                    cAAACTGCTGACATACGGCGGTTCTTTCATGGGACTGGCGACTGTCAGCCCCACTCCGGCCCGCATT‑GCGCCCGCGAGGAAAACCCCCTGACCGTCAAAAATGCCGTTAATCCCACCGCCCCGCGGAACCCCaacccccc                                                                                >  2:842064/1‑134 (MQ=255)                                                        ctgctgACATCCGGCGTTTCTTTCAGGGTACGGGCGCCTGTCACCCAAACTGCGCCCAGCATT‑GCGCCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgccacgcc                                                                            <  2:1075391/140‑1 (MQ=255)                                                         tgctgACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑ACGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgccacgcca                                                                           <  1:966478/140‑1 (MQ=255)                                                                   cggcggTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAAAAATACCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggc                                                                 <  2:1326059/140‑1 (MQ=255)                                                                     gcggTTCTTTCATGGTACTGGCGACCTCCAGCCACACTGCGGCCCGCATT‑GCGCCAGCGAGAAAAACCCCCTGACCGTCCAACACGCCGGTAAACCCGCCCCCCCGCGAAACGCCCAACCCCACCCCCAGaaaatcgcct                                                               >  2:1330942/1‑133 (MQ=255)                                                                       gtttCTTTCATGGTCCGGGCGCTTGTCAGCCACACGGCGG‑CAACATTGGCGCCGGCGAGAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggt                                                             <  2:1376816/138‑1 (MQ=255)                                                                       ggTTCTTTCATGGTACTGGCGACTGTCACCCACACTGCGGCCAGCATT‑GCGCCAGCGAAAAATCCCCCCTGACCGCCAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggt                                                             <  2:1210338/140‑1 (MQ=255)                                                                        gTTCTTTCATGGGACTGGCGACTGTCAGCCACACTGCGGGCAGCATT‑TCGCCCGCGAGAAATAACCCCCGACCGTCAAAAATGCCGGTAATCCAGCCGCCCAGCCAACCCCCAATCGCCACGCCAAGAAActggctgggg                                                            >  2:638210/1‑138 (MQ=255)                                                                         ttCTTTTATTGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATT‑GCGCCCGCGAGAAAAACCCCCCGACCCTCAAAAAAGCCGTTAATCCAACCGCCCAGCGAACCGCCCAACCCCACGCCCAAAAAATGGCtggggg                                                           >  2:1005569/1‑140 (MQ=255)                                                                                    tccTGGCGCCTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAccc                                                <  2:1237104/138‑1 (MQ=255)                                                                                       tGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATc                                             >  2:586988/1‑140 (MQ=255)                                                                                       tGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATc                                             >  2:387758/1‑140 (MQ=255)                                                                                       tGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATc                                             >  1:1139900/1‑140 (MQ=255)                                                                                       tGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATc                                             >  1:122648/1‑140 (MQ=255)                                                                                         gCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                           >  1:1237104/1‑140 (MQ=255)                                                                                         gCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                           >  2:125461/1‑140 (MQ=255)                                                                                           gACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                         <  1:586988/140‑1 (MQ=255)                                                                                           gACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                         <  1:391642/140‑1 (MQ=255)                                                                                            aCTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                        >  2:911858/1‑140 (MQ=255)                                                                                            aCTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                        >  2:811402/1‑140 (MQ=255)                                                                                            aCTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                        >  1:396034/1‑140 (MQ=255)                                                                                              tgtaaGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                      >  1:1484600/5‑140 (MQ=255)                                                                                                ccacccccacTGCGGCCAGCTTG‑GCGCCGGCGAAAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCt                                    <  2:122648/133‑1 (MQ=255)                                                                                                    ccACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGt                                <  1:911858/140‑1 (MQ=255)                                                                                                        acTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAcc                            >  1:793369/1‑140 (MQ=255)                                                                                                          tGCGGCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                           <  1:31126/123‑1 (MQ=255)                                                                                                          tGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                           >  2:31126/1‑123 (MQ=255)                                                                                                          tGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                          >  2:1308519/1‑140 (MQ=255)                                                                                                            cGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                        >  2:486383/1‑140 (MQ=255)                                                                                                            cGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                        >  1:807123/1‑140 (MQ=255)                                                                                                            cGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                        >  1:714269/1‑140 (MQ=255)                                                                                                            cGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                        >  1:1452212/1‑140 (MQ=255)                                                                                                              gCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                       <  2:411713/123‑1 (MQ=255)                                                                                                              gCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                       >  1:411713/1‑123 (MQ=255)                                                                                                              gCCAGCATT‑GCGCCAGCGAAAAATCCCCCCTGCCCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                      <  2:1263664/140‑1 (MQ=255)                                                                                                                 aGCATT‑GCGCCGGGGAAAAATCCCCCCTGACCGACAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   <  2:1484600/140‑1 (MQ=255)                                                                                                                 aGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   <  1:387758/140‑1 (MQ=255)                                                                                                                 aGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   >  2:379373/1‑140 (MQ=255)                                                                                                                      tcgcgcCAGTAAGAAATCCCCCCTGACCGTCAAACAGGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct               <  1:1330942/138‑1 (MQ=255)                                                                                                                        gcgcCGGCAAAAAATCCCCCCTGCCCGTCAAACATGCCGTAAATCCACCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  2:396034/140‑1 (MQ=255)                                                                                                                        gcgcCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  1:638210/140‑1 (MQ=255)                                                                                                                        gcgcCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  1:1308519/140‑1 (MQ=255)                                                                                                                        gcccaaGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  2:807123/135‑1 (MQ=255)                                                                                                                             aGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa        <  1:486383/140‑1 (MQ=255)                                                                                                                              gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt       >  2:815408/1‑140 (MQ=255)                                                                                                                                 agaAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAg    >  1:664650/1‑140 (MQ=255)                                                                                                                                  gaAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGt   <  1:379373/140‑1 (MQ=255)                                                                                                                                  aaaaTCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCCCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGt   <  1:842064/139‑1 (MQ=255)                                                                                                                                   aaaTACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTg  >  2:792594/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                        ACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTG  >  NC_000913/445425‑445693

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