Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	547,203	A→G	7.8%	H18R (CAC→CGC)	fusA →	elongation factor G 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	547,203	0	A	G	7.8%	182.8 / 4.0	64	H18R (CAC→CGC)	fusA	elongation factor G 1
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base G (2/3);  total (29/35)
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

AAGAAGCGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCTCCGCTGCC  >  NC_002947/547055‑547351                                                                                                                                                     |                                                                                                                                                     aagaagCGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCa                                                                                                                                                      >  4:259737/1‑149 (MQ=255)  agaagCGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCAc                                                                                                                                                     >  1:215606/1‑149 (MQ=255)     agCGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGtt                                                                                                                                                  <  5:11154/149‑1 (MQ=255)       cGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGa                                                                                                                                                <  4:176336/149‑1 (MQ=255)         tGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGAcg                                                                                                                                              >  8:64025/1‑149 (MQ=255)          gAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGtcaacgt                                                                                                                                             >  1:9505/1‑143 (MQ=255)            aGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCgg                                                                                                                                           <  6:82628/149‑1 (MQ=255)                    aCCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGTGGGCacagat                                                                                                                                   <  1:272874/149‑6 (MQ=255)                      cGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGTGGGCacagatac                                                                                                                                 <  2:9505/149‑8 (MQ=255)                      cGTATGGCTGAAGCCAACAAAGCGTGCTCGCACTACCGCTTCTAAATCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGGGGGCacagatac                                                                                                                                 <  1:141861/149‑8 (MQ=255)                          tGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACcaccac                                                                                                                              >  5:246343/1‑148 (MQ=255)                            gCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTAc                                                                                                                           >  7:88320/1‑149 (MQ=255)                             cTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGAAGCGGGCAAGACCACCACTAcc                                                                                                                          >  3:138242/1‑149 (MQ=255)                                aaGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAg                                                                                                                       >  8:333483/1‑149 (MQ=255)                                  gCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCg                                                                                                                     >  2:169643/1‑149 (MQ=255)                                     aacaaAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGaa                                                                                                                  >  3:290149/1‑148 (MQ=255)                                      acaaAGCGTTCTCGCACTACCGCTTCTATTTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATc                                                                                                                 >  1:215038/1‑149 (MQ=255)                                      acaaAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATc                                                                                                                 <  3:134601/149‑1 (MQ=255)                                         aaGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTg                                                                                                              <  4:79917/149‑1 (MQ=255)                                          aGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGt                                                                                                             <  8:88320/149‑1 (MQ=255)                                           gCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTgat                                                                                                            <  2:215038/149‑3 (MQ=255)                                           gCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGtt                                                                                                            <  4:169507/149‑1 (MQ=255)                                            cGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTc                                                                                                           >  6:326627/1‑149 (MQ=255)                                              ttCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTa                                                                                                         <  3:259737/149‑1 (MQ=255)                                                ctcGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGTCGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTaca                                                                                                       >  5:320919/1‑149 (MQ=255)                                                   gCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGA‑GGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAgg                                                                                                   >  6:123366/1‑149 (MQ=255)                                                     aCTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGaccac                                                                                                                                 <  2:258135/118‑1 (MQ=255)                                                     aCTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGaccac                                                                                                                                 >  1:258135/1‑118 (MQ=255)                                                     aCTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGt                                                                                                  <  3:253954/149‑1 (MQ=255)                                                      cTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTc                                                                                                 <  6:147211/149‑1 (MQ=255)                                                            cTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGt                                                                                                             <  6:314392/131‑1 (MQ=255)                                                            cTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGt                                                                                                             >  5:314392/1‑131 (MQ=255)                                                             ttCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGcc                                                                                           >  3:5818/1‑148 (MQ=255)                                                                   ttCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAAt                                                                                    <  1:169643/149‑1 (MQ=255)                                                                   ttCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAAt                                                                                    <  3:77787/149‑1 (MQ=255)                                                                      aaGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGaccacc                                                                                                                                <  3:230456/102‑1 (MQ=255)                                                                      aaGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGaccacc                                                                                                                                >  4:230456/1‑102 (MQ=255)                                                                         cATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCTCCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGa                                                                              <  2:215606/149‑1 (MQ=255)                                                                          aTCAATCTTTTTGCGAGGGCTTTTTGGCTCTGACTACAGCAATTAACCGCTTCCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAg                                                                             <  4:290149/149‑1 (MQ=255)                                                                           tCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATgg                                                                                  <  5:375747/143‑1 (MQ=255)                                                                           tCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACCTTGACGCGGGCAAGACCACAACTACCGACCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATgt                                                                                  >  6:375747/1‑142 (MQ=255)                                                                             aaTCTTTTTGCGA‑GGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGc                                                                         <  5:123366/149‑1 (MQ=255)                                                                                 tttttGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCAt                                                                       >  2:135150/1‑148 (MQ=255)                                                                                   tttGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGAc                                                                    >  3:321198/1‑149 (MQ=255)                                                                                     tGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACg                                                                   <  7:333483/148‑1 (MQ=255)                                                                                       cGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGgc                                                                 <  3:187628/148‑1 (MQ=255)                                                                                        gAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTTCATGACGgcgc                                                               <  4:138242/149‑1 (MQ=255)                                                                                          gggCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGcgccg                                                             >  3:277070/1‑149 (MQ=255)                                                                                           ggCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGcgccgc                                                            >  1:289686/1‑149 (MQ=255)                                                                                            gCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCcgcg                                                           <  6:320919/149‑1 (MQ=255)                                                                                                tATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGacc                                                        <  2:289686/148‑1 (MQ=255)                                                                                                          actacAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGa                                                          >  4:293970/1‑136 (MQ=255)                                                                                                          actacAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGa                                                          <  3:293970/136‑1 (MQ=255)                                                                                                              cAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGt                                          <  4:5818/148‑1 (MQ=255)                                                                                                                 cAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTgcag                                      >  3:330065/1‑149 (MQ=255)                                                                                                                       aCCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTgcaggagcag                                <  1:29653/149‑1 (MQ=255)                                                                                                                              ccGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCggg                         <  4:277070/149‑1 (MQ=255)                                                                                                                              ccGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCg                           <  7:64025/147‑1 (MQ=255)                                                                                                                              ccGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCg                           <  1:135150/147‑1 (MQ=255)                                                                                                                                    cATTGGTATCAGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCa                   >  8:320929/1‑149 (MQ=255)                                                                                                                                             tcttgCCCGCGTCAACGTGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCt          >  8:296081/4‑149 (MQ=25)                                                                                                                                               gtgCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCgggg                        >  6:24710/1‑133 (MQ=255)                                                                                                                                               gtgCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCgggg                        <  5:24710/133‑1 (MQ=255)                                                                                                                                                   ccACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCTCCGCTg    >  3:167548/1‑149 (MQ=255)                                                                                                                                                   ccACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCTCCGCTg    >  6:1999/1‑149 (MQ=255)                                                                                                                                                     aCGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCTCCGCTGcc  >  5:384023/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AAGAAGCGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCTCCGCTGCC  >  NC_002947/547055‑547351

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