Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,061,435	T→G	7.7%	L107R (CTG→CGG) ‡	PP_3579 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,061,435	0	T	G	7.7%	187.4 / 6.8	65	L107R (CTG→CGG) ‡	PP_3579	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base G (2/3);  total (32/33)
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.77e-01

CCAGCCCGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGCGAGCGC  >  NC_002947/4061288‑4061577                                                                                                                                                    |                                                                                                                                               ccagccCGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTgg                                                                                                                                              >  1:833233/1‑150 (MQ=255)     cccGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGt                                                                                                                                          <  1:651789/150‑1 (MQ=255)      ccGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTg                                                                                                                                         >  1:693019/1‑150 (MQ=255)       cGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCAGGGCCGTGc                                                                                                                                        <  1:314370/150‑1 (MQ=255)        ggAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGTTCGCCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCAGGGCCGTGcc                                                                                                                                       <  2:349149/150‑2 (MQ=255)        ggAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTTCa                                                                                                                                       >  2:250909/1‑150 (MQ=255)         gAAAAGGTCGGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGGGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCAGGGCCGTGcct                                                                                                                                      <  1:1092834/150‑3 (MQ=255)            aaGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCCTGCCGCGGTAAACCTCCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTAt                                                                                                                                   <  1:812226/150‑1 (MQ=255)              ggTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTAtct                                                                                                                                 <  1:389721/150‑1 (MQ=255)                  cGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTac                                                                                                                             <  2:833233/150‑1 (MQ=255)                    cAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTacac                                                                                                                           <  1:1338209/150‑1 (MQ=255)                      aGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTg                                                                                                                         <  2:1246625/150‑1 (MQ=255)                         tCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGcc                                                                                                                      >  1:580150/1‑150 (MQ=255)                            ccGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGa                                                                                                                   <  1:1307247/150‑1 (MQ=255)                              ggCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGAcc                                                                                                                 >  1:705037/1‑150 (MQ=255)                                     aaGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCGCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCg                                                                                                          <  2:1228950/150‑1 (MQ=255)                                          gagaACGCCGGCCTTGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTcc                                                                                                     >  2:1422001/1‑150 (MQ=255)                                            gaACGCCGGCCTGGTCCTGCTGACAGTGTGAATCGCCGCCATCGCCAGCCTTGCCGCGGGCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTg                                                                                                   >  2:1085889/1‑150 (MQ=255)                                                      cTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATc                                                                                         >  2:391406/1‑150 (MQ=255)                                                      cTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATc                                                                                         >  2:913591/1‑150 (MQ=255)                                                         gTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGc                                                                                      <  1:114979/150‑1 (MQ=255)                                                         gTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGCCCGTGCCCTATCTCTACTCTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGc                                                                                      <  2:565375/150‑1 (MQ=255)                                                         gTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCATGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGc                                                                                      >  2:114979/1‑150 (MQ=255)                                                          tCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCt                                                                                     >  2:670711/1‑150 (MQ=255)                                                          tCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCt                                                                                     <  2:693019/150‑1 (MQ=255)                                                              gctgACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCat                                                                                 >  2:347970/1‑150 (MQ=255)                                                                 gACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCg                                                                                                          >  1:999061/1‑122 (MQ=255)                                                                 gACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCg                                                                                                          <  2:999061/122‑1 (MQ=255)                                                                  aCCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGATCATGGGCAGCATATTc                                                                             >  2:748083/1‑150 (MQ=255)                                                                      tgtgCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCAccc                                                                         <  1:347970/150‑1 (MQ=255)                                                                        tgCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTg                                                                       <  1:670711/150‑1 (MQ=255)                                                                          cATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTTCa                                                                     <  1:1422001/150‑1 (MQ=255)                                                                               ccgccATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTg                                                                                                                         <  1:1053573/93‑1 (MQ=255)                                                                               ccgccATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTg                                                                                                                         >  2:1053573/1‑93 (MQ=255)                                                                               ccgccATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTTGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTAcg                                                                <  2:219566/150‑1 (MQ=255)                                                                                     tCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTAcgcccgcc                                                          >  1:1309194/1‑150 (MQ=255)                                                                                     tCGCCAGCCTTGCCGCGGTAACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCGCCCTGCACTAcgcccgcc                                                          <  2:1108415/150‑1 (MQ=255)                                                                                         cAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGtt                                                      >  1:975256/1‑150 (MQ=255)                                                                                              ttGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTaca                                                 >  1:787828/1‑150 (MQ=255)                                                                                               tGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTacac                                                <  1:1263274/150‑1 (MQ=255)                                                                                                  cgcgGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCt                                                                                                    <  2:994751/95‑1 (MQ=255)                                                                                                  cgcgGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCt                                                                                                    >  1:994751/1‑95 (MQ=255)                                                                                                  cgcgGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCgg                                             <  1:1323454/150‑1 (MQ=255)                                                                                                  cgcgGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCgg                                             >  2:250996/1‑150 (MQ=255)                                                                                                           ccTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGgctgc                                                                                              <  1:92156/92‑1 (MQ=255)                                                                                                           ccTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGgctgc                                                                                              >  2:92156/1‑92 (MQ=255)                                                                                                           ccTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCa                                    <  2:153246/150‑1 (MQ=255)                                                                                                           ccTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCa                                    >  1:1361645/1‑150 (MQ=255)                                                                                                           ccTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCa                                    <  1:250996/150‑1 (MQ=255)                                                                                                               cAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGaa                                <  2:787828/150‑1 (MQ=255)                                                                                                                        tcCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTg                       >  2:484714/1‑150 (MQ=255)                                                                                                                             ccGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCtt                  <  1:268065/150‑1 (MQ=255)                                                                                                                                     tGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTTACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAAcccc                            >  1:1306299/1‑132 (MQ=255)                                                                                                                                     tGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTTACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAAcccc                            <  2:1306299/132‑1 (MQ=255)                                                                                                                                     tGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTTTCGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATg          <  1:1085889/150‑1 (MQ=255)                                                                                                                                     tGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATg          <  2:1050393/150‑1 (MQ=255)                                                                                                                                     tGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATg          >  1:1050393/1‑150 (MQ=255)                                                                                                                                     tGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATg          >  1:1254762/1‑150 (MQ=255)                                                                                                                                     tGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATg          <  1:250909/150‑1 (MQ=255)                                                                                                                                     tGCAAGGCAAGCCCATGGCCGTGAACTCTCTCTAGACTGGGCTGATCGTCGCCGGTTCCTTGCTGCTGTTCGGCGGCATATTCACCCTGCACTACGCCCGCCTCTTCTACACCTGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATg          >  2:15081/1‑150 (MQ=255)                                                                                                                                       tagtgCACGGCCAGGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGc        >  2:1353480/5‑150 (MQ=255)                                                                                                                                          ggCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGCGAg     >  2:1311763/1‑150 (MQ=255)                                                                                                                                           gCACGGCCAGGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGCGAgc    >  2:909194/1‑150 (MQ=255)                                                                                                                                            cACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCCGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGCGAgcg   >  2:486097/1‑150 (MQ=255)                                                                                                                                             aCGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGCGAgcgc  >  1:863847/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                               CCAGCCCGGAAAAGGTCCGCAAGGTCGCCGGCGTCGAAGACGAGAACGCCGGCCTGGTCCTGCTGACCGTGTGCATCGCCGCCATCGCCAGCCTTGCCGCGGTCACCCTGCAACTGGTCTCCAGCCGCGGCCTGCAAGGCACGGCCCTGGCCGTGCACTATCTCTACACTGGCCTGACCGTCGCCGGTTCCTGGCTGCTGATCGGCTGCATATTCACCCTGCACTACGCCCGCCTGTTCTACACCGGGCAAAACCACGAACCCCCACTGCGCTTCGCTGATGGCGAGCGC  >  NC_002947/4061288‑4061577

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