Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	714,958	G→A	11.0%	intergenic (+24/+53)	ispH → / ← PP_0607	4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	714,958	0	G	A	11.0%	142.5 / 5.9	54	intergenic (+24/+53)	ispH/PP_0607	4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base A (3/3);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCGTGGCGTGATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTC  >  NC_002947/714811‑715106                                                                                                                                                    |                                                                                                                                                       ggTGCGTGGCGTGATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGACCATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAgg                                                                                                                                                       >  6:270901/1‑149 (MQ=255)      gtggcgtgATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGgctg                                                                                                                                                  <  6:110761/149‑1 (MQ=255)       tggcgtgATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGgctgc                                                                                                                                                 <  4:127803/149‑1 (MQ=255)             gATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTg                                                                                                                                            >  7:59196/1‑147 (MQ=255)               tCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGc                                                                                                                                        <  2:242209/149‑1 (MQ=255)                cGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGcc                                                                                                                                       >  8:104843/1‑149 (MQ=255)                 gAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGccc                                                                                                                                      <  8:59196/149‑1 (MQ=255)                   aCAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCtt                                                                                                                                    >  1:117756/1‑149 (MQ=255)                       cTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTgcgc                                                                                                                                >  1:119332/1‑149 (MQ=255)                        tCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTgcgc                                                                                                                                >  3:193267/1‑148 (MQ=255)                          aaGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTgcgc                                                                                                                                <  4:193267/146‑1 (MQ=255)                            ggCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGacac                                                                                                                           <  2:119332/149‑1 (MQ=255)                            ggCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGaca                                                                                                                            >  2:70017/1‑148 (MQ=255)                            ggCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAAGGGCTGC‑AAAGCAGCCCCTgcattggcg                                                                                                                           <  2:112808/149‑8 (MQ=255)                               cTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGacac                                                                                                                           >  4:146571/1‑146 (MQ=255)                                  gggCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGa                                                                                                                       <  3:30833/147‑1 (MQ=255)                                       ccACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGc                                                                                                                                           >  3:35236/1‑122 (MQ=255)                                       ccACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGc                                                                                                                                           <  4:35236/122‑1 (MQ=255)                                        cACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAGGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAgggg                                                                                                                                                     >  8:100062/1‑112 (MQ=255)                                                   aaGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACcg                                                                                                     <  2:89311/148‑1 (MQ=255)                                                      agCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACcgctgc                                                                                                 >  3:67982/1‑146 (MQ=255)                                                         tGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTag                                                                                              >  3:57057/1‑148 (MQ=255)                                                            aCGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTagcagc                                                                                          <  1:238516/149‑1 (MQ=255)                                                            aCGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTagcagc                                                                                          <  5:205502/149‑1 (MQ=255)                                                              gggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCg                                                                                         <  1:67307/149‑1 (MQ=255)                                                               ggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGt                                                                                       <  4:57057/149‑1 (MQ=255)                                                                  gcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCAt                                                                                    <  5:147215/149‑1 (MQ=255)                                                                  gcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCAt                                                                                    <  8:2071/149‑1 (MQ=255)                                                                    gAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAggg                                                                                                                                                      >  6:68961/1‑83 (MQ=255)                                                                    gAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAggg                                                                                                                                                      <  5:68961/83‑1 (MQ=255)                                                                    gAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAggg                                                                                                                                                      >  4:249278/1‑83 (MQ=255)                                                                    gAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAggg                                                                                                                                                      <  3:249278/83‑1 (MQ=255)                                                                     aaGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGAAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCt                                                                                                >  7:53178/1‑134 (MQ=255)                                                                      agagAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGc                                                                                 <  4:67982/148‑1 (MQ=255)                                                                      agagAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCt                                                                                <  7:104843/149‑1 (MQ=255)                                                                         gaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGCCACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAg                                                                             >  2:206074/1‑149 (MQ=255)                                                                               cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcaggagc                                                                           >  5:211172/1‑139 (MQ=39)                                                                               cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcaggagc                                                                           <  6:211172/145‑7 (MQ=39)                                                                                  cTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCg                                                                     >  1:156890/1‑148 (MQ=255)                                                                                     ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                  <  3:146571/148‑1 (MQ=255)                                                                                               aaaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                        >  4:74792/1‑148 (MQ=255)                                                                                                            ttCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATg                                          >  5:88722/1‑149 (MQ=255)                                                                                                               gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                       >  4:21930/1‑133 (MQ=255)                                                                                                               gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                       <  3:21930/133‑1 (MQ=255)                                                                                                                ctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAg                                      <  1:63482/149‑1 (MQ=255)                                                                                                                ctcgctGATCTGACGGGCAGTACCGCCGACAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAg                                      <  3:4910/149‑1 (MQ=255)                                                                                                                 tcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGc                                     >  8:217257/1‑149 (MQ=255)                                                                                                                        tCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCa                                       >  4:128864/1‑140 (MQ=255)                                                                                                                        tCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCa                                       <  3:128864/140‑1 (MQ=255)                                                                                                                              gggCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATc                        <  3:74792/149‑1 (MQ=255)                                                                                                                                gCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                      >  1:226800/1‑149 (MQ=255)                                                                                                                                         ccgccAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                       >  2:68789/1‑107 (MQ=255)                                                                                                                                         ccgccAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                       <  1:68789/107‑1 (MQ=255)                                                                                                                                           tctcgcgCAAGGGCTGC‑AAAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagcag            <  1:114268/142‑1 (MQ=255)                                                                                                                                              aaTTCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacc                                  <  6:35087/123‑1 (MQ=255)                                                                                                                                              aaTTCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacc                                  >  5:35087/1‑123 (MQ=255)                                                                                                                                               gcgCAAGGGCTGC‑AGAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAgcggcg       >  7:228188/3‑149 (MQ=255)                                                                                                                                               gcgCAAGGGCTGC‑AAAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagcag            >  2:114268/3‑144 (MQ=255)                                                                                                                                                 gCAAGGGCTGC‑AAAGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGgcgc      <  8:64877/148‑1 (MQ=255)                                                                                                                                                 gCAAGGGCTGC‑AAAGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGcgcg     >  7:64877/1‑149 (MQ=255)                                                                                                                                                    ggggCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTc  >  7:143172/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GGTGCGTGGCGTGATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTC  >  NC_002947/714811‑715106

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