Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,522,499	C→G	5.5%	A116A (GCC→GCG)	idh →	isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,522,499	0	C	G	5.5%	238.1 / 3.4	73	A116A (GCC→GCG)	idh	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base G (2/2);  total (39/34)
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.03e-01

GGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGGCGCCT  >  NC_002947/4522355‑4522647                                                                                                                                                 |                                                                                                                                                      ggAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGccgccg                                                                                                                                                   >  1:399175/1‑149 (MQ=255)  gAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACcgc                                                                                                                                                  <  1:248718/149‑1 (MQ=255)          ccAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGgaagaa                                                                                                                                           >  4:229633/1‑148 (MQ=255)             gCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCGCCGCGGAAGAAAAAg                                                                                                                                       <  2:399175/149‑1 (MQ=255)               tGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGaa                                                                                                                                     >  4:167043/1‑149 (MQ=255)               tGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGaa                                                                                                                                     >  1:167120/1‑149 (MQ=255)                 gCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATc                                                                                                                                   <  2:167120/149‑1 (MQ=255)                        ccccTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGaaaaa                                                                                                                                        <  7:118254/137‑1 (MQ=255)                        ccccTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGaaaaa                                                                                                                                        >  8:118254/1‑137 (MQ=255)                        ccccTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGcc                                                                                                                            <  5:387345/149‑1 (MQ=255)                           cTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGc                                                                                                                         >  5:84785/1‑149 (MQ=255)                                gCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGgaa                                                                                                                                              >  1:200640/1‑123 (MQ=255)                                gCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGgaa                                                                                                                                              <  2:200640/123‑1 (MQ=255)                                gCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGa                                                                                                                    >  1:153795/1‑149 (MQ=255)                                 ccAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGAc                                                                                                                   >  7:165819/1‑149 (MQ=255)                                    aCATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGaaaaa                                                                                                                                        <  7:83194/125‑1 (MQ=255)                                    aCATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGaaaaa                                                                                                                                        >  8:83194/1‑125 (MQ=255)                                        catcaAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATc                                                                                                             >  2:42081/1‑148 (MQ=255)                                         atcaAGCTGCCGAACATCAGGGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATc                                                                                                             >  8:279122/1‑147 (MQ=255)                                          tcaAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          <  8:165819/149‑1 (MQ=255)                                          tcaAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          >  3:195841/1‑149 (MQ=255)                                          tcaAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          >  3:203725/1‑149 (MQ=255)                                            aaGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAggg                                                                                                        <  1:356844/149‑1 (MQ=255)                                              gCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCt                                                                                                      >  8:127325/1‑149 (MQ=255)                                                  ccGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCGACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTccgc                                                                                                  >  7:60852/1‑149 (MQ=255)                                                  ccGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCGACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTccgc                                                                                                  >  3:278906/1‑149 (MQ=255)                                                     aaCATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGt                                                                                               >  8:157300/1‑149 (MQ=255)                                                        aTCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGaa                                                                                            >  1:281027/1‑149 (MQ=255)                                                        aTCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGaa                                                                                            >  1:237685/1‑149 (MQ=255)                                                        aTCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACAACCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGaa                                                                                            >  1:328822/1‑149 (MQ=255)                                                             cgcCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCgg                                                                                       <  1:162204/149‑1 (MQ=255)                                                                 tCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCTGTGAACCCGGtgct                                                                                   >  8:146833/1‑149 (MQ=255)                                                                   ggTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          >  3:284357/1‑124 (MQ=255)                                                                   ggTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAg                                                                                                          <  4:284357/124‑1 (MQ=255)                                                                       ccGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAAcc                                                                                          <  6:302555/136‑1 (MQ=255)                                                                       ccGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAAcc                                                                                          >  5:302555/1‑136 (MQ=255)                                                                        cGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCTGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtg                                                                                     <  7:110435/140‑1 (MQ=255)                                                                        cGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCTGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtg                                                                                     >  8:110435/1‑140 (MQ=255)                                                                        cGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGaa                                                                            <  3:167043/149‑1 (MQ=255)                                                                         gCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAg                                                                           <  1:42081/149‑1 (MQ=255)                                                                          cAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtgct                                                                                   <  3:372270/140‑1 (MQ=255)                                                                          cAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtgct                                                                                   >  4:372270/1‑140 (MQ=255)                                                                          cAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAgg                                                                          >  4:116808/1‑149 (MQ=255)                                                                                aaaGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTc                                                                   >  2:285219/1‑149 (MQ=255)                                                                                 aaGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTTCGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTcc                                                                  >  2:200771/1‑149 (MQ=255)                                                                                 aaGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTcc                                                                  <  3:116808/149‑1 (MQ=255)                                                                                  aGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGtg                                                                                     >  8:234274/1‑130 (MQ=255)                                                                                    ccGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGAc                                                               <  7:279122/149‑1 (MQ=255)                                                                                              ggAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑c                                                                        >  4:2461/1‑130 (MQ=255)                                                                                              ggAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑c                                                                        <  3:2461/130‑1 (MQ=255)                                                                                                  cTGCAAGGCAAGGGCTTCAACATCCATGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACAGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑AAACTCCGACAGACGCGCGCAGCt                                                 >  7:151968/1‑149 (MQ=255)                                                                                                   tGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTg                                                >  4:344697/1‑149 (MQ=255)                                                                                                     caaggcaaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTc                                              >  5:246881/1‑149 (MQ=255)                                                                                                         gcaaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTc                                          <  2:153795/149‑1 (MQ=255)                                                                                                         gcaaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTc                                          <  6:246881/149‑1 (MQ=255)                                                                                                           aaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTc                                          >  3:179850/1‑147 (MQ=255)                                                                                                           aaggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCaa                                        <  4:203725/149‑1 (MQ=255)                                                                                                            aggGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAg                                       <  2:328822/149‑1 (MQ=255)                                                                                                             gggCTTCAACATCCCTGACTACGCCGACGAGCCGGCGACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACAGCGGCGCGCGGCGGTCggagttgc                                      <  8:60852/149‑7 (MQ=255)                                                                                                             gggCTTCAACATCCCTGACTACGCCGACGAGCCGGCGACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACAGCGGCGCGCGGCGGTCggagttgc                                      <  4:278906/149‑7 (MQ=255)                                                                                                                   cAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACg                                >  4:234157/1‑149 (MQ=255)                                                                                                                       aTCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGGCCAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACTccc                             <  1:285219/149‑1 (MQ=255)                                                                                                                          ccTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCaa                         <  3:229633/149‑1 (MQ=255)                                                                                                                           cTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAg                        <  7:127325/149‑1 (MQ=255)                                                                                                                                  cGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCAccc                   <  2:237685/147‑1 (MQ=255)                                                                                                                                   gCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCa                      <  4:81005/143‑1 (MQ=255)                                                                                                                                   gCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCa                      >  3:81005/1‑143 (MQ=255)                                                                                                                                   gCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCcaca                >  3:397112/1‑149 (MQ=255)                                                                                                                                   gCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCcac                 >  3:390835/1‑148 (MQ=255)                                                                                                                                      gacgaGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACaa             <  7:157300/149‑1 (MQ=255)                                                                                                                                       acgaGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCcacac               >  3:280720/1‑146 (MQ=255)                                                                                                                                           gCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATgg        <  6:84785/149‑1 (MQ=255)                                                                                                                                               gCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGgtgc    <  4:397112/149‑1 (MQ=255)                                                                                                                                               gCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGgtgc    <  4:390835/149‑1 (MQ=255)                                                                                                                                                 cACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGGCGCCt  >  5:261575/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAACATCATCAAGCTGCCGAACATCAGCGCCTCGGTACCGCAGCTCAAAGCCGCGATCAAGGAACTGCAAGGCAAGGGCTTCAACATCCCTGACTACGCCGACGAGCCGGCCACCGCGGAAGAAAAAGAATCCCGCGCCCGCTACGACCGCATCAAGGGCTCCGCCGTGAACCCGGTGCTGCGCGAAGG‑CAACTCCGACCGCCGCGCGCCGCTGTCGGTCAAGAACTACGCCCGCAAGCACCCACACAAGATGGGCGCCT  >  NC_002947/4522355‑4522647

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