Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,332,965	G→A	5.2%	A257T (GCC→ACC)	PP_3804 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,332,965	0	G	A	5.2%	242.1 / 2.0	77	A257T (GCC→ACC)	PP_3804	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base A (2/2);  total (37/40)
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.87e-01

GGTCGAGCAGCCACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGC  >  NC_002947/4332823‑4333107                                                                                                                                               |                                                                                                                                                ggTCGAGCAGCCACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGc                                                                                                                                           >  7:61138/1‑149 (MQ=255) ggTCGAGCAGCCACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGc                                                                                                                                           <  6:142803/149‑1 (MQ=255)  gTCGAGCAGCCACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCt                                                                                                                                          >  3:123328/1‑149 (MQ=255)   tCGAGCAGCCACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCtg                                                                                                                                         <  7:158666/149‑1 (MQ=255)      agcagcCACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCtggtg                                                                                                                                      <  7:165850/149‑1 (MQ=255)          gcCACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCtggtggtgg                                                                                                                                  <  4:201913/149‑1 (MQ=255)            cACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCtggtggtggtg                                                                                                                                <  2:95925/149‑1 (MQ=255)              cTGCGGGCGGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACGTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGga                                                                                                                              <  8:46964/149‑1 (MQ=255)              cTGCCGGCCGAGTGGGATGCGGCTGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGga                                                                                                                              <  8:40836/149‑1 (MQ=255)                 ccggccgAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGgaga                                                                                                                            <  8:15475/148‑1 (MQ=255)                 ccggccgAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGgaga                                                                                                                            <  3:14840/148‑1 (MQ=255)                 ccggccgAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGgaga                                                                                                                            <  2:334014/148‑1 (MQ=255)                    gccgAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGa                                                                                                                        >  2:259245/1‑149 (MQ=255)                      cgAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGGGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGAcc                                                                                                                      <  4:123328/149‑1 (MQ=255)                         gTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTg                                                                                                                   >  3:342290/1‑149 (MQ=255)                            ggATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTg                                                                                                                         <  4:143393/140‑1 (MQ=255)                            ggATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTg                                                                                                                         >  3:143393/1‑140 (MQ=255)                              aTGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGAt                                                                                                              >  8:162657/1‑149 (MQ=255)                               tGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATg                                                                                                             >  3:302705/1‑149 (MQ=255)                                gCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCtggtg                                                                                                                                      <  8:199576/123‑1 (MQ=255)                                gCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCtggtg                                                                                                                                      >  7:199576/1‑123 (MQ=255)                                         cTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCggg                                                                                                   >  4:184205/1‑149 (MQ=255)                                         cTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCggg                                                                                                   >  8:297492/1‑149 (MQ=255)                                         cTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCggg                                                                                                   >  1:208634/1‑149 (MQ=255)                                         cTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCggg                                                                                                   >  2:264129/1‑149 (MQ=255)                                         cTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGACCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCggg                                                                                                   >  8:221377/1‑149 (MQ=255)                                         cTCAACGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGACCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCggg                                                                                                   >  8:215695/1‑149 (MQ=255)                                             aGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCt                                                                                                >  6:138724/1‑148 (MQ=255)                                               gcgcTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCAt                                                                                             <  6:238257/149‑1 (MQ=255)                                                 gcTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCt                                                                                                >  2:387248/1‑144 (MQ=255)                                                 gcTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCt                                                                                                <  1:387248/144‑1 (MQ=255)                                                    gggggAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGACCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCa                                                                                        <  7:221377/149‑1 (MQ=255)                                                     ggggAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGACCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCa                                                                                        <  7:215695/148‑1 (MQ=255)                                                       ggAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCg                                                                                     >  6:35025/1‑149 (MQ=255)                                                       ggAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCg                                                                                     <  1:12351/149‑1 (MQ=255)                                                        gagaACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGt                                                                                    <  6:43370/149‑1 (MQ=255)                                                         agaACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCGTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCGAAGCTGGGGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGtt                                                                                   <  8:97129/149‑1 (MQ=255)                                                           aaCCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGa                                                                                 <  5:35025/149‑1 (MQ=255)                                                             ccTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGAcc                                                                               <  5:228999/149‑1 (MQ=255)                                                              cTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGAcc                                                                               >  4:276142/1‑148 (MQ=255)                                                                gAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGAGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAg                                                                            >  4:61690/1‑149 (MQ=255)                                                                     cgcAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGtt                                                                       <  2:5349/149‑1 (MQ=255)                                                                      gcAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATGGCCAGCGTTGACCAAGTGGTTg                                                                      >  7:161581/1‑149 (MQ=255)                                                                       cAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATGGCCAGCGTTGACCAAGTGGTTgg                                                                     <  8:161581/149‑1 (MQ=255)                                                                        aGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTggg                                                                    <  3:276142/149‑1 (MQ=255)                                                                                        gCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCt                                                   >  7:299358/1‑149 (MQ=255)                                                                                        gCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCt                                                   >  7:115520/1‑149 (MQ=255)                                                                                          tGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGCCCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGa                                                 >  8:62517/1‑149 (MQ=255)                                                                                           ggACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTgg                                                            <  4:141048/137‑1 (MQ=255)                                                                                           ggACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTgg                                                            >  3:141048/1‑137 (MQ=255)                                                                                            gACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATc                                               >  3:280638/1‑149 (MQ=255)                                                                                            gACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATc                                               <  2:208634/149‑1 (MQ=255)                                                                                             aCCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCg                                              <  5:282828/149‑1 (MQ=255)                                                                                             aCCATCGTCAGCCTGAATCGGCAGGGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGGTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCg                                              <  5:65632/149‑1 (MQ=255)                                                                                              ccATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTgggg                                                                   <  2:248389/128‑1 (MQ=255)                                                                                              ccATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTgggg                                                                   >  1:248389/1‑128 (MQ=255)                                                                                                aTCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGAc                                           <  3:61690/149‑1 (MQ=255)                                                                                                       gCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCcg                                    >  8:71043/1‑149 (MQ=255)                                                                                                        ccTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCcgc                                   >  4:328097/1‑149 (MQ=255)                                                                                                          tGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGCGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCcgcgc                                 >  6:114600/1‑149 (MQ=255)                                                                                                            aaTCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCgg                               >  2:143522/1‑149 (MQ=255)                                                                                                              tCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTg                             <  1:259245/149‑1 (MQ=255)                                                                                                                ggCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGaa                                                                                                                 >  6:47576/1‑64 (MQ=255)                                                                                                                ggCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGaa                                                                                                                 <  5:47576/64‑1 (MQ=255)                                                                                                                 gCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGAt                          >  3:35284/1‑149 (MQ=255)                                                                                                                  cAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGAt                          <  1:264129/149‑1 (MQ=255)                                                                                                                    gTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGc                       <  5:138724/149‑1 (MQ=255)                                                                                                                         cgtggTGATCAAGGCGGCCGGGGCCAAGCTGGGGGGGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGt                              <  8:170535/134‑1 (MQ=255)                                                                                                                         cGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGCTCGGACAGGCCCGCGCGGt                              >  7:170535/1‑137 (MQ=255)                                                                                                                           aGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCt                <  5:159311/149‑1 (MQ=255)                                                                                                                                 tCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGAt                          <  6:144025/133‑1 (MQ=255)                                                                                                                                 tCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGAt                          >  5:144025/1‑133 (MQ=255)                                                                                                                                   aaGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCgt        <  3:184205/149‑1 (MQ=255)                                                                                                                                     ggcggcCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCgtgt      >  3:59106/1‑149 (MQ=255)                                                                                                                                      gcggcCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGtt     <  1:143522/149‑1 (MQ=255)                                                                                                                                        ggcCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTcg   >  1:381667/1‑149 (MQ=255)                                                                                                                                         gcCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTcgc  >  3:195235/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                GGTCGAGCAGCCACTGCCGGCCGAGTGGGATGCGGCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCA‑TTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGC  >  NC_002947/4332823‑4333107

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