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%	256.1 / 3.7	79	R1057G (CGA→GGA)	PP_4410	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base G (2/2);  total (38/41)
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.31e-01

TTGCACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATA  >  NC_002947/5006669‑5006965                                                                                                                                                      |                                                                                                                                                      ttGCACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTcc                                                                                                                                                       >  3:317604/1‑149 (MQ=255)   gCACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGa                                                                                                                                                     >  7:208220/1‑149 (MQ=255)    cACCTATCCAGTCTCCCTGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCCTGGCGCC‑GCACTTATGCCATCAAATATCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGa                                                                                                                                                     <  2:102996/148‑1 (MQ=255)     aCCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGat                                                                                                                                                    >  2:281731/1‑148 (MQ=255)            cAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTg                                                                                                                                            <  8:208220/149‑1 (MQ=255)             aGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTg                                                                                                                                            >  1:151856/1‑148 (MQ=255)               aCTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGATCGGAGCATCCCtacaa                                                                                                                                         <  7:211141/149‑5 (MQ=255)               aCTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑ACACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGATCGGAGCATCCCtacaa                                                                                                                                         <  7:58785/149‑5 (MQ=255)                 tCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑tt                                                                                                                                      >  6:26097/1‑149 (MQ=255)                  ccATGTCGAGGGAGTTTGAGTAGAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGAGCATTCAAACTTGTAGGGCTGCTCCGATCCTCCCTGAAC‑ttt                                                                                                                                     >  8:59341/1‑149 (MQ=255)                   cATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑ttt                                                                                                                                     <  2:303158/148‑1 (MQ=255)                    aTGTCGAGGGAGTTTGAGTACAAGATTGTTCGTGCAACTCGCCTCCTCTTGATTAACGACGGAAATACGCTTGGCTGC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGatcatgc                                                                                                                                               <  2:342294/138‑3 (MQ=255)                    aTGTCGAGGGAGTTTGAGTACAAGATTGATCGAGCAACTCGCCTCCTCATGATTAACGACGGAAATACTCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGAGCATTCAAACTTGTAGGGATGCTCCGatcatac                                                                                                                                               >  1:342294/1‑136 (MQ=255)                         gAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGAc                                                                                                                                  >  1:122047/1‑145 (MQ=255)                         gAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGAc                                                                                                                                  >  8:124611/1‑145 (MQ=255)                            ggAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCg                                                                                                                           <  5:26097/149‑1 (MQ=255)                             gAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGa                                                                                                                          <  3:1960/149‑1 (MQ=255)                              aGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGAc                                                                                                                         <  2:122047/149‑1 (MQ=255)                                    aGTTCAAGATTGATCGTGCAACTCGCCTCATCATTATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGGGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTaa                                                                                                                    <  7:73431/148‑1 (MQ=255)                                    aGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTTAAC‑TTTGACAAGATCGACATTAAg                                                                                                                   <  2:183458/149‑1 (MQ=255)                                      tACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGAt                                                                                                                             >  2:194713/1‑137 (MQ=255)                                      tACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGAt                                                                                                                             <  1:194713/137‑1 (MQ=255)                                      tACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGa                                                                                                                 >  3:322627/1‑149 (MQ=255)                                        cAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGAtt                                                                                                               >  4:8720/1‑149 (MQ=255)                                        cAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGAtt                                                                                                               >  1:167418/1‑149 (MQ=255)                                             ttGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGt                                                                                                           <  4:317604/148‑1 (MQ=255)                                             ttGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTc                                                                                                          <  1:281731/149‑1 (MQ=255)                                                aTCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCAc                                                                                                       <  7:1703/149‑1 (MQ=255)                                                 tCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCAcc                                                                                                      >  7:29799/1‑149 (MQ=255)                                                 tCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCAcc                                                                                                      >  4:115084/1‑149 (MQ=255)                                                       aaCTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTtctc                                                                                                 >  4:133992/1‑148 (MQ=255)                                                       aaCTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCa                                                                                                <  2:119817/149‑1 (MQ=255)                                                         cTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTtctcta                                                                                               >  2:121595/1‑146 (MQ=255)                                                              ctcctcATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcga                                                                                         <  1:134252/149‑1 (MQ=255)                                                              ctcctcATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcga                                                                                         <  2:167418/149‑1 (MQ=255)                                                              ctcctcATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcga                                                                                         <  7:26329/149‑1 (MQ=255)                                                              ctcctcATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcg                                                                                          >  8:26329/1‑148 (MQ=255)                                                                 ctcATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcgacga                                                                                      >  2:229356/1‑149 (MQ=255)                                                                      gATTAACGACGGAAATACGCTTGGCGCCGGCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATa                                                                                    >  1:21365/1‑147 (MQ=255)                                                                      gATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAg                                                                                 >  6:258750/1‑149 (MQ=255)                                                                      gATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCACTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAg                                                                                 >  5:62983/1‑149 (MQ=255)                                                                           acgacgGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCTATCTCGACGATAAAGGCTgg                                                                            <  1:121595/149‑1 (MQ=255)                                                                           acgacgGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTgg                                                                            <  3:8720/149‑1 (MQ=255)                                                                           acgacgGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTg                                                                             <  2:13737/148‑1 (MQ=255)                                                                            cgacgGAAATACGCTTGGCGAC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTggg                                                                           <  5:258750/149‑1 (MQ=255)                                                                              acgGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGag                                                                         <  2:21365/149‑1 (MQ=255)                                                                                ggAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTgg                                                                            >  8:168925/1‑144 (MQ=255)                                                                                ggAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTgg                                                                            <  7:168925/144‑1 (MQ=255)                                                                                     tACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCAt                                                                  >  2:163869/1‑149 (MQ=255)                                                                                        gCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTAc                                                               >  7:126861/1‑149 (MQ=255)                                                                                          ttGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATaaa                                                                                  >  5:21162/1‑128 (MQ=255)                                                                                          ttGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATaaa                                                                                  <  6:21162/128‑1 (MQ=255)                                                                                             gcgcC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCAc                                                           <  7:59341/148‑1 (MQ=255)                                                                                             gcgcC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCAc                                                           <  3:115084/148‑1 (MQ=255)                                                                                                cc‑gcACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTtat                                                       <  4:146670/149‑1 (MQ=255)                                                                                                         aTGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcgac                                                                                        >  7:109071/1‑108 (MQ=255)                                                                                                         aTGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcgac                                                                                        <  8:109071/108‑1 (MQ=255)                                                                                                         aTGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGa                                                <  7:124611/148‑1 (MQ=255)                                                                                                          tGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGAcc                                              <  8:126861/149‑1 (MQ=255)                                                                                                            ccATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCaa                                            <  4:322627/149‑1 (MQ=255)                                                                                                             cATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAg                                           >  7:44056/1‑149 (MQ=255)                                                                                                              aTCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGt                                          >  7:151945/1‑149 (MQ=255)                                                                                                                cAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTc                                         >  6:6641/1‑148 (MQ=255)                                                                                                                  aaTTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGAc                                      <  1:229356/149‑1 (MQ=255)                                                                                                                   aTTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGagca                                                                       <  1:163177/115‑1 (MQ=255)                                                                                                                   aTTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGagca                                                                       >  2:163177/1‑115 (MQ=255)                                                                                                                       cccTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCt                                  <  2:259131/148‑1 (MQ=255)                                                                                                                        ccTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAACGTTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCt                                  >  1:256731/1‑148 (MQ=255)                                                                                                                        ccTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGa                                <  6:62983/149‑1 (MQ=255)                                                                                                                          tAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATg                              >  8:3677/1‑149 (MQ=255)                                                                                                                            gcgcATTCAAACTTGTAGGGATGCTCCGAGCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATggg                            >  7:230247/1‑148 (MQ=255)                                                                                                                                aTTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCt                           >  7:23116/1‑146 (MQ=255)                                                                                                                                aTTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCt                           <  8:23116/146‑1 (MQ=255)                                                                                                                                         ttcaGGGATGATCGGAGCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAAtt                >  3:348405/4‑148 (MQ=255)                                                                                                                                           tAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcga                                                                                         >  2:2299/1‑73 (MQ=255)                                                                                                                                           tAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcga                                                                                         <  1:2299/73‑1 (MQ=255)                                                                                                                                             gggATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCaaaa           <  3:133992/149‑1 (MQ=255)                                                                                                                                             gggATGATCGGAGCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCaaaa           >  3:83516/1‑149 (MQ=255)                                                                                                                                               gATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATc         >  8:69496/1‑149 (MQ=255)                                                                                                                                                 tGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCa       <  8:29799/149‑1 (MQ=255)                                                                                                                                                      cGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCata  <  5:6641/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TTGCACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATA  >  NC_002947/5006669‑5006965

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