Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	159,376	T→A	5.9%	Q338L (CAA→CTA)	PP_0149 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	159,376	0	T	A	5.9%	209.4 / 2.9	68	Q338L (CAA→CTA)	PP_0149	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/29);  new base A (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

TCTCGATTTCGTCGCGGTTGAGCTTCAAGGGCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGATTGGCTT  >  NC_002947/159230‑159523                                                                                                                                                   |                                                                                                                                                         tctcGATTTCGTCGCGGTTGAGCTTCAAGGGCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGc                                                                                                                                                        <  3:207942/149‑1 (MQ=255)        ttCGTCGCGGTTGAGCTTCAAGGGCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGTCAGTTGCAcgtcg                                                                                                                                                  <  4:255153/148‑1 (MQ=255)            tcgCGGTTGAGCTTCAAGGGCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGAc                                                                                                                                             >  8:158769/1‑149 (MQ=255)                        ttCAAGGGCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTa                                                                                                                                 >  4:289766/1‑149 (MQ=255)                           aaGGGCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGa                                                                                                                              >  1:48214/1‑149 (MQ=255)                           aaGGGCACCGAGGCCAGGCGGCGCAGCGCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGa                                                                                                                              <  3:137883/149‑1 (MQ=255)                             gggCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGa                                                                                                                              >  3:116430/1‑147 (MQ=255)                             gggCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGa                                                                                                                              >  6:255854/1‑147 (MQ=255)                                  ccGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCAcc                                                                                                                       <  5:255854/149‑1 (MQ=255)                                       gCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGACTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCt                                                                                                                  >  6:271686/1‑149 (MQ=255)                                         cAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCttt                                                                                                                >  3:126802/1‑149 (MQ=255)                                         cAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCtt                                                                                                                 <  4:126802/148‑1 (MQ=255)                                          aGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCTTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACCGCACCAGGCttt                                                                                                                >  8:113176/1‑148 (MQ=255)                                            gCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGAt                                                                                                             >  3:307047/1‑149 (MQ=255)                                                cgcAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAg                                                                                                         <  5:71326/149‑1 (MQ=255)                                                     cTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTc                                                                                                      >  1:173647/1‑147 (MQ=255)                                                            ttGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAggcgg                                                                                               >  1:55131/1‑147 (MQ=255)                                                              ggTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCa                                                                                             <  6:279494/147‑1 (MQ=255)                                                               gtgtATTCGTCGCTCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGAAAGTTGCACGTCGTTGACGAAGGCGCCGTACGACAGAACCAGGCTTTGATCGAGGTCCAGGCGGCAGTc                                                                                          >  3:191943/1‑149 (MQ=255)                                                               gtgtATTCGTCGATCTCCTTGCCCAGCGTCAGCACGAATAGCCTCGTCGGGTAGTTGCAGACAAGGCCGTCCCAGGTGGAGTTTTGTACGTCGTTGAGGAAGGCGGGGTTCGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTc                                                                                          <  2:210413/149‑1 (MQ=255)                                                                  tATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑g                                                                                       <  8:152150/148‑1 (MQ=255)                                                                   aTTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCCCCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGt                                                                                     >  6:190637/1‑149 (MQ=255)                                                                    ttCCTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑c                                                                                 >  5:79534/1‑149 (MQ=255)                                                                       gtcgATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑g                                                                                       >  4:276898/1‑143 (MQ=255)                                                                       gtcgATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑g                                                                                       <  3:276898/143‑1 (MQ=255)                                                                       gtcgATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑Ccgc                                                                              >  2:165778/1‑149 (MQ=255)                                                                       gtcgATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑Ccg                                                                               <  7:158769/148‑1 (MQ=255)                                                                          gATCCCCGGGCCAAGCGGCAGCACGAATAGCCTCGACGGGTAGTGTCCTACAAGGCCGTCCCAGCTGGACAGGTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCATGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTg                                                                           <  3:181707/149‑1 (MQ=255)                                                                                cTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGggcggc                                                                     >  3:181971/1‑149 (MQ=255)                                                                                cTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTAGACGAAGGCGCGGTACGACAGCCCCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGggcggc                                                                     >  2:178298/1‑149 (MQ=255)                                                                                   gCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCa                                                                                                    >  4:72557/1‑119 (MQ=255)                                                                                   gCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCa                                                                                                    <  3:72557/119‑1 (MQ=255)                                                                                        gCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTg                                                             >  3:93584/1‑149 (MQ=255)                                                                                        gCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTg                                                             >  3:98039/1‑149 (MQ=255)                                                                                          ggCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCa                                                           <  4:116430/149‑1 (MQ=255)                                                                                              gcaCGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGCCAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATc                                                       >  7:205137/1‑149 (MQ=255)                                                                                                aCGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑c                                                    >  4:296591/1‑149 (MQ=255)                                                                                                 cGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑cc                                                   <  4:181971/149‑1 (MQ=255)                                                                                                  gAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCtt                                                                                                                 <  3:144444/91‑1 (MQ=255)                                                                                                  gAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑c                                                    >  2:94156/1‑147 (MQ=255)                                                                                                  gAATAGCCTCGACGGGTACTTGCCGACAAGGCCGTCCCAGCTCGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCtt                                                                                                                 >  4:144444/1‑91 (MQ=255)                                                                                                   aaTAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAg                                                 <  2:173647/149‑1 (MQ=255)                                                                                                       gCCTCGTCGGGTAGTTGCCGTCAAGGGCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGgcg                                              <  8:63720/148‑1 (MQ=255)                                                                                                       gCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGgcg                                              >  2:61926/1‑148 (MQ=255)                                                                                                        ccTCGACGGGTAGTTGCAGACAAGGCCGTCCCAGCTGGACAGTTGCACGTAGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCACGGCTTGCAGATCA‑ACACGCgca                                            >  2:281162/1‑149 (MQ=255)                                                                                                         cTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCgcag                                           >  2:222602/1‑149 (MQ=255)                                                                                                         cTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCgcag                                           >  3:127330/1‑149 (MQ=255)                                                                                                            gACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCgcagc                                          <  1:94156/147‑1 (MQ=255)                                                                                                            gACGGGTAGTTGCCCACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCgcagc                                          <  3:217015/147‑1 (MQ=255)                                                                                                               gggTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCt                                     >  1:163483/1‑149 (MQ=255)                                                                                                                 gTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCtgt                                   >  2:166217/1‑149 (MQ=255)                                                                                                                 gTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTAGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCtgt                                   >  4:283537/1‑149 (MQ=255)                                                                                                                 gTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTAGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCtgt                                   >  4:278877/1‑149 (MQ=255)                                                                                                                       gCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGcccc                              <  6:79534/148‑1 (MQ=255)                                                                                                                       gCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCACAGCATGCTGCGcc                                <  3:289766/146‑1 (MQ=255)                                                                                                                         cgtgaaGGCCGTCGCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGcc                                <  6:28341/140‑1 (MQ=255)                                                                                                                         cGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCTGGGTGAGCCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGccc                               >  1:265240/1‑149 (MQ=255)                                                                                                                         cGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGcc                                >  5:28341/1‑144 (MQ=255)                                                                                                                              aGGCCGTCCCTGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCt                      <  5:36046/149‑1 (MQ=255)                                                                                                                              aGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCt                      <  5:271686/149‑1 (MQ=255)                                                                                                                                      ccAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGt              >  2:182736/1‑149 (MQ=255)                                                                                                                                        aGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTc            <  1:166217/149‑1 (MQ=255)                                                                                                                                             gACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCACCCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGa         <  1:61926/148‑1 (MQ=255)                                                                                                                                             gACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGAtt       <  4:93584/149‑1 (MQ=255)                                                                                                                                             gACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGAtt       <  4:98039/149‑1 (MQ=255)                                                                                                                                               cAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGATTgg     >  8:103697/1‑149 (MQ=255)                                                                                                                                                aGTAGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGATTGGc    <  3:283537/149‑1 (MQ=255)                                                                                                                                                aGTAGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGATTGGc    <  3:278877/149‑1 (MQ=255)                                                                                                                                                  ttGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGATTGGCtt  >  5:286319/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                         TCTCGATTTCGTCGCGGTTGAGCTTCAAGGGCACCGAGGCCAGGCTGCGCAGCTCGGTCTTGGTGTATTCGTCGATCACCTGGCCCAGCGGCAGCACGAATAGCCTCGACGGGTAGTTGCCGACAAGGCCGTCCCAGCTGGACAGTTGCACGTCGTTGACGAAGGCGCGGTACGACAGCACCAGGCTTTGATCGAGGTCCAGGCGGCAGTCT‑GGT‑‑‑CCGCGTGGGCGGCCCGGCTTGCAGATCA‑CCAGGCGCAGCATGCTGTGCCCCCAGCGGCTGACCAGGTTCTGATTGGCTT  >  NC_002947/159230‑159523

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