Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	196,759	T→C	19.3%	T755T (ACT→ACC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	196,759	0	T	C	19.3%	109.7 / 10.5	54	T755T (ACT→ACC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base C (5/4);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGG  >  AE015451/196535‑196980                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              accGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACtt                                                                                                                                                                                                                    <  1:17131/236‑1 (MQ=33)      gACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGaa                                                                                                                                                                                                                >  1:526553/1‑235 (MQ=2)                                        ggcgTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCg                                                                                                                                                                                              <  2:220045/219‑1 (MQ=9)                                                   cTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGaa                                                                                                                                                                                                                <  1:146736/190‑1 (MQ=35)                                                              cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAGACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGaa                                                                                                                                                                                                                >  2:211938/1‑179 (MQ=21)                                                              cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAATTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAAc                                                                                                                                                                      >  1:1082344/1‑221 (MQ=14)                                                                    gCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAGACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCa                                                                                                                                                                                                                        >  1:1338542/1‑165 (MQ=14)                                                                    gCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCa                                                                                                                                                                                   <  1:1020322/202‑1 (MQ=255)                                                                                     ccGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCg                                                                                                                                  >  1:552145/1‑234 (MQ=40)                                                                                       ggTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGcc                                                                                                                              <  1:977425/236‑1 (MQ=40)                                                                                          gACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAccc                                                                                                                                          >  1:115336/1‑221 (MQ=40)                                                                                             aGTGACCCTGTCCAACGGCCAGACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACt                                                                                                                                                                  <  2:1024078/194‑1 (MQ=16)                                                                                                gACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGc                                                                                                                                                        <  2:939274/201‑1 (MQ=255)                                                                                                gACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGc                                                                                                                                                        >  1:939276/1‑201 (MQ=255)                                                                                                  cccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAAATTCGAACAGCTGACCCCGAACCCGACCCCGGCTCAGACCACg                                                                                                                                                                            >  2:606958/1‑179 (MQ=16)                                                                                                        ccAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACATCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                                                                                                       <  1:822447/226‑1 (MQ=255)                                                                                                        ccAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                                                                                                       >  2:822445/1‑226 (MQ=255)                                                                                                        ccAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGa                                                                                                                 >  1:1384080/1‑232 (MQ=40)                                                                                                        ccAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACACCGAACCCGTCGCCGGCACAGACCACTGTCTACGAGACGCTGGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAg                                                                                                               <  2:526551/234‑1 (MQ=11)                                                                                                          aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAACTTCGaa                                                                                                                                                                                                                >  1:1296579/1‑135 (MQ=11)                                                                                                             ggCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                          <  2:1384077/234‑1 (MQ=255)                                                                                                                    ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAAATTCGAACAACTGACCCCGAACCCGa                                                                                                                                                                                             >  2:456093/1‑144 (MQ=255)                                                                                                                                aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                          <  1:1503079/215‑1 (MQ=35)                                                                                                                                aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                          >  2:1503076/1‑215 (MQ=35)                                                                                                                                aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                                  >  1:1108722/1‑223 (MQ=35)                                                                                                                                aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCGACAACAACGGTTCGACCGTCAGGGTCACCAACGAGAGCGCCACGGGCGGCaa                                                                                                                                                                                                                       >  2:1248321/1‑106 (MQ=11)                                                                                                                                  gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAACTTCGAAc                                                                                                                                                                                                               >  2:303316/1‑112 (MQ=21)                                                                                                                                   ccGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGa                                                                                                                                                                                                   <  2:414792/123‑1 (MQ=17)                                                                                                                                                       gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGc                                                               <  1:1382596/235‑1 (MQ=34)                                                                                                                                                       gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAGCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                              <  1:1350846/236‑1 (MQ=14)                                                                                                                                                           aTTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                                  <  2:552143/196‑1 (MQ=33)                                                                                                                                                              tCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACa                                                       >  1:502355/1‑236 (MQ=0)                                                                                                                                                                           gAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggc                                                                                                             >  2:831200/1‑169 (MQ=32)                                                                                                                                                                             gtGACGGCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACt                                                                                                                                                                  <  2:814351/113‑1 (MQ=17)                                                                                                                                                                                   tCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCg                                                                                                                                                                                        <  2:1209229/86‑1 (MQ=18)                                                                                                                                                                                             cGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGcc                                                                                                                              <  1:198297/134‑1 (MQ=32)                                                                                                                                                                                             cGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGcc                                                                                                                              >  2:198296/1‑134 (MQ=32)                                                                                                                                                                                                ttCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGccaccc                                                                                                                              >  2:1075778/1‑126 (MQ=17)                                                                                                                                                                                                ttCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                        <  1:1199720/217‑1 (MQ=38)                                                                                                                                                                                                ttCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                        >  2:1199717/1‑217 (MQ=37)                                                                                                                                                                                                  cGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                                                                                            >  1:1192830/1‑147 (MQ=32)                                                                                                                                                                                                  cGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                       >  1:1247493/1‑216 (MQ=37)                                                                                                                                                                                                  cGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGAc                                                                                      >  1:916839/1‑169 (MQ=17)                                                                                                                                                                                                  cGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                                                                                            <  2:1192827/147‑1 (MQ=17)                                                                                                                                                                                                       gTCAGCGTCCCCATCGAGAACGCCACAGGCGGTAACTGCGAACAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCa                                                                    <  2:202619/182‑1 (MQ=255)                                                                                                                                                                                                                 ccATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAc                                                                                                                                            <  2:751395/100‑1 (MQ=18)                                                                                                                                                                                                                 ccATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTg                                                    <  2:1247490/188‑1 (MQ=38)                                                                                                                                                                                                                  cATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                                                                                    <  2:322185/139‑1 (MQ=17)                                                                                                                                                                                                                  cATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAAc                                <  1:487647/207‑1 (MQ=35)                                                                                                                                                                                                                    tCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAggg  >  1:1307760/1‑235 (MQ=37)                                                                                                                                                                                                                          aCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTacac                                                                                            >  1:1213416/1‑139 (MQ=33)                                                                                                                                                                                                                              cACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCg                                                                                                                  >  2:76625/1‑113 (MQ=32)                                                                                                                                                                                                                              cACTGGCGGCAACTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCg                                                                                                                  <  1:76625/113‑1 (MQ=17)                                                                                                                                                                                                                              cACTGGCGGCAACTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                          <  1:916817/121‑1 (MQ=17)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              ACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGG  >  AE015451/196535‑196980

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑