Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,111,977	T→A	7.1%	H118Q (CAT→CAA)	PP_5362 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,111,977	0	T	A	7.1%	191.9 / 2.8	66	H118Q (CAT→CAA)	PP_5362	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base A (2/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAACC  >  NC_002947/6111831‑6112125                                                                                                                                                   |                                                                                                                                                      gCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTg                                                                                                                                                     <  1:181559/149‑1 (MQ=255) gCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTg                                                                                                                                                     <  2:112982/149‑1 (MQ=255)  cacaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTg                                                                                                                                                     >  2:255589/1‑148 (MQ=255)  cacaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGc                                                                                                                                                    <  7:1796/149‑1 (MQ=255)   acaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCa                                                                                                                                                   <  3:300927/149‑1 (MQ=255)     aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGCCCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGc                                                                                                                                                 >  8:21472/1‑149 (MQ=255)     aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGc                                                                                                                                                 >  6:44803/1‑149 (MQ=255)     aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGATCCAGCGGGGAGCGCTTCAAGATGCATTGCAGc                                                                                                                                                 >  4:7719/1‑149 (MQ=255)        ccAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTg                                                                                                                                              <  2:199079/149‑1 (MQ=255)          aaCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGAt                                                                                                                                            >  3:207518/1‑149 (MQ=255)          aaCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                            <  8:39328/149‑2 (MQ=255)                cGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAggg                                                                                                                                      <  3:305513/149‑1 (MQ=255)                  cTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGt                                                                                                                                    >  8:66170/1‑149 (MQ=255)                  cTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGt                                                                                                                                    <  8:334724/149‑1 (MQ=255)                          caTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑a                                                                                                                           <  3:185682/149‑1 (MQ=255)                           aTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑Ac                                                                                                                          <  1:255589/149‑1 (MQ=255)                                cTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGcctc                                                                                                                     <  4:207518/149‑1 (MQ=255)                                 tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGcctcc                                                                                                                    >  7:49616/1‑149 (MQ=255)                                 tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGcc                                                                                                                       >  7:197932/1‑146 (MQ=255)                                 tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGcc                                                                                                                       <  8:197932/146‑1 (MQ=255)                                   tttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGcctcct                                                                                                                   <  2:129384/148‑1 (MQ=255)                                            gCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTg                                                                                                          >  3:341049/1‑148 (MQ=255)                                               tGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                            >  2:370114/1‑111 (MQ=255)                                               tGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                            <  1:370114/112‑2 (MQ=255)                                                  ccAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGttt                                                                                                   <  5:172068/149‑1 (MQ=255)                                                      cGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTtc                                                                                               >  3:77356/1‑149 (MQ=255)                                                      cGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTtc                                                                                               <  8:49616/149‑1 (MQ=255)                                                         tGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCACTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTctct                                                                                            >  1:203245/1‑149 (MQ=255)                                                          ggTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTctct                                                                                            <  4:77356/148‑1 (MQ=255)                                                                  cTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAt                                                                                   >  5:137309/1‑149 (MQ=255)                                                                     ggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAtttt                                                                                <  3:203526/149‑1 (MQ=255)                                                                     ggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAtttt                                                                                <  2:66858/149‑1 (MQ=255)                                                                              cAGCCGCCGCGTAGGCTGGACCACCTTGAACCTCGTATCCAGTTCCAGCGTGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTc                                                                         <  3:315133/147‑1 (MQ=255)                                                                              cAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGt                                                                          >  8:229843/1‑146 (MQ=255)                                                                               aGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑CGGGCTCCTTTCAGTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGc                                                                      >  6:95312/1‑149 (MQ=255)                                                                                gccgccGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCtt                                                                                                                  <  7:64068/104‑1 (MQ=255)                                                                                gccgccGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCtt                                                                                                                  >  8:64068/1‑104 (MQ=255)                                                                                 ccgccgCGTAGGCTGGACCACCTTGAACCACGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTtc                                                                                               <  1:10389/122‑1 (MQ=255)                                                                                 ccgccgCGTAGGCTGGACCACCTTGAACCACGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTtc                                                                                               >  2:10389/1‑122 (MQ=255)                                                                                   gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTc                                                                   >  4:194184/1‑148 (MQ=255)                                                                                   gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCACTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCg                                                                  <  2:203245/149‑1 (MQ=255)                                                                                      gcgTAGGCTGGCCCCCCTTGCACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCGTGCTCGAGGGGTGCCGTTG‑ACGCCTCCTTTGATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCgcgc                                                               >  2:219501/1‑149 (MQ=255)                                                                                       cgTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCa                                                              <  5:44803/149‑1 (MQ=255)                                                                                        gTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCAc                                                             >  8:82653/1‑149 (MQ=255)                                                                                         tAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCAcg                                                            <  7:66170/149‑1 (MQ=255)                                                                                              tGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGtttt                                                                                                  >  6:13071/1‑106 (MQ=255)                                                                                                  ccaccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaa                                                    >  2:298048/1‑148 (MQ=255)                                                                                                   caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaaa                                                   >  4:351464/1‑148 (MQ=255)                                                                                                   caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTAGATCGAGGGGTGCCGTTG‑ACGcct                                                                                                                      >  7:299742/1‑81 (MQ=255)                                                                                                      cTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACAcc                                               >  6:337437/1‑149 (MQ=255)                                                                                                         gAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAg                                             <  3:78002/148‑1 (MQ=255)                                                                                                            ccGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCtt                                         >  3:37771/1‑149 (MQ=255)                                                                                                                 tATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGAcc                                     >  5:235498/1‑148 (MQ=255)                                                                                                                  aTCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCaa                                   >  8:173112/1‑149 (MQ=255)                                                                                                                   tCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCaa                                   <  6:137309/148‑1 (MQ=255)                                                                                                                     cAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAAccc                                <  7:173112/149‑1 (MQ=255)                                                                                                                     cAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAAcc                                 <  7:229843/148‑1 (MQ=255)                                                                                                                      aGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCt                               >  7:62943/1‑149 (MQ=255)                                                                                                                        ttCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGAAGGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCtt                              >  1:285999/1‑149 (MQ=255)                                                                                                                        ttCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCttt                             <  1:298048/149‑1 (MQ=255)                                                                                                                           cAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGt                           >  2:164673/1‑148 (MQ=255)                                                                                                                                gggAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTATCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGAcc                     <  4:341049/149‑1 (MQ=255)                                                                                                                                  ccctcgatCAAGCTGCAATGCAGCTGGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCg                    >  5:255120/8‑148 (MQ=255)                                                                                                                                            aGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATATTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGc         >  1:200040/1‑149 (MQ=255)                                                                                                                                             gCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGcc        >  1:142452/1‑149 (MQ=255)                                                                                                                                              cTGCATTGCAGCTTGATCGAGGGGTTCCGTTG‑ACTCCTCCTTTAATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCg       <  1:108033/149‑1 (MQ=255)                                                                                                                                                   ttGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAAcc  >  3:116325/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTG‑ACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAACC  >  NC_002947/6111831‑6112125

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