Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,111,977	T→A	8.7%	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	8.7%	223.6 / 8.1	77	H118Q (CAT→CAA)	PP_5362	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base A (4/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGT  >  NC_002947/6111832‑6112121                                                                                                                                                    |                                                                                                                                                 cacaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CAATg                                                                                                                                                >  8:85547/1‑148 (MQ=255)  acaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCa                                                                                                                                              <  5:381335/149‑1 (MQ=255)      aCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCtt                                                                                                                                          >  3:240599/1‑149 (MQ=255)         aaCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CAATGCAGCTTGaa                                                                                                                                       <  7:85547/149‑2 (MQ=255)          aCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATc                                                                                                                                      >  5:277717/1‑149 (MQ=255)          aCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCACCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATc                                                                                                                                      >  3:79180/1‑149 (MQ=255)            aGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGa                                                                                                                                    <  4:240599/149‑1 (MQ=255)             gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAg                                                                                                                                   <  6:241507/149‑1 (MQ=255)              ccGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAgg                                                                                                                                  <  7:118220/149‑1 (MQ=255)               cGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAggg                                                                                                                                 <  7:332432/149‑1 (MQ=255)               cGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAggg                                                                                                                                 <  4:362757/149‑1 (MQ=255)               cGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAggg                                                                                                                                 <  5:266636/149‑1 (MQ=255)                      gCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGt                                                                                                                          >  8:21842/1‑149 (MQ=255)                          aTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGAc                                                                                                                      <  8:213757/149‑1 (MQ=255)                             gACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGTGGTGCCGATGACGcc                                                                                                                   >  3:312597/1‑149 (MQ=255)                               cTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctc                                                                                                                 <  2:13235/149‑1 (MQ=255)                                     ttACCGGCCTGGCCAGCGATGGTGAGGCGTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTc                                                                                                            <  4:79180/148‑1 (MQ=255)                                           gCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAt                                                                                                          >  3:133165/1‑144 (MQ=255)                                           gCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAt                                                                                                          <  4:133165/144‑1 (MQ=255)                                                gCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGCCGCCTCCTTTCATTTTGAACGtt                                                                                                >  4:221382/1‑149 (MQ=255)                                                 ccAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttt                                                                                               >  4:163900/1‑149 (MQ=255)                                                    gCGATGGTGAGGCCTGGGTGGGCTACGAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttttt                                                                                             >  1:4011/1‑149 (MQ=255)                                                     cGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtttt                                                                                                       >  4:222180/1‑137 (MQ=255)                                                     cGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtttt                                                                                                       <  3:222180/137‑1 (MQ=255)                                                           tGAGGCCTTGGTGGGCTAA‑AGCCGCATCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCATCGGGGAGCGCTTCAAGTTG‑AAGTGCAGCTTGATCGAGGGGTGCCGTTGCCGCCTCCTTTAAGTTTGAACGTTTTTTCTctgct                                                                                     >  8:334760/1‑149 (MQ=255)                                                           tGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGTTG‑CATTGCAGCTTGCTCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctgct                                                                                     >  8:329688/1‑149 (MQ=255)                                                               gCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGtt                                                                                                                         <  5:279072/109‑1 (MQ=255)                                                               gCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGtt                                                                                                                         >  6:279072/1‑109 (MQ=255)                                                                 cTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAt                                                                               >  2:79239/1‑149 (MQ=255)                                                                  tgggtgggCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtt                                                                                                         >  1:212078/1‑122 (MQ=255)                                                                  tgggtgggCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtt                                                                                                         <  2:212078/122‑1 (MQ=255)                                                                        ggCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTaa                                                                        <  6:277717/149‑1 (MQ=255)                                                                          cTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGAGGCCGCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGt                                                                      >  4:82505/1‑149 (MQ=255)                                                                          cTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAg                                                                       >  1:121087/1‑148 (MQ=255)                                                                           tAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGt                                                                      >  1:351056/1‑148 (MQ=255)                                                                               aGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGc                                                                  >  3:6527/1‑149 (MQ=255)                                                                                 ccgccgCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCtt                                                                >  4:23364/1‑149 (MQ=255)                                                                                 ccgccgCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCtt                                                                >  6:28/1‑149 (MQ=255)                                                                                  cgccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTc                                                               <  7:21842/149‑1 (MQ=255)                                                                                   gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCg                                                              <  2:121087/149‑1 (MQ=255)                                                                                    ccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCt                                                                                 >  4:329847/1‑129 (MQ=255)                                                                                    ccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCt                                                                                 <  3:329847/129‑1 (MQ=255)                                                                                    ccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                             <  3:23364/149‑1 (MQ=255)                                                                                      gcgTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGc                                                           <  3:327325/149‑1 (MQ=255)                                                                                        gTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCAc                                                         <  2:46934/149‑1 (MQ=255)                                                                                        gTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCAc                                                         <  3:163900/149‑1 (MQ=255)                                                                                         tAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTg                                                                                                      <  1:274735/104‑1 (MQ=255)                                                                                         tAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTg                                                                                                      >  2:274735/1‑103 (MQ=255)                                                                                            gCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACgcgc                                                     >  7:107801/1‑149 (MQ=255)                                                                                            gCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACgcgc                                                     >  2:156773/1‑149 (MQ=255)                                                                                                  ccaccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaaa                                               >  4:221464/1‑149 (MQ=255)                                                                                                   caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaaa                                               >  6:144077/1‑148 (MQ=255)                                                                                                   caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaaa                                               >  1:189127/1‑148 (MQ=255)                                                                                                           aCCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGcc                                       <  8:107801/148‑1 (MQ=255)                                                                                                           aCCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCt                                      <  2:189127/149‑1 (MQ=255)                                                                                                             cgcgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTg                                    <  5:28/149‑1 (MQ=255)                                                                                                              gcgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGa                                   <  3:221382/149‑1 (MQ=255)                                                                                                               cgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGAc                                  <  2:351056/149‑1 (MQ=255)                                                                                                               cgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGAc                                  <  1:156773/149‑1 (MQ=255)                                                                                                                   accAGTTCCAGCGGGTAGCGCTTCAAGCTG‑CATTGCAGCTTTATCTAGGGGGGCCGTTTACGCCACCTTTCATTTTGAACGTTTTTTCTCTCCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCaa                               <  1:79239/147‑1 (MQ=255)                                                                                                                           cAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGCTGCCGTTGACGCCTCCTTTCAATTTCAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCCCTGAAAACACCAGCCTTGACCAACCCTTTGt                       <  6:260754/148‑1 (MQ=255)                                                                                                                                   aGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                          <  8:296743/121‑1 (MQ=255)                                                                                                                                   aGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                          >  7:296743/1‑121 (MQ=255)                                                                                                                                        ttCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGaaaa            >  7:108833/1‑146 (MQ=255)                                                                                                                                        ttCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGaaaa            <  8:108833/146‑1 (MQ=255)                                                                                                                                        atCAAGCTG‑CAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCttt                         <  8:100822/132‑1 (MQ=255)                                                                                                                                        atCAAGCTG‑CAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCttt                         >  7:100822/2‑133 (MQ=255)                                                                                                                                         tCAAGCTG‑CAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCt         <  7:114771/148‑1 (MQ=255)                                                                                                                                         tCAAGCTG‑CAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCt         >  8:114771/1‑148 (MQ=255)                                                                                                                                          cAAGCTG‑CAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAATACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTAt       >  5:330474/1‑149 (MQ=255)                                                                                                                                           aaGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGAc                  >  3:184567/1‑137 (MQ=255)                                                                                                                                           aaGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGAc                  <  4:184567/137‑1 (MQ=255)                                                                                                                                           aaGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATg      <  2:4011/149‑1 (MQ=255)                                                                                                                                               tG‑CATTGCAGCTTGATCGAGGGGTCCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCg   >  4:170975/1‑148 (MQ=255)                                                                                                                                                g‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCa                                                                    <  5:148781/82‑1 (MQ=255)                                                                                                                                                g‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCa                                                                    >  6:148781/1‑82 (MQ=255)                                                                                                                                                g‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGt  >  4:19707/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTAC‑AGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTG‑CATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGT  >  NC_002947/6111832‑6112121

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