Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,855,829	T→G	6.7%	N64T (AAC→ACC)	folA ←	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,855,829	0	T	G	6.7%	230.7 / 7.0	75	N64T (AAC→ACC)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base G (2/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GTACAGCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCACGCGGTTCTC  >  NC_002947/5855684‑5855977                                                                                                                                                  |                                                                                                                                                     gTACAGCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGCCCACGATGtt                                                                                                                                                     >  5:155951/1‑147 (MQ=255)    cAGCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGTCTGGCGGCTGACCACGATGTTCAgcc                                                                                                                                                <  8:325731/149‑1 (MQ=255)    cAGCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAgcc                                                                                                                                                <  1:200069/149‑1 (MQ=255)    cAGCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAgcc                                                                                                                                                <  7:359498/149‑1 (MQ=255)      gCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAgccgc                                                                                                                                              >  4:268282/1‑149 (MQ=255)      gCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAgccgc                                                                                                                                              <  8:218341/149‑1 (MQ=255)      gCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGAACATCGTGGTCAgccgc                                                                                                                                              <  8:356523/149‑1 (MQ=255)        tGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGACTGACCACGATGTTCAgccgcc                                                                                                                                             >  7:102028/1‑148 (MQ=255)           gcgcCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAgg                                                                                                                                         <  1:129021/149‑1 (MQ=255)           gcgcCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAgg                                                                                                                                         <  6:166913/149‑1 (MQ=255)           gcgcCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAgg                                                                                                                                         <  4:334407/149‑1 (MQ=255)               ccccAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGc                                                                                                                                     >  3:164607/1‑149 (MQ=255)                 ccAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGAACATCGTGGTCAgccgc                                                                                                                                              >  5:248197/1‑138 (MQ=255)                 ccAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGAACATCGTGGTCAgccgc                                                                                                                                              <  6:248197/138‑1 (MQ=255)                     tcagcatcagcTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCACGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGgccgc                                                                                                                               >  3:68556/1‑149 (MQ=255)                         catcagcTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCCGCCGCCCAGGCAGCGGCCGccc                                                                                                                             >  6:212781/1‑147 (MQ=255)                          atcagcTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAg                                                                                                                          >  4:226623/1‑149 (MQ=255)                              gcTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAg                                                                                                                      <  6:126754/149‑1 (MQ=255)                              gcTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAg                                                                                                                      >  4:342460/1‑149 (MQ=255)                                  aCCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTcc                                                                                                                   <  3:226623/148‑1 (MQ=255)                                    cGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTCGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAg                                                                                                                >  2:302896/1‑149 (MQ=255)                                         cccTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGa                                                                                                                       <  8:243099/137‑1 (MQ=255)                                         cccTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGa                                                                                                                       >  7:243099/1‑137 (MQ=255)                                              cTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGcc                                                                                                      >  8:352288/1‑149 (MQ=255)                                                    cgcCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCatga                                                                                                <  4:22421/149‑1 (MQ=255)                                                    cgcCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCatga                                                                                                <  7:37787/149‑1 (MQ=255)                                                       ccACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCatgatg                                                                                              >  7:324416/1‑148 (MQ=255)                                                       ccACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCAtgatga                                                                                             >  5:216751/1‑149 (MQ=255)                                                        cACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATgatgat                                                                                            >  2:77802/1‑149 (MQ=255)                                                             cTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCt                                                                                       >  6:2541/1‑149 (MQ=255)                                                                  cccGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGcc                                                                                   >  1:110590/1‑148 (MQ=255)                                                                   ccGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTCGCGCCCCATGAGGATCGGCTTGcc                                                                                   >  4:77607/1‑147 (MQ=255)                                                                     gCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAgg                                                                               <  2:306659/149‑1 (MQ=255)                                                                       aCCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGggt                                                                             >  1:172737/1‑149 (MQ=255)                                                                         cAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTgc                                                                                                         <  5:4338/119‑1 (MQ=255)                                                                         cAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTgc                                                                                                         >  6:4338/1‑119 (MQ=255)                                                                         cAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGggtgg                                                                           >  7:334222/1‑149 (MQ=255)                                                                           aCGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGgtggtg                                                                         <  3:80493/149‑1 (MQ=255)                                                                           aCGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGgtggtg                                                                         <  2:110590/149‑1 (MQ=255)                                                                              ccTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGcc                                                                      <  5:2541/149‑1 (MQ=255)                                                                                tctTAAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACCATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGcc                                                                      >  4:353532/1‑147 (MQ=255)                                                                                  ttCAAGCGAGGCGAACACTTCGGCACCTTCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGACGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCtt                                                                    <  4:383602/147‑1 (MQ=255)                                                                                  ttCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCtt                                                                    <  1:302896/147‑1 (MQ=255)                                                                                  ttCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTg                                                                   >  1:82113/1‑148 (MQ=255)                                                                                   acAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCACGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGCCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCACGGTGGTGGCCTTg                                                                   >  1:340695/2‑147 (MQ=255)                                                                                    cAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGaa                                                                 <  6:216751/148‑1 (MQ=255)                                                                                    cAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGaa                                                                 <  4:164607/148‑1 (MQ=255)                                                                                         gAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGGCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACtt                                                           <  8:334222/149‑1 (MQ=255)                                                                                         gAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACtt                                                           <  3:268282/149‑1 (MQ=255)                                                                                         gAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACtt                                                           <  3:145998/149‑1 (MQ=255)                                                                                         gAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACtt                                                           <  7:352288/149‑1 (MQ=255)                                                                                            gCGAACACTTCTGCACCTGCAAGTTCAAGCCCGGGCTGGCGGCTGACCACTATGTTCAGCCGCCCAGGCATCGGCCGCCCCAGCGAGTCCCAGGTCATGCGCCCCATGATGTTCGGCTTGCCCAGGGTGGTGGCCTTGAATTACTTGaa                                                        <  2:131591/149‑1 (MQ=255)                                                                                             cGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAg                                                       <  1:77802/149‑1 (MQ=255)                                                                                              gAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACtt                                                           <  8:50063/144‑1 (MQ=255)                                                                                              gAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACtt                                                           >  7:50063/1‑144 (MQ=255)                                                                                                 cacTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTccc                                                   <  5:267871/149‑1 (MQ=255)                                                                                                    ttCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCg                                                >  8:173715/1‑149 (MQ=255)                                                                                                        gCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTg                                                                   >  7:172044/1‑126 (MQ=255)                                                                                                        gCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTg                                                                   <  8:172044/126‑1 (MQ=255)                                                                                                              gCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGc                                        >  1:146382/1‑147 (MQ=255)                                                                                                                aaGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAgg                                    <  2:172737/149‑1 (MQ=255)                                                                                                                  gTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAggg                                   >  4:112498/1‑148 (MQ=255)                                                                                                                   ttCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCa                                 >  7:126090/1‑149 (MQ=255)                                                                                                                   ttCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCa                                 >  8:197587/1‑149 (MQ=255)                                                                                                                       aGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGa                             <  8:5308/149‑1 (MQ=255)                                                                                                                             ggCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTc                       <  1:273050/149‑1 (MQ=255)                                                                                                                                 ggcggcTGAACATCGTGGTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGAt                    >  2:359637/1‑148 (MQ=255)                                                                                                                                 ggcggcTGAACATCGTGGTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGAt                    <  1:359637/148‑1 (MQ=255)                                                                                                                                  gcggcTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGc                  <  3:112498/149‑1 (MQ=255)                                                                                                                                   cggcTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGcc                 <  2:324820/149‑1 (MQ=255)                                                                                                                                           cACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCAcgc         >  3:42771/1‑149 (MQ=255)                                                                                                                                           cACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCAcgc         >  3:76128/1‑149 (MQ=255)                                                                                                                                           cACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGGGGCCTTGAAGTACTTCAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCAcgc         <  8:324416/149‑1 (MQ=255)                                                                                                                                            aCGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCAcgcg        >  3:229258/1‑149 (MQ=255)                                                                                                                                            aCGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCAcgcg        <  5:212781/149‑1 (MQ=255)                                                                                                                                                 gTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCACGCGGTtct   >  3:72466/1‑149 (MQ=255)                                                                                                                                                  ttCAGCCGCGCAGGCAGCGGCCGCCCCAGCGTGTGCCTGGTCTTGCGGCCCATGTTGATCGGGTTGCCGAGGGTGGTGGCCTTGAAGTACTTGAAGTGGCGCGTCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCACGCGGTtctc  <  2:340695/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTACAGCTGCGCGCCCCCAATCAGCATCAGCTCACCGACGCCCTGCTCGCGCGCCCACTGCTCGGCCCGCACCAACGCCTCTTCAAGCGAGGCGAACACTTCGGCACCTGCAAGTTCAAGCCCCGGCTGGCGGCTGACCACGATGTTCAGCCGCCCAGGCAGCGGCCGCCCCAGCGAGTCCCAGGTCTTGCGCCCCATGATGATCGGCTTGCCCAGGGTGGTGGCCTTGAAGTACTTGAAGTCCCCCGGCAAATGCCAGGGCATGGAGTTGTCGATGCCGATCACGCGGTTCTC  >  NC_002947/5855684‑5855977

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