Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	564,393	T→C	7.1%	intergenic (+45/‑128)	rplQ → / → katA	50S ribosomal protein L17/catalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	564,393	0	T	C	7.1%	198.4 / 5.0	70	intergenic (+45/‑128)	rplQ/katA	50S ribosomal protein L17/catalase
Reads supporting (aligned to +/- strand):  ref base T (35/30);  new base C (2/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

TGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCA  >  NC_002947/564245‑564539                                                                                                                                                      |                                                                                                                                                   tGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGttttt                                                                                                                                                    >  4:203554/1‑149 (MQ=255)    gcATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGt                                                                                                                                                 <  2:983/149‑1 (MQ=255)     cATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCACTAG‑GCCCGGTTTTCTGta                                                                                                                                                <  4:24905/149‑2 (MQ=255)            aaGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACAt                                                                                                                                         <  3:17112/149‑1 (MQ=255)                  ggTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTg                                                                                                                                   >  4:256349/1‑149 (MQ=255)                   gTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGt                                                                                                                                  <  1:228812/149‑1 (MQ=255)                   gTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGt                                                                                                                                  <  1:233429/149‑1 (MQ=255)                    tttCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGtt                                                                                                                                 >  7:53087/1‑149 (MQ=255)                          cgcTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTc                                                                                                                           >  3:167228/1‑149 (MQ=255)                             tGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAttt                                                                                                                            >  6:89238/1‑145 (MQ=255)                             tGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAttt                                                                                                                            <  5:89238/145‑1 (MQ=255)                               gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCTTCCGGTCGGCGGTTCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAt                                                                                                                        <  7:139515/147‑1 (MQ=255)                               gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTTCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAtt                                                                                                                       <  5:24303/148‑1 (MQ=255)                               gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAtt                                                                                                                       <  6:232965/148‑1 (MQ=255)                                 gACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGAt                                                                                                                    <  5:152186/149‑1 (MQ=255)                                       gcgcCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGACGGCGGTGCTGGAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAg                                                                                                               >  1:221221/1‑148 (MQ=255)                                        cgcCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGt                                                                                                              >  7:229329/1‑148 (MQ=255)                                                gCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAAt                                                                                                      >  6:303412/1‑148 (MQ=255)                                                         gAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTcc                                                                                            <  5:303412/149‑1 (MQ=255)                                                         gAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTcc                                                                                            <  8:53087/149‑1 (MQ=255)                                                         gAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAG‑GCCCGGTTTTCTGta                                                                                                                                                >  3:59545/1‑96 (MQ=37)                                                         gAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAG‑GCCCGGTTTTCTGta                                                                                                                                                <  4:59545/97‑2 (MQ=37)                                                           gCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATGTATTTTTTATTTATTATTATGTTATTCATTTAATGAATTTGtacaa                                                                                          >  7:257747/1‑145 (MQ=255)                                                           gCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCt                                                                                          <  4:167228/149‑1 (MQ=255)                                                            gtGGTTGCTCGGCCGGTCGGCGGGGCTGTAGAATCTGCTGAGTAAGACGTTCAGTCTGCTACCGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTAATGGTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCt                                                                                          <  5:117872/147‑1 (MQ=255)                                                                  gATCGTCCGGTCGTCGGTGCTGTAGAAGCTGCTTAGTAAGACGTTCTGTCTGCTACAGAAAACCGTGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATc                                                                                    <  7:166560/148‑1 (MQ=255)                                                                   aTCGTCCGGTAGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCtt                                                                                  <  5:130507/149‑1 (MQ=255)                                                                           gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATc                                                                                    <  7:191370/139‑1 (MQ=255)                                                                           gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATc                                                                                    >  8:191370/1‑139 (MQ=255)                                                                              ggcggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATa                                                                       >  3:79460/1‑149 (MQ=255)                                                                                   tGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctt                                                                  <  8:221134/149‑2 (MQ=255)                                                                                      tGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctc                                                                >  4:64481/1‑148 (MQ=255)                                                                                       gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctct                                                               >  3:91372/1‑148 (MQ=255)                                                                                       gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctct                                                                 <  6:309730/146‑1 (MQ=255)                                                                                       gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctct                                                                 >  5:309730/1‑146 (MQ=255)                                                                                        tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTa                                                              >  7:172569/1‑148 (MQ=255)                                                                                         aGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTa                                                              >  2:108047/1‑147 (MQ=255)                                                                                          gAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAg                                                           >  2:235365/1‑149 (MQ=255)                                                                                              ctgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGa                                                       <  7:275937/149‑1 (MQ=255)                                                                                               tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGAt                                                      >  5:325224/1‑149 (MQ=255)                                                                                               tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGAt                                                      >  7:256262/1‑149 (MQ=255)                                                                                                gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGAt                                                      >  3:63291/1‑148 (MQ=255)                                                                                                gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACACTCTCTACAGCCGAt                                                      >  4:251117/1‑148 (MQ=255)                                                                                                      tAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCgg                                               >  4:98079/1‑149 (MQ=255)                                                                                                      tAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCgg                                               >  4:93552/1‑149 (MQ=255)                                                                                                            gTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCtt                                          >  8:205413/1‑148 (MQ=255)                                                                                                             ttCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGa                                        <  1:235365/149‑1 (MQ=255)                                                                                                               ctttctgCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGaa                                       <  3:226996/145‑1 (MQ=255)                                                                                                                  tctgCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGAcc                                   >  4:59797/1‑149 (MQ=255)                                                                                                                      cTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGaaa                               <  2:221221/149‑1 (MQ=255)                                                                                                                        aCAAAAAACCGGGCACTAG‑GCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg                             >  1:39429/1‑149 (MQ=255)                                                                                                                         cAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAAt                              <  7:205413/147‑1 (MQ=255)                                                                                                                          aaaaaaCCGGGCACTAG‑GCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAAt                                                                                                      <  3:164533/72‑1 (MQ=25)                                                                                                                          aaaaaaCCGGGCACTAG‑GCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAAt                                                                                                      >  4:164533/3‑74 (MQ=25)                                                                                                                          aGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg                             >  7:251663/1‑147 (MQ=255)                                                                                                                          aGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg                             >  5:138882/1‑147 (MQ=255)                                                                                                                           gAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgag                          <  4:63291/149‑1 (MQ=255)                                                                                                                           aaaaaCCGGGCACTAG‑GCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg                             >  6:199768/2‑146 (MQ=255)                                                                                                                            aaaaCCGGGCACTAG‑GCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTAATGTAATGAATTTGTCACTGTAATCTTGCTCGTCAATACTCTCTCTACAGAAGATTAGCCGGCCCTTGAAGACAGAAATGgaga                         >  1:128446/1‑149 (MQ=255)                                                                                                                            aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc                                                                                <  7:86617/94‑1 (MQ=255)                                                                                                                            aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc                                                                                >  8:86617/1‑94 (MQ=255)                                                                                                                            aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCTCTTGAAGACCGAAATGga                           >  6:16193/1‑147 (MQ=255)                                                                                                                            aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga                           >  7:251278/1‑147 (MQ=255)                                                                                                                            aaaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga                           >  4:164586/1‑147 (MQ=255)                                                                                                                                 cGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga                           <  3:164988/142‑1 (MQ=255)                                                                                                                                 cGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga                           >  4:164988/1‑142 (MQ=255)                                                                                                                                  gggC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCAt                   <  4:91372/149‑1 (MQ=255)                                                                                                                                    gC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGa                 <  6:325224/149‑1 (MQ=255)                                                                                                                                             cccGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGAtt         >  6:259216/1‑149 (MQ=255)                                                                                                                                                gtttttttGTTTATACATCTATTTTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATAATCTCTCTTCAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATtctc      <  3:251117/147‑1 (MQ=255)                                                                                                                                                ggTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATtctc      <  2:39429/149‑1 (MQ=255)                                                                                                                                                ggTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATtctc      <  5:199768/149‑1 (MQ=255)                                                                                                                                                   tttttGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCTCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTcac    <  5:16193/148‑1 (MQ=255)                                                                                                                                                    ttttGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTcacca  >  2:36766/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCA  >  NC_002947/564245‑564539

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