Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,006,817	C→G	5.1%	R1057G (CGA→GGA)	PP_4410 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,006,817	0	C	G	5.1%	272.2 / 3.0	87	R1057G (CGA→GGA)	PP_4410	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (43/39);  new base G (2/2);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TTGCACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATC  >  NC_002947/5006669‑5006958                                                                                                                                                     |                                                                                                                                                ttGCACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTcc                                                                                                                                                 >  4:182968/1‑149 (MQ=255)   gCACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGa                                                                                                                                               >  8:315489/1‑149 (MQ=255)    cACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTccca                                                                                                                                               >  7:91084/1‑146 (MQ=255)     aCCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGatc                                                                                                                                             <  3:265818/149‑1 (MQ=255)       cTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGatcat                                                                                                                                           <  8:40511/149‑1 (MQ=255)                 tCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATcc                                                                                                                                         >  1:75200/1‑141 (MQ=255)                 tCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATcc                                                                                                                                         <  2:75200/141‑1 (MQ=255)                 tCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACtt                                                                                                                                >  1:99169/1‑149 (MQ=255)                 tCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACtt                                                                                                                                >  7:204391/1‑149 (MQ=255)                    aTGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGa                                                                                                                             <  1:246382/149‑1 (MQ=255)                    aTGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGa                                                                                                                             >  6:66565/1‑149 (MQ=255)                     tGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGa                                                                                                                             >  3:60646/1‑148 (MQ=255)                      gTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGAc                                                                                                                            >  6:115067/1‑148 (MQ=255)                         gAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGAc                                                                                                                            >  4:219473/1‑145 (MQ=255)                         gAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAg                                                                                                                         <  3:182968/148‑1 (MQ=255)                          aGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGAc                                                                                                                            >  7:42089/1‑144 (MQ=255)                                tttGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACAt                                                                                                                 >  3:209089/1‑149 (MQ=255)                                     gTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAgg                                                                                                            <  5:381815/149‑1 (MQ=255)                                       aCAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCTCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGa                                                                                                           <  5:371435/148‑1 (MQ=255)                                       aCAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGAt                                                                                                          >  3:340499/1‑149 (MQ=255)                                        cAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATTAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATcc                                                                                                                                         >  7:200386/1‑118 (MQ=255)                                        cAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATTAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATcc                                                                                                                                         <  8:200386/118‑1 (MQ=255)                                        cAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGAtt                                                                                                         >  2:315102/1‑149 (MQ=255)                                         aaGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATTCGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGAttt                                                                                                        >  7:298253/1‑149 (MQ=255)                                         aaGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGAttt                                                                                                        >  8:209980/1‑149 (MQ=255)                                         aaGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGAttt                                                                                                        >  3:38397/1‑149 (MQ=255)                                         aaGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGAttt                                                                                                        >  2:78627/1‑149 (MQ=255)                                            aTTGATCGTGCAACTCGCATCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAAATTGTAGGGATGCTCCGATCATCCCA‑GAACTTTGACAAGATCGACATTAAGGATTTAGt                                                                                                     >  3:59573/1‑149 (MQ=255)                                             ttGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTc                                                                                                    <  5:115067/149‑1 (MQ=255)                                              tGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATTCGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTc                                                                                                    <  8:298253/148‑1 (MQ=255)                                                   gTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCtt                                                                                             >  3:256432/1‑149 (MQ=255)                                                   gTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCt                                                                                              >  1:289460/1‑148 (MQ=255)                                                     gCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTACGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGAt                                                                                                                       <  5:417613/122‑1 (MQ=255)                                                     gCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTACGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGAt                                                                                                                       >  6:417613/1‑122 (MQ=255)                                                       aaCTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTtctc                                                                                          >  4:330481/1‑148 (MQ=255)                                                       aaCTCGCCTCCTCATCATTCACCACGGAAATACGCTAGCCGCCGCACTTATGCCATCAAATTACCCTAGCGCCTTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCa                                                                                         <  6:129731/149‑1 (MQ=255)                                                          tCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAAtc                                                                                      >  3:41199/1‑149 (MQ=255)                                                            gCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAAtctc                                                                                    <  4:38397/149‑1 (MQ=255)                                                            gCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAAtct                                                                                     <  4:340499/148‑1 (MQ=255)                                                              ctcctcATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGTCAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTcga                                                                                  <  1:78627/149‑1 (MQ=255)                                                               tcctcATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTcgac                                                                                 >  2:199832/1‑149 (MQ=255)                                                                 ctcATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTcgac                                                                                 <  7:315489/147‑1 (MQ=255)                                                                    aTGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATa                                                                             >  2:194793/1‑148 (MQ=255)                                                                    aTGATTAACGACGGAAATACGCTTGGCGCCGCACATATGCCATCAAATTTCCCAAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGAAAAGATCGACAGTAAGGATTTTCTCCA‑CCTTCTCAATCACGACGATa                                                                             >  2:348735/1‑148 (MQ=255)                                                                           acgacgGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTg                                                                      <  3:330481/148‑1 (MQ=255)                                                                           acgacgGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTAGACGA‑CATACTCTTTCTGCACGATAAAgccaa                                                                      >  2:160509/1‑144 (MQ=255)                                                                           acgaagGAAATACTCTTGGCGCCGCACTTATGCCATAAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGATCCGATCATCCCT‑GAAATTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGc                                                                        <  1:160509/146‑1 (MQ=255)                                                                            cgacgGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTggg                                                                    <  4:194150/149‑1 (MQ=255)                                                                             gacgGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTTGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTgg                                                                     >  1:287235/1‑147 (MQ=255)                                                                             gacgGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAAATTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTgg                                                                     >  5:334431/1‑147 (MQ=255)                                                                                 gAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGagcaa                                                               <  3:304202/149‑1 (MQ=255)                                                                                   aaTACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTtc                                                                                            <  4:19631/118‑1 (MQ=255)                                                                                   aaTACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTtc                                                                                            >  3:19631/1‑118 (MQ=255)                                                                                   aaTACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGagcaagc                                                             <  5:66565/149‑1 (MQ=255)                                                                                             gcgcCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCAc                                                    <  7:209980/148‑1 (MQ=255)                                                                                             gcgcCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCAc                                                    <  2:99169/148‑1 (MQ=255)                                                                                               gccgcACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCAt                                                           >  8:159528/1‑139 (MQ=255)                                                                                               gccgcACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCAt                                                           <  7:159528/139‑1 (MQ=255)                                                                                               gccgcACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTta                                                 <  3:219473/149‑1 (MQ=255)                                                                                                ccgcACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTtat                                                <  4:209089/149‑1 (MQ=255)                                                                                                  gcACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGatcat                                                                                                                                           <  2:246394/58‑1 (MQ=255)                                                                                                  gcACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGatcat                                                                                                                                           >  1:246394/1‑58 (MQ=255)                                                                                                         tGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGAc                                        <  1:199832/148‑1 (MQ=255)                                                                                                         tGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCACT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGTTATAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGAcc                                       <  1:103692/149‑1 (MQ=255)                                                                                                           ccATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCATACTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGAcc                                       <  1:315102/148‑1 (MQ=255)                                                                                                            cATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCGGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAg                                    <  5:348955/149‑1 (MQ=255)                                                                                                            cATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAg                                    >  7:322740/1‑149 (MQ=255)                                                                                                             aTCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCAc                                                    <  3:344884/132‑1 (MQ=255)                                                                                                             aTCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCAc                                                    >  4:344884/1‑132 (MQ=255)                                                                                                              tCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTc                                  >  8:243251/1‑149 (MQ=255)                                                                                                                 aaTTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGAc                               >  6:362765/1‑149 (MQ=255)                                                                                                                      gagtAGGGCTTTGAAAGTTTTTGGGATGCTCTGATTATCCCT‑TTAGTTTGACAAGATTGACATTATGGATTTTGTCCA‑CCTTCTCTATCTCGACGTTATAGGCTGGGAGCAAGCATTAGTCAGTTATTAGCCGACCAAGTCGACGTCt                           <  8:322740/145‑1 (MQ=255)                                                                                                                          aGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATgg                      >  1:166509/1‑149 (MQ=255)                                                                                                                           gcgcATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATg                       <  1:194793/147‑1 (MQ=255)                                                                                                                             gcATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTc                   >  4:26533/1‑149 (MQ=255)                                                                                                                              cATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCg                  >  3:919/1‑149 (MQ=255)                                                                                                                               aTTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGc                 >  3:398158/1‑149 (MQ=255)                                                                                                                                  cAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGccc               <  2:287235/148‑1 (MQ=255)                                                                                                                                   aaaCTTGTAGGGATGCCCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCt              >  6:323989/1‑148 (MQ=255)                                                                                                                                    aaCTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTaa            <  5:199106/149‑1 (MQ=255)                                                                                                                                     aagttcaGGGATGATCGGAGCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGACGACCAAGTCGACGTCTGATGGATCGCCCTaaa           >  8:351483/7‑149 (MQ=255)                                                                                                                                     aCTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTaaa           <  4:60646/149‑1 (MQ=255)                                                                                                                                       ttGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAAtt         <  4:41199/149‑1 (MQ=255)                                                                                                                                       gttcaGGGATGATCGGAGCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACAATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAAt          >  8:251131/5‑148 (MQ=255)                                                                                                                                           aGGGATGATCGGAGCATCCCTACAAGTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCaa      <  8:273597/149‑1 (MQ=255)                                                                                                                                            gggATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCaaaa    >  2:133964/1‑149 (MQ=255)                                                                                                                                            gggATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCaaaa    <  2:289460/149‑1 (MQ=255)                                                                                                                                              gATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATc  <  4:398158/149‑1 (MQ=255)                                                                                                                                              gATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTAACATATAAGCCGACCAAGTCGACTTCTGATGGCTCGCCCTAAATTCAAAATc  >  3:18739/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                TTGCACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCCGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCT‑GAACTTTGACAAGATCGACATTAAGGATTTTGTCCA‑CCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATC  >  NC_002947/5006669‑5006958

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