Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	496,096	T→A	9.2%	intergenic (‑63/‑33)	PP_0408 ← / → PP_0409	hypothetical protein/two‑component system sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	496,096	0	T	A	9.2%	172.6 / 6.4	62	intergenic (‑63/‑33)	PP_0408/PP_0409	hypothetical protein/two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base A (4/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCT  >  NC_002947/495949‑496244                                                                                                                                                    |                                                                                                                                                      tttCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgct                                                                                                                                                       <  3:284237/148‑1 (MQ=255)   tCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCTATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgctgc                                                                                                                                                     <  5:133578/148‑1 (MQ=255)      gCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATgg                                                                                                                                                 <  2:157217/149‑1 (MQ=255)        aTCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc                                                                                                                                               >  7:245551/1‑149 (MQ=255)         tCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc                                                                                                                                               >  6:276128/1‑148 (MQ=255)         tCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc                                                                                                                                               <  5:276128/148‑1 (MQ=255)              ggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTGGCCCGTGACGGCCATGCTGCATGGCCGTCttt                                                                                                                                         <  1:45022/149‑1 (MQ=255)               gcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTcacgg                                                                                                                                        <  3:272121/149‑5 (MQ=255)               gcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcacgg                                                                                                                                        <  4:115407/149‑5 (MQ=255)                    ggCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAg                                                                                                                                   <  8:245551/149‑1 (MQ=255)                         cGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGAGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAg                                                                                                                              <  5:295696/149‑1 (MQ=255)                             tGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGta                                                                                                                          >  6:243582/1‑149 (MQ=255)                               gAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGta                                                                                                                          >  6:181486/1‑147 (MQ=255)                                aaGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGta                                                                                                                          >  1:76043/1‑146 (MQ=255)                                  ggACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGtata                                                                                                                        <  6:233595/146‑1 (MQ=255)                                  ggACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGtata                                                                                                                        >  5:233595/1‑146 (MQ=255)                                    aCCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGGATTTGCCCGTGACGGCCATGCTGCATGGCCGTTTTTGCCAAGTCGAGCGTATAGATg                                                                                                                    >  5:96450/1‑148 (MQ=255)                                     ccGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATg                                                                                                                    >  2:274601/1‑147 (MQ=255)                                      cGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAg                                                                                                                 >  2:105804/1‑149 (MQ=255)                                         ggCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGa                                                                                                              <  4:51859/149‑1 (MQ=255)                                         ggCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGa                                                                                                              <  1:274601/149‑1 (MQ=255)                                             aGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTAGTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGc                                                                                                          >  7:242657/1‑149 (MQ=255)                                                      aTTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGACTTTGCCAAGTCGAGCGTATAGATGAAGAGAGTGGGTCGCCTg                                                                                                  >  7:307388/1‑148 (MQ=255)                                                       ttGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGtt                                                                                                <  6:131675/149‑1 (MQ=255)                                                       ttGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGtt                                                                                                <  3:70142/149‑1 (MQ=255)                                                             gCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑t                                                                                         >  2:29668/1‑149 (MQ=255)                                                                  gtGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGgct                                                                                    >  7:255283/1‑149 (MQ=255)                                                                        tAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGgcagccac                                                                               >  4:48764/1‑145 (MQ=255)                                                                             tgtgGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTg                                                                         >  3:193093/1‑149 (MQ=255)                                                                               tgGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGcc                                                                       >  5:49220/1‑149 (MQ=255)                                                                                  aCATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCtt                                                                    <  5:181486/149‑1 (MQ=255)                                                                                        gTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGt                                                                                                                                  <  6:117199/82‑1 (MQ=255)                                                                                        gTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGt                                                                                                                                  >  5:117199/1‑82 (MQ=255)                                                                                         tCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGgc                                                                                     >  8:4563/1‑125 (MQ=255)                                                                                         tCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGgc                                                                                     <  7:4563/125‑1 (MQ=255)                                                                                            aCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTTTTTGCCAAGGCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTgg                                                          <  6:96450/149‑1 (MQ=255)                                                                                            aCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCaa                                                                                                                                    >  3:127681/1‑76 (MQ=255)                                                                                            aCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCaa                                                                                                                                    <  4:127681/76‑1 (MQ=255)                                                                                              ccATGGGCATAATGAATCCGGCTGCCTAATTGTATTTGCCCGTGTCGGCGATGCTGCTTGGCCGTCTTTGCCTAGTCGAGCGTATAGATGAAGCGAGTGGGTCGCCTGTTGGTTT‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGc                                                        <  8:307388/149‑1 (MQ=255)                                                                                                aaGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGACATGCAGCATGGCCCTCTTTCCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTAA‑TCGGCAGCTTGTGGCTGCCCTTGATTGCCTGGGcc                                                       >  1:83061/1‑148 (MQ=255)                                                                                                 aGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCaa                                                     <  8:242657/149‑1 (MQ=255)                                                                                                  ggCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAg                                                    >  6:309677/1‑149 (MQ=255)                                                                                                   gCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGc                                                   >  8:97671/1‑149 (MQ=255)                                                                                                       gAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGc                                                   <  4:20689/145‑1 (MQ=255)                                                                                                       gAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGc                                                   >  3:20689/1‑145 (MQ=255)                                                                                                         aTGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAgg                                             <  5:243582/149‑1 (MQ=255)                                                                                                             aaCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCt                                         <  7:97671/149‑1 (MQ=255)                                                                                                                 ggCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGgc                                            <  5:367952/142‑1 (MQ=255)                                                                                                                 ggCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGgc                                            >  6:367952/1‑142 (MQ=255)                                                                                                                    aGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCctct                                  <  2:76043/149‑1 (MQ=255)                                                                                                                      ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCctct                                  <  5:319649/147‑1 (MQ=255)                                                                                                                      ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCctct                                  >  6:319649/1‑147 (MQ=255)                                                                                                                      ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgt                                >  8:284121/1‑149 (MQ=255)                                                                                                                       cTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgtg                               <  4:193093/149‑1 (MQ=255)                                                                                                                          aTAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATTCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTgtggt                             <  1:105804/149‑1 (MQ=255)                                                                                                                                tttGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTg                         <  3:353954/146‑1 (MQ=255)                                                                                                                                tttGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTg                         >  4:353954/1‑146 (MQ=255)                                                                                                                                tttGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGa                       <  1:29668/148‑1 (MQ=255)                                                                                                                                      cGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGGCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc                <  8:255283/149‑1 (MQ=255)                                                                                                                                      aaagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTgtggtg                            >  4:142895/4‑137 (MQ=255)                                                                                                                                      aaagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTgtggtg                            <  3:142895/134‑1 (MQ=255)                                                                                                                                      aaagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc                >  6:253638/4‑149 (MQ=255)                                                                                                                                          aCGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcag             >  2:202035/1‑148 (MQ=255)                                                                                                                                            ggCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcag          >  5:185320/1‑149 (MQ=255)                                                                                                                                                  gctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTgg    >  8:110738/1‑149 (MQ=255)                                                                                                                                                   ctgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTgg    <  5:309677/148‑1 (MQ=255)                                                                                                                                                    tgcATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCt  >  8:107029/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TTTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTA‑TCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCT  >  NC_002947/495949‑496244

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑