Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,566,181	A→C	8.2%	Y62D (TAC→GAC) ‡	PP_3148 ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,566,181	0	A	C	8.2%	220.9 / 10.0	73	Y62D (TAC→GAC) ‡	PP_3148	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (38/29);  new base C (2/4);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

TGCCTTGCTTGTCGAAGGTGTCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATCCGGTGTTCTT  >  NC_002947/3566034‑3566327                                                                                                                                                    |                                                                                                                                                   tGCCTTGCTTGTCGAAGGTGTCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAg                                                                                                                                                   >  3:282143/1‑149 (MQ=255)       gcttgTCGAAGGTGTCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGc                                                                                                                                               <  8:53618/147‑1 (MQ=255)       gcttgTCGAAGGTGTCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCg                                                                                                                                             >  6:202765/1‑149 (MQ=255)            tCGAAGGTGTCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGa                                                                                                                                         <  2:52818/148‑1 (MQ=255)            tCGAAGGTGTCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGAc                                                                                                                                        >  6:119968/1‑149 (MQ=255)              gAAGGTGTCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCACCGGCGACTACGTCGACGATGAca                                                                                                                                       <  1:150499/148‑2 (MQ=255)                  gtgtCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCACCGGCGACTACGTCGACGATGAcatctat                                                                                                                                  <  1:176481/149‑7 (MQ=255)                   tgtCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTgc                                                                                                                                 <  5:152925/149‑1 (MQ=255)                     tCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTgca                                                                                                                                <  5:98406/148‑1 (MQ=255)                        tGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCACCGGCGACTACGTCGACGATGAcatct                                                                                                                                    >  1:54/1‑137 (MQ=255)                        tGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCACCGGCGACTACGTCGACGATGAcatct                                                                                                                                    <  2:54/141‑5 (MQ=255)                            tCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACg                                                                                                                        >  2:261282/1‑149 (MQ=255)                              aCGATCTCTGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGTGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCt                                                                                                                      <  6:170570/149‑1 (MQ=255)                                gATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACt                                                                                                                                       <  1:169701/130‑1 (MQ=255)                                gATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACt                                                                                                                                       >  2:169701/1‑130 (MQ=255)                                gATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCtt                                                                                                                     >  3:31477/1‑148 (MQ=255)                                  tCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTTGGATCGGGGCGGAAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCt                                                                                                                  >  8:189949/1‑149 (MQ=255)                                  tCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCt                                                                                                                  <  5:202765/149‑1 (MQ=255)                                  tCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCt                                                                                                                  >  5:149465/1‑149 (MQ=255)                                  tCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCt                                                                                                                  <  4:282143/149‑1 (MQ=255)                                   cgcgGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCt                                                                                                                  <  8:232711/148‑1 (MQ=255)                                    gcgGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTgg                                                                                                                <  7:144272/149‑1 (MQ=255)                                            cTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGCCTGTCTGCAGCAGCACGCTTTCTGGCAGACGCa                                                                                                        >  6:316745/1‑149 (MQ=255)                                                     ccACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGc                                                                                                >  3:20368/1‑148 (MQ=255)                                                     ccACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCt                                                                                               >  5:129870/1‑149 (MQ=255)                                                      cACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTgg                                                                                                                <  4:232549/131‑1 (MQ=255)                                                      cACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTgg                                                                                                                >  3:232549/1‑131 (MQ=255)                                                       aCGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCa                                                                                             >  1:156851/1‑149 (MQ=255)                                                         gggATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATg                                                                                           <  1:195095/149‑1 (MQ=255)                                                            atttGGTAAACGGCGGTGGTATAGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTTTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATggag                                                                                        <  7:290257/146‑1 (MQ=255)                                                            aTCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATggag                                                                                        >  1:19209/1‑149 (MQ=255)                                                               tGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATggagg                                                                                       >  8:329942/1‑147 (MQ=255)                                                               tGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATggagg                                                                                       <  7:329942/147‑1 (MQ=255)                                                                    aaCGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCACCGGCGACTACGTCGACGATGAcat                                                                                                                                      >  4:220033/1‑93 (MQ=37)                                                                    aaCGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCACCGGCGACTACGTCGACGATGAcat                                                                                                                                      <  3:220033/95‑3 (MQ=37)                                                                         cggTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTgcag                                                                                                                               <  2:75491/97‑1 (MQ=255)                                                                            tGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGgc                                                                        >  5:50109/1‑149 (MQ=255)                                                                                aTCGGGGCGGCAGACCATGGCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCaa                                                                    <  4:31477/149‑1 (MQ=255)                                                                                 tCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTg                                                                                 >  6:321742/1‑135 (MQ=255)                                                                                 tCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTg                                                                                 <  5:321742/135‑1 (MQ=255)                                                                                 tCGGGGCGGCAGACCATCTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGAAGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGAAGCATGCCTCGCTCATGGAGGAACTTGGAGGTAGGCGCaa                                                                    >  2:15058/1‑148 (MQ=255)                                                                                        ggcAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCcg                                                            <  2:19209/149‑1 (MQ=255)                                                                                         gcAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCcgc                                                           >  1:151058/1‑149 (MQ=255)                                                                                               cATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGa                                                     >  4:220374/1‑149 (MQ=255)                                                                                                 tgtcgatgtcgGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATg                                                   >  6:274225/1‑149 (MQ=255)                                                                                                    cgctgtcgGCGGCGTCGAGCAGGTGGTAATAGATGTCATCGTCGACGTAGTCGCCGGGGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCg                                                <  8:338277/146‑1 (MQ=255)                                                                                                     gatgtcgGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCgg                                               >  3:18299/1‑149 (MQ=255)                                                                                                     gatgtcgGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCgg                                               <  2:156851/149‑1 (MQ=255)                                                                                                        gtcgGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCgg                                                         >  7:300962/1‑136 (MQ=255)                                                                                                        gtcgGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCgg                                                         <  8:300962/136‑1 (MQ=255)                                                                                                           ggcggcGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGt                                         <  7:340770/149‑1 (MQ=255)                                                                                                               gcgTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGACGATGCCGGTCAAGTCGCt                                     <  7:189949/149‑1 (MQ=255)                                                                                                                 gTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGa                                   >  3:248931/1‑149 (MQ=255)                                                                                                                  tCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGAt                                  <  4:20368/149‑1 (MQ=255)                                                                                                                       cAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACg                              >  1:106509/1‑148 (MQ=255)                                                                                                                          gTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCAtt                          >  7:287267/1‑149 (MQ=255)                                                                                                                          gTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCAtt                          >  8:11443/1‑149 (MQ=255)                                                                                                                          gTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTAATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCAtt                          >  3:66104/1‑149 (MQ=255)                                                                                                                           ttGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTc                         >  4:307458/1‑149 (MQ=255)                                                                                                                           ttGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTc                         <  7:11443/149‑1 (MQ=255)                                                                                                                              tAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGAc                      >  2:272952/1‑149 (MQ=255)                                                                                                                                  aGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCg                  >  7:66181/1‑149 (MQ=255)                                                                                                                                   gATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCg                  <  5:274225/148‑1 (MQ=255)                                                                                                                                   gATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCg                  <  7:339851/148‑1 (MQ=255)                                                                                                                                    aTGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGt                >  3:284482/1‑149 (MQ=255)                                                                                                                                         aTCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATcc           >  6:116754/1‑149 (MQ=255)                                                                                                                                         aTCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATcc           >  6:266615/1‑149 (MQ=255)                                                                                                                                         aTCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGAAGAAGCATGACTCGCTCATTGAGGAACTTGTCGGTAGGAGCAATCTTGACGCAGGCGATGCCGGTCAAGTCGCTGATCCGGCAGTCGACTTCGGTGATcc           >  4:286619/1‑149 (MQ=255)                                                                                                                                            gtcgACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATCCGgt        <  8:66181/149‑1 (MQ=255)                                                                                                                                            gtcgACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATCCGgt        >  8:253374/1‑149 (MQ=255)                                                                                                                                            gtcgACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATCCGgt        <  5:116754/149‑1 (MQ=255)                                                                                                                                                aCGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATCCGGTGTTc    >  4:13358/1‑149 (MQ=255)                                                                                                                                                  gTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATCCGGTGTTCtt  >  4:44569/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TGCCTTGCTTGTCGAAGGTGTCGTGGATCACGATCGCGGTTGGCTCGATCGCCCACGGGATCTGGTAAACGGCGGTGGGATCGGGGCGGCAGACCATGTCGATGTCGGCGGCGTCGAGCAGGTTGTAATAGATGTCATCGTCGACGTAGTCGCCGGTGACTGTCTGCAGCAGCACGCTTTCTGGCAGACGCATGCCTCGCTCATGGAGGAACTTGGCGGTAGGCGCAATCTTGCCGCGGGCGATGCCGGTCAAGTCGCTGATCACGCATTCGACTTCGGTGATCCGGTGTTCTT  >  NC_002947/3566034‑3566327

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