Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,013,766	A→C	6.0%	V206G (GTG→GGG) ‡	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,766	0	A	C	6.0%	254.6 / 6.5	84	V206G (GTG→GGG) ‡	potF‑I	putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (38/41);  new base C (2/3);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

CCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGCTCATGTACTC  >  NC_002947/1013623‑1013914                                                                                                                                                            |                                                                                                                                                      ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtc                                                                                                                                                  <  1:9143/149‑1 (MQ=255) ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtc                                                                                                                                                  <  6:259063/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                                                                                                                                   <  2:387528/136‑1 (MQ=255)                             ggCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                                   >  1:387528/1‑136 (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                                          <  5:355361/149‑1 (MQ=255)                                                                                                                   cGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCa                                          <  2:305118/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                            <  1:246668/149‑1 (MQ=255)                                                                                                                                  aGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca                            >  3:8365/1‑149 (MQ=255)                                                                                                                                   gTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca                            >  4:225200/1‑148 (MQ=255)                                                                                                                                      gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCa                                                                                     <  5:200888/89‑1 (MQ=255)                                                                                                                                      gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCa                                                                                     >  6:200888/1‑89 (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:306675/149‑1 (MQ=255)                                                                                                                                            gCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCa                  <  3:157662/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)                                                                                                                                                          gCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATc                                                                   >  7:336584/1‑86 (MQ=255)                                                                                                                                                          gCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATc                                                                   <  8:336584/86‑1 (MQ=255)                                                                                                                                                           cACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGCTCATGTACt   >  5:107970/1‑149 (MQ=255)                                                                                                                                                            aCTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGCTCATGTACTc  <  8:8920/149‑1 (MQ=255)                                                                                                                                                            |                                                                                                                                                      CCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCC‑‑‑CT‑GAGAACCC‑‑‑‑‑‑CACTACCACGCAGATATCGCCGTTGG‑CCAGGTCGCCCGTGTACT‑TGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTA‑GTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCTTGCTCATGTACTC  >  NC_002947/1013623‑1013914

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