Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	195,406	C→T	12.1%	N304N (AAC→AAT)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	195,406	0	C	T	12.1%	180.8 / 12.8	65	N304N (AAC→AAT)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (34/23);  new base T (3/5);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AGTGGTAGCGTCAACGTCGCAGTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCA  >  AE015451/195173‑195632                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                   aGTGGTAGCGTCAACGTCGCAGTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAAcc                                                                                                                                                                                                                                   >  1:674963/1‑235 (MQ=255)      tAGCGTCAACGTCGCAGTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGGGACCCTGAGCAATCCTGccc                                                                                                                                                                                                                             >  1:289589/1‑236 (MQ=255)            cAACGTCGCAGTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACg                                                                                                                                                                                                                        <  1:593290/235‑1 (MQ=255)              aCGTCGCAGTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGc                                                                                                                                                                                                                               >  1:284108/1‑226 (MQ=255)                    cAGTTCCGGCAAACGACACACCTTACACCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCt                                                                                                                                                                                                                                 <  1:1210910/218‑1 (MQ=255)                      gTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccg                                                                                                                                                                                                              <  1:657840/235‑1 (MQ=255)                             cAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTACCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtaaccgtgac                                                                                                                                                                                                          >  2:670672/1‑232 (MQ=255)                                   acacacCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                <  1:728597/236‑1 (MQ=255)                                       acCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                <  2:986841/232‑1 (MQ=255)                                           tACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGt                                                                                                                                                                                        <  1:604638/236‑1 (MQ=255)                                               tCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                      >  2:1129747/1‑202 (MQ=255)                                                      tGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTg                                                                                                                                                                             >  1:554684/1‑236 (MQ=255)                                                           aGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAcccc                                                                                                                                                                                                       >  1:588326/1‑204 (MQ=255)                                                           aGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGa                                                                                                                                                                            <  2:1156689/232‑1 (MQ=255)                                                             aTTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                       <  1:28775/235‑1 (MQ=255)                                                             aTTCCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGAAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                <  1:368791/210‑1 (MQ=255)                                                              tttCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                       <  2:1181852/234‑1 (MQ=255)                                                                                           ggtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAc                                                                                                                                        >  1:87229/1‑236 (MQ=255)                                                                                             tggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCCGCCCAGACGcc                                                                                                                                                                                                                      >  1:975620/1‑158 (MQ=255)                                                                                             tggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCa                                                                                                                                                                                                  >  1:733013/1‑176 (MQ=255)                                                                                              ggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCCGCCCac                                                                                                                                                                                                                           >  1:339619/1‑149 (MQ=255)                                                                                                   ccTGACGCTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCATGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCaaac                                                                                                                                                                                        <  2:654089/179‑4 (MQ=18)                                                                                                       aCGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCggag                                                                                                                                                                                                                  >  1:189290/1‑148 (MQ=255)                                                                                                            gACATCGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                      <  2:588321/141‑1 (MQ=255)                                                                                                                 tGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                   <  2:1213127/219‑1 (MQ=255)                                                                                                                  gCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCCACGGCCAgg                                                                                                                                                                                         >  1:292422/1‑164 (MQ=255)                                                                                                                   ccGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACCCTACAACCGCGCCCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGc                                                                                                                                                                         <  1:405148/179‑1 (MQ=255)                                                                                                                           gCCAGCGGTTCCCCAAGTGCCCGATACAATCG‑‑CCACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCt                                                                                                                                                                                                       <  1:283088/141‑1 (MQ=17)                                                                                                                                 ggTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGCGCCCGTGACCCTGTCCCACGGCCAg                                                                                                                                                                                          >  2:399766/1‑148 (MQ=255)                                                                                                                                        cAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGGAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccgtg                                                                                                                                                                                                            <  2:690081/123‑1 (MQ=255)                                                                                                                                        cAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGt                                                                                                                                                                                                     >  1:1021142/1‑130 (MQ=255)                                                                                                                                             gACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGa                                                                                                                                                                                                                          >  1:1122383/1‑104 (MQ=255)                                                                                                                                             gACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGc                                                                                                                                                                                                                       <  2:1021137/107‑1 (MQ=255)                                                                                                                                              aCCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCa                                                                                                  >  1:94603/1‑223 (MQ=11)                                                                                                                                                   tACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGcc                                                                                                                                                                                                                              >  2:524015/1‑94 (MQ=255)                                                                                                                                                    aCAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGc                                                                               >  1:261451/1‑236 (MQ=255)                                                                                                                                                      aaTCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCa                                                                                                                                                                                                                            <  2:320250/93‑1 (MQ=255)                                                                                                                                                        tCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                           <  1:300291/188‑1 (MQ=255)                                                                                                                                                                  aCAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt                                                                                                                         >  2:10009/1‑182 (MQ=255)                                                                                                                                                                  aCAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACCGTTGAAGcc                                                                                                                                                                        >  2:1031071/1‑135 (MQ=255)                                                                                                                                                                  aCAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGCGACCGTGACCCTGCCCAACGGCCAGGTCATCACTGTTGCAGCCGGCAAAACCCAGGGCAGCGTGGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTGCGACCGTCAGc                                                                                                >  1:423528/1‑207 (MQ=255)                                                                                                                                                                    aaCCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACCGTTGAAGCGGGCAAAACCCCGGGCac                                                                                                                                                       >  1:1101423/1‑149 (MQ=255)                                                                                                                                                                       cgcgACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGt                                                                                                    <  2:485633/198‑1 (MQ=255)                                                                                                                                                                          gACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                           <  1:258814/172‑1 (MQ=255)                                                                                                                                                                              cTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                      >  1:694147/1‑77 (MQ=39)                                                                                                                                                                              cTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATCAcc                                                                                                                                                                                 >  1:482705/1‑113 (MQ=255)                                                                                                                                                                                    gCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                       <  2:758186/118‑1 (MQ=255)                                                                                                                                                                                    gCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                       >  1:758191/1‑118 (MQ=255)                                                                                                                                                                                         gCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACg                                                                                                                                                                                               >  1:1112927/1‑89 (MQ=255)                                                                                                                                                                                                          gTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAcc                           <  1:1129752/236‑1 (MQ=255)                                                                                                                                                                                                           tGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                          <  1:344160/220‑1 (MQ=255)                                                                                                                                                                                                             gCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                      <  1:1254512/46‑1 (MQ=38)                                                                                                                                                                                                                   tCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGA‑CCCGGCTCAGACCACGATCAACg                 <  1:1078025/236‑1 (MQ=255)                                                                                                                                                                                                                     aCCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGc                                                                               >  1:528451/1‑173 (MQ=33)                                                                                                                                                                                                                               tGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGc                                                           >  1:149671/1‑183 (MQ=38)                                                                                                                                                                                                                                gACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAATGATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCg        >  1:108612/1‑233 (MQ=255)                                                                                                                                                                                                                                  cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                >  1:1056095/1‑47 (MQ=255)                                                                                                                                                                                                                                  cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCATGTCTTAACCATTGATGCGGGCCCACCCCAGGGCAGCGTCGATTGCCAGACCCCGTCCAATTACGt                                                                                                                         >  1:1051131/1‑118 (MQ=14)                                                                                                                                                                                                                                  cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGCCCAGAtcctca                                                                                                                                                                                   >  1:342999/1‑60 (MQ=14)                                                                                                                                                                                                                                   ccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGattcca  >  1:894033/1‑231 (MQ=255)                                                                                                                                                                                                                                    cTGAGCAACCCTGCCCAGACGCCGgtga                                                                                                                                                                                                                 >  2:545896/1‑28 (MQ=255)                                                                                                                                                                                                                                    cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTg                                                                                                                                                                                                      >  1:1182145/1‑39 (MQ=255)                                                                                                                                                                                                                                      gAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg            <  2:894028/223‑1 (MQ=255)                                                                                                                                                                                                                                        gCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACgggcg                                                                       >  1:903096/1‑162 (MQ=21)                                                                                                                                                                                                                                          aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTc         >  1:1158345/1‑222 (MQ=255)                                                                                                                                                                                                                                           aCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACtggcggca                                                                    <  1:564760/162‑3 (MQ=9)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                   AGTGGTAGCGTCAACGTCGCAGTTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCG‑‑ACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCA  >  AE015451/195173‑195632

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