Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,013,757	T→A	5.9%	Y209F (TAC→TTC)	potF‑I ←	putrescine‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,013,757	0	T	A	5.9%	256.1 / 4.7	84	Y209F (TAC→TTC)	potF‑I	putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (39/40);  new base A (2/3);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

TATCGAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGC  >  NC_002947/1013612‑1013904                                                                                                                                                              |                                                                                                                                                     tATCGAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAg                                                                                                                                                   >  1:280211/1‑149 (MQ=255)   tCGAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑ACAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGc                                                                                                                                                 >  7:27565/1‑149 (MQ=255)     gAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTATGGCCg                                                                                                                                               >  4:331629/1‑149 (MQ=255)            ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtc                                                                                                                                        <  6:259063/149‑1 (MQ=255)            ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtc                                                                                                                                        <  1:9143/149‑1 (MQ=255)             cACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtca                                                                                                                                       >  7:318673/1‑149 (MQ=255)              aCCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCCCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtcat                                                                                                                                      >  6:168941/1‑149 (MQ=255)                cTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGAAAGTGCGGCCTTacatca                                                                                                                                    <  1:96982/149‑1 (MQ=255)                   cACGCCTTTGTTCGCCTCGTCGGGGGGGTTCTTCGCCTGCAGCACATCCCTTTCTCTTCATCCCTTTTAACTCTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCa                                                                                                                                             <  1:334597/147‑1 (MQ=16)                        cATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACTCTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTc                                                                                                                             <  1:299876/149‑1 (MQ=255)                        cATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCACGTcc                                                                                                                            >  2:394520/1‑149 (MQ=255)                         aTTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGAAAGTGCGGCCTTACATcagctatttcc                                                                                                                           <  4:354127/149‑8 (MQ=38)                            gTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         <  1:326040/148‑1 (MQ=255)                            gTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATCTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         <  1:34315/148‑1 (MQ=255)                                gCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCtt                                                                                                                    >  3:10924/1‑149 (MQ=255)                                 ccTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTtc                                                                                                                   >  1:195625/1‑149 (MQ=255)                                 ccTCGTCGGCGCGGTTCTTCGCCTGCAGCACAAACC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGAGGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCCGGTCCTTGGCCTtc                                                                                                                   >  6:136518/1‑149 (MQ=255)                                     gtcgGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTtctc                                                                                                                 <  1:81382/147‑1 (MQ=255)                                     gtcgGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTtctc                                                                                                                 <  4:400854/147‑1 (MQ=255)                                     gtcgGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTAACACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCg                                                                                                                <  1:310840/148‑1 (MQ=255)                                      tcgGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTa                                                                                                              <  3:1871/149‑1 (MQ=255)                                        ggCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         >  1:387528/1‑136 (MQ=255)                                        ggCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         <  2:387528/136‑1 (MQ=255)                                              gTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATGGCTGATGTAAGGCCGCACTTTCATTAGCAGGTCCTTGGCCTTCTCGTAGTCCTTgg                                                                                                      >  7:222488/1‑149 (MQ=255)                                              gTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         <  7:133244/130‑1 (MQ=255)                                              gTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                         >  8:133244/1‑130 (MQ=255)                                              gTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCATGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTgg                                                                                                      >  1:191904/1‑149 (MQ=255)                                               ttcttcGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTggg                                                                                                     <  8:318673/149‑1 (MQ=255)                                                      ccTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGGACCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTg                                                                                               >  1:305118/1‑149 (MQ=255)                                                            gcaCATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGa                                                                                        <  5:244972/149‑1 (MQ=255)                                                              acaTCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATg                                                                                      >  4:306675/1‑149 (MQ=255)                                                               caTCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTGCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATgg                                                                                     >  8:155826/1‑149 (MQ=255)                                                                aTCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATggg                                                                                    <  2:195625/149‑1 (MQ=255)                                                                aTCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACAACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATggg                                                                                    >  8:138047/1‑149 (MQ=255)                                                                   cc‑‑‑cT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGc                                                                                  >  4:157662/1‑148 (MQ=255)                                                                           gagaACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGcc                                                                             >  6:355361/1‑149 (MQ=255)                                                                               aCCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAgg                                                                         <  4:68296/149‑1 (MQ=255)                                                                               aCCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCATGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAgg                                                                         <  2:191904/149‑1 (MQ=255)                                                                                ccc‑‑‑‑‑‑cACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACTAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGTTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGt                                                                        <  1:210433/149‑1 (MQ=255)                                                                                 cc‑‑‑‑‑‑cACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTa                                                                       >  2:246668/1‑149 (MQ=255)                                                                                  c‑‑‑‑‑‑cACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑g                                                                     >  3:317647/1‑149 (MQ=255)                                                                                           cTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGc                                                                  <  6:390894/149‑1 (MQ=255)                                                                                              ccACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATGGCTGATGTAAGGCCGCACTTTCATTAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAg                                                                <  8:222488/148‑1 (MQ=255)                                                                                               cACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAggg                                                              <  4:10924/149‑1 (MQ=255)                                                                                               cACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAggg                                                              <  5:168941/149‑1 (MQ=255)                                                                                                aCGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAg                                                                          >  2:377681/1‑137 (MQ=255)                                                                                                aCGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAg                                                                          <  1:377681/137‑1 (MQ=255)                                                                                                aCGCAGATATCGCCGTTCG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTAGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑CTGCAGGGc                                                             >  7:106501/1‑149 (MQ=255)                                                                                                    aGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATTGGGCAAGCCCAGGTA‑GTGCAgagatcc                                                          >  1:143120/1‑142 (MQ=255)                                                                                                    aGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCa                                                                 <  2:143120/141‑1 (MQ=255)                                                                                                    aGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATc                                                         <  7:246205/149‑1 (MQ=255)                                                                                                       tatCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCgg                                                       >  6:379801/1‑148 (MQ=255)                                                                                                        atCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCgg                                                       >  6:50163/1‑147 (MQ=255)                                                                                                         tCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGgcag                                                    <  4:317647/149‑1 (MQ=255)                                                                                                           gCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGgcag                                                    >  6:29087/1‑147 (MQ=255)                                                                                                             cGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTg                                                                                               <  6:107621/103‑1 (MQ=255)                                                                                                             cGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTg                                                                                               >  5:107621/1‑103 (MQ=255)                                                                                                              gTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGtt                                               <  5:29087/149‑1 (MQ=255)                                                                                                               ttGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTc                                              <  2:280211/149‑1 (MQ=255)                                                                                                               ttGG‑ACAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTc                                              <  8:27565/149‑1 (MQ=255)                                                                                                                    ccAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCgggg                                          >  2:2724/1‑149 (MQ=255)                                                                                                                           gCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGTTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGt                                   >  7:8920/1‑149 (MQ=255)                                                                                                                              cGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCa                                <  2:305118/149‑1 (MQ=255)                                                                                                                              cGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCa                                <  5:355361/149‑1 (MQ=255)                                                                                                                               gtgtACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGGGATGTGGCCAGGCCAGGTA‑GTGCAGGGGGATCGGCAGCAGTTCGGGGCCATTGtgaag                               >  6:337366/1‑145 (MQ=255)                                                                                                                                     t‑tGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTACGTGCAGGGCGATCGGCAGCAGTTCGGGGCCAtt                                     <  1:259556/138‑1 (MQ=255)                                                                                                                                     t‑tGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCAtt                                     >  2:259556/1‑137 (MQ=255)                                                                                                                                             aGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca                  >  3:8365/1‑149 (MQ=255)                                                                                                                                             aGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca                  <  1:246668/149‑1 (MQ=255)                                                                                                                                              gTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca                  >  4:225200/1‑148 (MQ=255)                                                                                                                                                 gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCa                                                                           >  6:200888/1‑89 (MQ=255)                                                                                                                                                 gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCa                                                                           <  5:200888/89‑1 (MQ=255)                                                                                                                                                 gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGGAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGc              >  2:276205/1‑149 (MQ=255)                                                                                                                                                 gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGc              >  2:167334/1‑149 (MQ=255)                                                                                                                                                    aTAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGc                                         >  4:65738/1‑119 (MQ=255)                                                                                                                                                    aTAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGc                                         <  3:65738/119‑1 (MQ=255)                                                                                                                                                    aTAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCttt           <  5:379801/149‑1 (MQ=255)                                                                                                                                                     ctgCTGATGAAAGTGCGGCCTTACATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCttt           >  8:104510/3‑148 (MQ=255)                                                                                                                                                     ctgCTGATGAAAGTGCGGCCTTACATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCtt            >  7:372513/3‑147 (MQ=255)                                                                                                                                                     ctgCTGATGAAAGTGCGGCCTTACATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCtt            <  8:372513/145‑1 (MQ=255)                                                                                                                                                       gCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCa        <  3:157662/149‑1 (MQ=255)                                                                                                                                                       gCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCa        <  3:306675/149‑1 (MQ=255)                                                                                                                                                       gCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCa        <  5:136518/149‑1 (MQ=255)                                                                                                                                                        cTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAg       >  6:261174/1‑149 (MQ=255)                                                                                                                                                              tAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGc  <  5:41620/148‑1 (MQ=255)                                                                                                                                                              |                                                                                                                                                     TATCGAATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGC  >  NC_002947/1013612‑1013904

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