Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,798,380	C→A	5.5%	D495Y (GAC→TAC) ‡	PP_5077 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,798,380	0	C	A	5.5%	234.2 / 3.2	72	D495Y (GAC→TAC) ‡	PP_5077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base A (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

TTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCACCACATAA  >  NC_002947/5798234‑5798525                                                                                                                                                   |                                                                                                                                                  ttAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgc                                                                                                                                                 >  4:47125/1‑149 (MQ=255) ttAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgc                                                                                                                                                 >  3:302880/1‑149 (MQ=255) ttAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCATGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgc                                                                                                                                                 >  5:170437/1‑149 (MQ=255)        cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                                                                                                                          <  4:172200/149‑1 (MQ=255)         tGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                                                                                                                          >  7:172446/1‑148 (MQ=255)           ggcggcGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCGACGCAGCGGTTgcgg                                                                                                                                         <  8:217143/147‑4 (MQ=255)            gcggcgGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGaaa                                                                                                                                      >  2:196802/1‑149 (MQ=255)               gcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTg                                                                                                                                   >  4:176108/1‑149 (MQ=255)                         gccgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCa                                                                                                                         >  8:118853/1‑149 (MQ=255)                          ccgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCac                                                                                                                        >  1:140856/1‑149 (MQ=255)                           cgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCac                                                                                                                        >  1:26630/1‑148 (MQ=255)                                ggAAGGTATCAGCGGGCCGTGGCCAGCTCTTGGTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTa                                                                                                                  <  6:323988/149‑1 (MQ=255)                                       aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTgg                                                                                                                                           <  6:127234/117‑1 (MQ=255)                                       aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTgg                                                                                                                                           >  5:127234/1‑117 (MQ=255)                                        tCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCt                                                                                                          >  7:100244/1‑149 (MQ=255)                                          aGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTg                                                                                                        <  2:26630/149‑1 (MQ=255)                                           gCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGc                                                                                                       <  6:107923/149‑1 (MQ=255)                                                ccGTGGCCAGCTCTTGTTGGACGCTGCCTACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGc                                                                                                  <  3:176108/149‑1 (MQ=255)                                                ccGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTg                                                                                                   >  2:264269/1‑148 (MQ=255)                                                      ccAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGttt                                                                                            >  7:232726/1‑149 (MQ=255)                                                             ttgttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCtt                                                                                       <  8:100244/147‑1 (MQ=255)                                                             ttgttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                                      >  5:74494/1‑148 (MQ=255)                                                             ttgttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                                      >  6:29713/1‑148 (MQ=255)                                                               gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                                    <  4:302880/148‑1 (MQ=255)                                                               gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                                    <  7:169144/148‑1 (MQ=255)                                                                ttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAg                                                                                   >  4:173795/1‑148 (MQ=255)                                                                   gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAg                                                                                <  1:196802/148‑1 (MQ=255)                                                                   gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAg                                                                                <  7:118853/148‑1 (MQ=255)                                                                   gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAg                                                                                <  1:264269/148‑1 (MQ=255)                                                                      gCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGAt                                                                            >  6:62020/1‑149 (MQ=255)                                                                          ccGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTc                                                                        >  4:250377/1‑149 (MQ=255)                                                                           cGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTct                                                                       <  4:174688/149‑2 (MQ=255)                                                                           cGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCg                                                                       <  4:174689/149‑1 (MQ=255)                                                                                    gTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTgg                                                              >  8:278964/1‑149 (MQ=255)                                                                                       gCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                                      >  5:61785/1‑122 (MQ=255)                                                                                       gCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                                      <  6:61785/122‑1 (MQ=255)                                                                                         cAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                         >  2:104641/1‑149 (MQ=255)                                                                                         cAAGGTTTACCAGCCTGGACATTAGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                         >  2:221973/1‑149 (MQ=255)                                                                                             gTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGcc                                                                     <  8:61970/133‑1 (MQ=255)                                                                                             gTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGcc                                                                     >  7:61970/1‑133 (MQ=255)                                                                                                 aCCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTaa                                                  <  4:118015/148‑1 (MQ=255)                                                                                                 aCCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTaa                                                  <  3:47125/148‑1 (MQ=255)                                                                                                   cAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAg                                                 >  5:40998/1‑147 (MQ=255)                                                                                                   cAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAg                                                 <  6:40998/147‑1 (MQ=255)                                                                                                     gCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTg                                             <  6:74494/149‑1 (MQ=255)                                                                                                     gCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTg                                             >  8:291076/1‑149 (MQ=255)                                                                                                      ccTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTgg                                            <  7:31967/149‑1 (MQ=255)                                                                                                       cTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTggg                                                             >  8:246000/1‑131 (MQ=255)                                                                                                       cTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTggg                                                             <  7:246000/131‑1 (MQ=255)                                                                                                        tGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACTCTGGGCCTTGACGTAAGCCTGGGc                                          >  8:314457/1‑149 (MQ=255)                                                                                                                cGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCt                                  >  4:233757/1‑149 (MQ=255)                                                                                                                 gCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTc                                 <  5:62020/149‑1 (MQ=255)                                                                                                                      cAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGc                            >  6:80401/1‑149 (MQ=255)                                                                                                                       aaGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGCGCt                           >  8:301537/1‑149 (MQ=255)                                                                                                                       aaGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCt                           <  6:170437/149‑1 (MQ=255)                                                                                                                       aaGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCt                           <  8:172446/149‑1 (MQ=255)                                                                                                                        aGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTg                          <  3:233757/149‑1 (MQ=255)                                                                                                                          ttCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTgg                         >  8:154249/1‑148 (MQ=255)                                                                                                                             ttGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCg                                <  5:52493/138‑1 (MQ=255)                                                                                                                             ttGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCg                                >  6:52493/1‑138 (MQ=255)                                                                                                                               gATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTTGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTAcc                     <  4:232747/147‑1 (MQ=255)                                                                                                                                aTTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCaa                   <  5:80401/148‑1 (MQ=255)                                                                                                                                 ttGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCa                    >  2:82276/1‑146 (MQ=255)                                                                                                                                 ttGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCa                    <  1:82276/146‑1 (MQ=255)                                                                                                                                 ttGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAg                  >  4:205453/1‑148 (MQ=255)                                                                                                                                 ttGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAg                  >  3:93282/1‑148 (MQ=255)                                                                                                                                   gCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGa                <  1:104641/148‑1 (MQ=255)                                                                                                                                    ccGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGAt               <  3:250377/148‑1 (MQ=255)                                                                                                                                     gccAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCt             >  7:290051/3‑149 (MQ=255)                                                                                                                                      ccAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                         <  8:75426/103‑1 (MQ=255)                                                                                                                                      ccAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                         >  7:75426/2‑104 (MQ=255)                                                                                                                                        aaCCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGca          >  1:135577/1‑149 (MQ=255)                                                                                                                                         aCCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGcac         <  7:278964/149‑1 (MQ=255)                                                                                                                                                cGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCACCACAt    <  8:232726/147‑1 (MQ=255)                                                                                                                                                 gTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGTCTCGGAGCTGGTACCAAGGATCTGCACCACATaa  <  1:221973/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCACCACATAA  >  NC_002947/5798234‑5798525

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