Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,653,239	A→C	6.4%	K12Q (AAG→CAG)	PP_5659 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,653,239	0	A	C	6.4%	223.9 / 6.7	79	K12Q (AAG→CAG)	PP_5659	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base C (2/3);  total (38/41)
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.43e-01

CAAGGCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAAT  >  NC_002947/4653091‑4653377                                                                                                                                                     |                                                                                                                                           cacggCTTCCACCAAAGGCCTTGAACCGTACTGCCAAGTCGCATTAGACCAAACATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACCGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCaa                                                                                                                                           <  2:696014/147‑1 (MQ=255)  aaGACTTCGAGCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGGTGTAGTAAGGACAAAAAAGCACTACAATACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCTCTAATCAGGGGCCCTGGCAAg                                                                                                                                          <  2:173080/150‑1 (MQ=255)   aGGCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAgg                                                                                                                                         <  1:906237/150‑1 (MQ=255)    ggCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAggg                                                                                                                                        <  1:1918058/150‑1 (MQ=255)     gCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGc                                                                                                                                       <  2:1647286/150‑1 (MQ=255)       ttCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTTGCCAGGGccc                                                                                                                                     <  1:376221/150‑2 (MQ=255)       ttCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCt                                                                                                                                     >  2:1863871/1‑150 (MQ=255)         cGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTc                                                                                                                                   <  2:113456/150‑1 (MQ=255)           aCCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCgg                                                                                                                                 >  2:1872009/1‑150 (MQ=255)            ccAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCggg                                                                                                                                <  2:796782/150‑1 (MQ=255)            ccAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCggg                                                                                                                                <  2:1949603/150‑1 (MQ=255)             cAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCgggg                                                                                                                               >  2:964239/1‑150 (MQ=255)                 ggCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTCCAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGgcgcg                                                                                                                           <  1:1484061/150‑1 (MQ=255)                 ggCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGgcgcg                                                                                                                           >  2:147908/1‑150 (MQ=255)                    cTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTc                                                                                                                        >  1:1020899/1‑150 (MQ=255)                     ttGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCCCCCGCTCATCAGGCGTCCTGGCAAGGGCCGTCGGGGCGCGCTcc                                                                                                                       >  2:146103/1‑150 (MQ=255)                       gAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCtt                                                                                                                     <  1:346236/150‑1 (MQ=255)                               ctgCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCa                                                                                                             <  2:1565165/150‑1 (MQ=255)                                tgCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCag                                                                                                            >  2:2361762/1‑150 (MQ=255)                                  ccAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAg                                                                                                                                          >  2:692020/1‑118 (MQ=255)                                  ccAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAg                                                                                                                                          <  1:692020/118‑1 (MQ=255)                                  ccAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTACCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCagga                                                                                                          >  2:826400/1‑150 (MQ=255)                                    aaGTCGCATTAGTCCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACGTTTACTTCCATCATATTATCGTTGTTACGATGGCCTCGCTCCCGCTACTCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCaggagg                                                                                                        <  1:768249/150‑1 (MQ=255)                                         gCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCt                                                                                                   >  2:1558169/1‑150 (MQ=255)                                         gCATTAGACCAAAGACGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGACTTATTCCCAAGCGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCt                                                                                                   >  1:381158/1‑150 (MQ=255)                                           aTTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAANACATTTNACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCg                                                                                                 >  1:2096386/1‑150 (MQ=255)                                            ttAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCgt                                                                                                <  2:2418238/150‑1 (MQ=255)                                                  cAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAg                                                                                                                                          >  1:97187/1‑102 (MQ=255)                                                  cAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAg                                                                                                                                          <  2:97187/102‑1 (MQ=255)                                                  cAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTc                                                                                                  >  1:659114/1‑142 (MQ=255)                                                  cAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTc                                                                                                  <  2:659114/142‑1 (MQ=255)                                                     aGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACg                                                                                       <  1:198220/150‑1 (MQ=255)                                                     aGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGCTTACTCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACg                                                                                       <  1:147908/150‑1 (MQ=255)                                                         gtaggaAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAAAGGCCCTGACGAAGTATTGCAGCACACACTTTTGCCACCCTAGGAGGCGCCTCGTGTCTATACGGccc                                                                                   <  1:1863871/150‑1 (MQ=255)                                                            gtaAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCcaca                                                                                >  2:77569/1‑150 (MQ=255)                                                                ggACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGc                                                                            >  1:597309/1‑150 (MQ=255)                                                                    aaaaaaGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCa                                                                                                                  <  2:182768/108‑1 (MQ=255)                                                                    aaaaaaGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCa                                                                                                                  >  1:182768/1‑108 (MQ=255)                                                                    aaaaaaGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGtt                                                                        <  1:1558169/150‑1 (MQ=255)                                                                    aaaaaaGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGtt                                                                        >  1:911221/1‑150 (MQ=255)                                                                         aGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACg                                                                   <  1:2279454/150‑1 (MQ=255)                                                                              tACAAAACATTTGACAGGGTCTTTTTTTTTTTTTCATGATGGCATCGCTCCCGCTAATCAGGTGCCCTGGCAAGGGCCTTCGGTGCGCGCTCCTTCCACCCCAGTACGCGCCTTGTGTCTATTCGGCCCACAAATCGGTTCGACGGGAtt                                                              >  2:1004716/1‑150 (MQ=255)                                                                              tACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTTGCCAGGGcccc                                                                                                                                    <  1:162167/80‑3 (MQ=255)                                                                              tACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTTGCCAGGGcccc                                                                                                                                    >  2:162167/1‑78 (MQ=255)                                                                              tACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGAtt                                                              <  1:77569/150‑1 (MQ=255)                                                                                cAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTgg                                                            >  1:1896099/1‑149 (MQ=255)                                                                                  aaaCATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAg                                                                              >  2:2390671/1‑130 (MQ=255)                                                                                  aaaCATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAg                                                                              <  1:2390671/130‑1 (MQ=255)                                                                                             aGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCggg                                                                                                                                <  2:679725/69‑1 (MQ=255)                                                                                             aGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCggg                                                                                                                                >  1:679725/1‑69 (MQ=255)                                                                                               ggTCTTTTTCCCAAGGGCATGATGGCATCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCTCAAGTCGGTTCGACGGGATTGCTACTGCAACgaagaa                                             >  2:520797/1‑150 (MQ=255)                                                                                                gTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGaagaag                                            >  2:1280447/1‑150 (MQ=255)                                                                                                   tttttCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCgg                                                                          <  2:667886/117‑1 (MQ=255)                                                                                                   tttttCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCgg                                                                          >  1:667886/1‑117 (MQ=255)                                                                                                   tttttCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGtt                                          >  2:130066/1‑149 (MQ=255)                                                                                                   tttttCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGtt                                          <  1:130066/149‑1 (MQ=255)                                                                                                   tttttCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTg                                         >  2:141320/1‑150 (MQ=255)                                                                                                     tttCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTc                                       <  2:1752498/150‑1 (MQ=255)                                                                                                         ccAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGa                                   <  1:964239/150‑1 (MQ=255)                                                                                                          cAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAg                                  <  1:2361762/150‑1 (MQ=255)                                                                                                          cAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAg                                  >  2:285749/1‑150 (MQ=255)                                                                                                          cAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAg                                  <  2:477769/150‑1 (MQ=255)                                                                                                              ggCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGAc                              <  1:1266321/150‑1 (MQ=255)                                                                                                                   gatgGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAg                         <  2:2096386/150‑1 (MQ=255)                                                                                                                        ccTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCAt                    >  1:528988/1‑150 (MQ=255)                                                                                                                         cTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGAc                              >  2:2272306/1‑139 (MQ=255)                                                                                                                         cTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGAc                              <  1:2272306/139‑1 (MQ=255)                                                                                                                         cTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTACAGAGTGAAGGAACCGCATc                   >  2:436677/1‑150 (MQ=255)                                                                                                                            gCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCAcc                <  1:1280447/150‑1 (MQ=255)                                                                                                                               cccGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATAGCGACTGCAACGAAGAAGTTGTCCCTAGGGACGGAAGCGCATCACCGCa             <  1:1872009/150‑1 (MQ=255)                                                                                                                                 cgatAATCAGGGGCCCTGGCAAGGGCCTCCGGGGAAGGCTCCTTCCACCCCAGGAAGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCAGCACCGCAGt           >  1:1143893/4‑150 (MQ=255)                                                                                                                                    tAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTAGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACt        >  2:1385301/1‑150 (MQ=255)                                                                                                                                      cccgaaggCCCTTGCCAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAg                                  <  1:238186/116‑1 (MQ=255)                                                                                                                                      cccgaaggCCCTTGCCAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAg                                  >  2:238186/7‑122 (MQ=255)                                                                                                                                      aTCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTgg      >  1:2309591/1‑150 (MQ=255)                                                                                                                                       tCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATAc                                                                                        >  1:168772/1‑67 (MQ=255)                                                                                                                                       tCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATAc                                                                                        <  2:168772/67‑1 (MQ=255)                                                                                                                                       tCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGc     >  1:1850662/1‑150 (MQ=255)                                                                                                                                          ggggCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAAt  <  2:2309591/150‑1 (MQ=255)                                                                                                                                          ggggCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAAt  <  2:1850662/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                           CAAGGCTTCGACCAAAGGCCTTGAAACTGACTGCCAAGTCGCATTAGACCAAAGATGTAGTAAGGACAAAAAAGCACTACAAAACATTTGACAGGGTCTTTTTCCCAAGGGCATGATGGCCTCGCTCCCGCTAATCAGGGGCCCTGGCAAGGGCCTTCGGGGCGCGCTCCTTCCACCCCAGGAGGCGCCTCGTGTCTATACGGCCCACAAGGCGGTTCGACGGGATTGCGACTGCAACGAAGAAGTTGTCCCGAGGGACGGAAGCGCATCACCGCAGTACTGGCAAT  >  NC_002947/4653091‑4653377

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