Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,348,043	A→G	6.2%	N173D (AAT→GAT)	PP_4704 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,348,043	0	A	G	6.2%	185.8 / 2.3	64	N173D (AAT→GAT)	PP_4704	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base G (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

AAACTACCGGCCAGCGAAGCCAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTGAAGGTCACT  >  NC_002947/5347896‑5348189                                                                                                                                                    |                                                                                                                                                    aaaCTACCGGCCAGCGAAGCCAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGaa                                                                                                                                                    >  8:33741/1‑149 (MQ=255)         ggCCAGCGAAGCCAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGc                                                                                                                                            >  4:52982/1‑149 (MQ=255)             aGCGAAGCCAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGa                                                                                                                                         >  1:336879/1‑148 (MQ=255)                 aaGCCAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAAc                                                                                                                                    >  7:45634/1‑149 (MQ=255)                 aaGCCAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAAc                                                                                                                                    <  7:18335/149‑1 (MQ=255)                   gCCAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCt                                                                                                                                  <  7:33741/149‑1 (MQ=255)                    ccAAACCCTACGAGCCGCCCAGCGACCCGTTCCAGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCTGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTc                                                                                                                                 <  8:297295/149‑1 (MQ=255)                     cAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑a                                                                                                                               <  8:190831/149‑1 (MQ=255)                         cccGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCCC‑AGCAc                                                                                                                           >  4:2199/1‑149 (MQ=255)                            gAAGAGCCGCCCAGCGACCCGTTCTCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCt                                                                                                                        >  3:178956/1‑149 (MQ=255)                                 gccgccCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTAAGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGccc                                                                                                                    >  8:271121/1‑148 (MQ=255)                                   cgccCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCCC‑AGCACCCTGCCCAAc                                                                                                                 <  3:2199/149‑1 (MQ=255)                                                    ccGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGa                                                                                                 <  3:52982/148‑1 (MQ=255)                                                           gTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTtggtg                                                                                         <  4:178956/149‑1 (MQ=255)                                                            tATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTTCGGCCCAGCGAGCCTATCGAGTTCCTCTACGACAACATTAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTtggtgg                                                                                        <  8:73387/149‑1 (MQ=255)                                                            tATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTtggtgg                                                                                        <  6:116357/149‑1 (MQ=255)                                                                  cagcaTCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCg                                                                                  >  8:96996/1‑149 (MQ=255)                                                                   agcaTCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTCAAGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCg                                                                                  <  1:76105/149‑1 (MQ=255)                                                                      aTCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTcc                                                                              <  8:45634/149‑1 (MQ=255)                                                                          aGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTtggt                                                                                          <  7:109932/133‑1 (MQ=255)                                                                          aGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTtggt                                                                                          >  8:109932/1‑133 (MQ=255)                                                                           gCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATg                                                                         >  5:179876/1‑149 (MQ=255)                                                                           gCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATCCCGACATGACCTGGGTGGCCAGCGGTCCCaaag                                                                         >  8:286749/1‑147 (MQ=255)                                                                               aTCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCg                                                                      >  2:19112/1‑148 (MQ=255)                                                                                      gTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATc                                                              <  2:336879/149‑1 (MQ=255)                                                                                      gTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATc                                                              <  4:107118/149‑1 (MQ=255)                                                                                            tGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTtggtg                                                                                         >  2:30475/1‑116 (MQ=255)                                                                                            tGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTtggtg                                                                                         <  1:30475/116‑1 (MQ=255)                                                                                                 ttCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAgg                                                    >  3:178612/1‑148 (MQ=255)                                                                                                 ttCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAAACTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAgg                                                    >  2:47871/1‑148 (MQ=255)                                                                                                  tCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGa                                                  >  5:50903/1‑149 (MQ=255)                                                                                                    gAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGaa                                                 <  4:178612/148‑1 (MQ=255)                                                                                                    gAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAAc                                                <  6:102301/149‑1 (MQ=255)                                                                                                    gAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAAc                                                >  6:305787/1‑149 (MQ=255)                                                                                                    gAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAAc                                                >  3:284231/1‑149 (MQ=255)                                                                                                     aGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCACAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAAcc                                               >  8:215303/1‑149 (MQ=255)                                                                                                      ggACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACTACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAAccc                                              <  7:286749/149‑1 (MQ=255)                                                                                                      ggACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGGATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGAc                                                           >  7:202435/1‑136 (MQ=255)                                                                                                      ggACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGGATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGAc                                                           <  8:197103/136‑1 (MQ=255)                                                                                                      ggACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGGATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGAc                                                           <  8:202435/136‑1 (MQ=255)                                                                                                      ggACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGGATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGAc                                                           >  7:197103/1‑136 (MQ=255)                                                                                                         ccTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGa                                           <  7:215303/149‑1 (MQ=255)                                                                                                             cGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGc                                         <  1:12787/147‑1 (MQ=255)                                                                                                               gCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGt                                     >  2:235907/1‑149 (MQ=255)                                                                                                                        gCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCg                            >  1:153985/1‑149 (MQ=255)                                                                                                                          cGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATccc                           >  4:101856/1‑146 (MQ=255)                                                                                                                            aTCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTc                        <  7:271121/149‑1 (MQ=255)                                                                                                                                 gTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAg                   >  2:30660/1‑149 (MQ=255)                                                                                                                                  ttCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAgg                  <  6:179876/149‑1 (MQ=255)                                                                                                                                    ccTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTTGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGca                >  4:139793/1‑149 (MQ=255)                                                                                                                                     ctctACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGc                 <  3:125030/147‑1 (MQ=255)                                                                                                                                     ctctACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGc                 >  4:125030/1‑147 (MQ=255)                                                                                                                                      tctACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGcaca              >  4:168535/1‑149 (MQ=255)                                                                                                                                       ctACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAg                   <  4:300676/143‑1 (MQ=255)                                                                                                                                       ctACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAg                   >  3:300676/1‑143 (MQ=255)                                                                                                                                       ctACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCacac             <  7:96996/149‑1 (MQ=255)                                                                                                                                          cGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTg           >  6:292352/1‑148 (MQ=255)                                                                                                                                          cGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTg           >  6:93669/1‑148 (MQ=255)                                                                                                                                           gACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTg           >  4:275258/1‑147 (MQ=255)                                                                                                                                           gACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTGaa         <  1:19112/149‑1 (MQ=255)                                                                                                                                            acaacaTGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTGAAg        <  5:292352/149‑1 (MQ=255)                                                                                                                                            acaacaTGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACACTGAGCCTGTCGCCGATCGCCTc                        >  3:245749/1‑133 (MQ=255)                                                                                                                                            acaacaTGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACACTGAGCCTGTCGCCGATCGCCTc                        <  4:245749/133‑1 (MQ=255)                                                                                                                                                 aTGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTGAAGGTCAc   >  1:181682/1‑149 (MQ=255)                                                                                                                                                   gAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTGAAGGTCACt  >  4:263526/1‑148 (MQ=255)                                                                                                                                                    aaTCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAgg                  <  5:197323/131‑1 (MQ=255)                                                                                                                                                    aaTCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAgg                  >  6:197323/1‑131 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAACTACCGGCCAGCGAAGCCAAACCCGAAGAGCCGCCCAGCGACCCGTTCCCGCTGCGTATTGGCAGCATCAAGCTCATCGACGGTTACCTGCACTTCGAGGACCTGCGGCCCAGCGAGCCGATCGAGTTCCTCTACGACAACATGAATCTGGAGCTGAAAAACCTC‑AGCACCCTGCCCAACGACAATGCCGACATGACCTTGGTGGCCAGCGGTCCCAATGGCGGACGCATCGACTGGAAGGGAACCCTGAGCCTGTCGCCGATCGCCTCCGAAGGCACACTGAAGGTCACT  >  NC_002947/5347896‑5348189

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