Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,006,820	T→G	6.4%	S1058A (TCA→GCA)	PP_4410 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,006,820	0	T	G	6.4%	241.5 / 5.1	76	S1058A (TCA→GCA)	PP_4410	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/39);  new base G (3/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGG  >  NC_002947/5006673‑5006968                                                                                                                                                     |                                                                                                                                                      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                                                                                                                                <  1:194713/137‑1 (MQ=255)                                  tACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGAt                                                                                                                                >  2:194713/1‑137 (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                                                                                            <  7:26329/149‑1 (MQ=255)                                                          ctcctcATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcga                                                                                            <  1:134252/149‑1 (MQ=255)                                                          ctcctcATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTcga                                                                                            <  2:167418/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                                                              <  3:115084/148‑1 (MQ=255)                                                                                         gcgcC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCAc                                                              <  7:59341/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                              <  8:23116/146‑1 (MQ=255)                                                                                                                            aTTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCt                              >  7:23116/1‑146 (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)                                                                                                                                                   gATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCatat    <  7:3677/149‑1 (MQ=255)                                                                                                                                                     tcatcCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATgg  <  8:151945/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ACCTATCCAGACTCCATGTCGAGGGAGTTTGAGTACAAGATTGATCGTGCAACTCGCCTCCTCATGATTAACGACGGAAATACGCTTGGCGCC‑GCACTTATGCCATCAAATTTCCCTAGCGCATTCAAACTTGTAGGGATGCTCCGATCATCCCTGAAC‑TTTGACAAGATCGACATTAAGGATTTTGTCCACCTTCTCAATCTCGACGATAAAGGCTGGGAGCAAGCATTACTCACTTATAAGCCGACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGG  >  NC_002947/5006673‑5006968

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