Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,239,055	G→T	6.1%	Q62H (CAG→CAT)	PP_2837 →	uronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,239,055	0	G	T	6.1%	211.9 / 3.0	66	Q62H (CAG→CAT)	PP_2837	uronate transporter
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base T (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TTCCTGCTGCTGATCGGCATCGTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGCTGGGGATTGGC  >  NC_002947/3238909‑3239202                                                                                                                                                   |                                                                                                                                                    ttCCTGCTGCTGATCGGCATCGTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGtt                                                                                                                                                   <  2:230720/149‑1 (MQ=255)       ctgctgATCGGCATCGTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCaa                                                                                                                                              <  1:119112/148‑1 (MQ=255)              tCGGCATCGTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATg                                                                                                                                      <  2:281571/149‑1 (MQ=255)                 gCATCGTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGatc                                                                                                                                   <  8:83727/149‑1 (MQ=255)                    tCGTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGatcatc                                                                                                                                <  5:152577/149‑1 (MQ=255)                    tCGTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGatcatc                                                                                                                                >  8:44146/1‑149 (MQ=255)                      gTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGa                                                                                                                              >  7:241539/1‑149 (MQ=255)                         aaCTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCg                                                                                                                           <  3:294343/149‑1 (MQ=255)                          aCTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCATTGGCAACTGGGCaaag                                                                                                                                     >  5:258866/1‑135 (MQ=255)                          aCTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCATTGGCAACTGGGCaaag                                                                                                                                     <  6:258866/138‑4 (MQ=255)                               cTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTggg                                                                                                                     >  5:112165/1‑149 (MQ=255)                               cTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTggg                                                                                                                     >  3:190223/1‑149 (MQ=255)                               cTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTggg                                                                                                                     >  8:272887/1‑149 (MQ=255)                               cTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTg                                                                                                                       >  2:295409/1‑147 (MQ=255)                               cTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTg                                                                                                                       <  1:295409/147‑1 (MQ=255)                                  gACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGgcag                                                                                                                  <  7:44146/149‑1 (MQ=255)                                   aCCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGgcag                                                                                                                  >  5:66555/1‑148 (MQ=255)                                    ccGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAAATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGgcagc                                                                                                                 >  1:264098/1‑148 (MQ=255)                                    ccGTAGTGCGCGGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGgcagc                                                                                                                 >  5:283553/1‑148 (MQ=255)                                                       gCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGcgtc                                                                                             <  2:193326/149‑1 (MQ=255)                                                       gCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGcgtc                                                                                             <  3:311452/149‑1 (MQ=255)                                                              ccTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGCGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTgctg                                                                                      <  5:282729/149‑1 (MQ=255)                                                               cTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTgctgg                                                                                     >  2:85040/1‑149 (MQ=255)                                                                tccatccaGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTgctggg                                                                                    <  8:166784/149‑1 (MQ=255)                                                                    tccaGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATc                                                                                                           <  4:127254/122‑1 (MQ=255)                                                                    tccaGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATc                                                                                                           >  3:127254/1‑122 (MQ=255)                                                                      caGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGt                                                                               <  4:190223/148‑1 (MQ=255)                                                                      caGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGt                                                                               <  6:66555/148‑1 (MQ=255)                                                                      caGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGt                                                                               <  5:327036/148‑1 (MQ=255)                                                                      caGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGt                                                                               >  6:327036/1‑148 (MQ=255)                                                                       aGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCg                                                                             >  5:218234/1‑149 (MQ=255)                                                                          aaGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGtt                                                                          <  1:85040/149‑1 (MQ=255)                                                                          aaGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGtt                                                                          <  1:196692/149‑1 (MQ=255)                                                                                tgatgaTCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAg                                                                    <  6:199085/149‑1 (MQ=255)                                                                                  atgatCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAgg                                                                   <  6:283553/148‑1 (MQ=255)                                                                                      tCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGgcggcg                                                              >  3:75932/1‑149 (MQ=255)                                                                                      tCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGgcggcg                                                              >  2:298622/1‑149 (MQ=255)                                                                                          cccGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTg                                                          <  6:112165/149‑1 (MQ=255)                                                                                                  aGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTc                                                  <  1:298622/149‑1 (MQ=255)                                                                                                       ggCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACa                                              <  7:272887/148‑1 (MQ=255)                                                                                                        gCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGt                                            >  2:122030/1‑149 (MQ=255)                                                                                                         cATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGtt                                           <  2:264098/149‑1 (MQ=255)                                                                                                          aTCCTGTTGTTGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGttt                                          >  3:232189/1‑149 (MQ=255)                                                                                                          aTCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTgctgggctg                                                                                 >  6:95835/1‑110 (MQ=255)                                                                                                          aTCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTgctgggctg                                                                                 <  5:95835/110‑1 (MQ=255)                                                                                                          aTCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGt                                                   <  3:264354/140‑1 (MQ=255)                                                                                                          aTCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGttt                                          >  7:99006/1‑149 (MQ=255)                                                                                                          aTCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCACAGTTGCCAATGGGCATGATCATCGAACGCCTGGGCAGCAAGATCGCCCTCGGAGCGTCGTTGCTGCGCTGGTAGGTTGAGCAGGCGGCGTATGGAATGGt                                                   >  4:264354/1‑140 (MQ=255)                                                                                                             cTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGc                                       <  6:218234/149‑1 (MQ=255)                                                                                                                       cGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATg                             <  8:241539/149‑1 (MQ=255)                                                                                                                         ttttCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATggg                           >  3:105945/1‑149 (MQ=255)                                                                                                                           ttCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGc                         >  3:119303/1‑149 (MQ=255)                                                                                                                           ttCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTACCGGTTTCATgggg                          <  1:122030/148‑1 (MQ=255)                                                                                                                            tCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCt                        >  5:236192/1‑149 (MQ=255)                                                                                                                                ttGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCgg                    <  6:236192/149‑1 (MQ=255)                                                                                                                                      aTGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTg              >  2:153946/1‑149 (MQ=255)                                                                                                                                          tcatGCCCATTGGCAACTGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATgg                            >  7:198060/4‑131 (MQ=37)                                                                                                                                          tcatGCCCATTGGCAACTGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATgg                            <  8:198060/128‑1 (MQ=37)                                                                                                                                          cTTTGCCCGGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGCTgg          >  6:80891/1‑149 (MQ=255)                                                                                                                                          cTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGCTgg          >  6:117777/1‑149 (MQ=255)                                                                                                                                               cccAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGCTGGGGAtt     >  3:206606/1‑149 (MQ=255)                                                                                                                                                ccAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGCTGGGGATTg    >  1:295371/1‑149 (MQ=255)                                                                                                                                                ccAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAAAAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGATGGGGATTg    >  6:148730/1‑149 (MQ=255)                                                                                                                                                 cAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGCTGGGGATTgg   <  4:75932/149‑1 (MQ=255)                                                                                                                                                  aGTTGCCAATTGGCATGATCATCGACCGCCTGGGCAGCATGATCGGCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAAGAGTTTTGCCGGTTTCATgg                            <  8:298846/123‑1 (MQ=255)                                                                                                                                                  aGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGCTGGGGATTGGc  <  8:99006/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TTCCTGCTGCTGATCGGCATCGTCAACTACCTGGACCGTAGTGCGCTGTCGATCGCCAACACCTCCATCCAGAAAGACATGATGATCAGCCCGTCGCAGATGGGCATCCTGTTGTCGGCGTTTTCCATTGCCTATGCCTTTGCCCAGTTGCCAATGGGCATGATCATCGACCGCCTGGGCAGCAAGATCGCCCTCGGCGCGTCGTTGCTGGGCTGGTCGGTTGCCCAGGCGGCGTTTGGCATGGTCAACAGTTTTGCCGGTTTCATGGGGCTGCGGGTGTTGCTGGGGATTGGC  >  NC_002947/3238909‑3239202

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