Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 216,994 C→T 9.5% S7500S (AGC→AGT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451216,9940CT 9.5% 156.1 / 5.0 52S7500S (AGC→AGTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base T (2/3);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

AACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAACCAGGTTCGGGAACTCCAGGCACCGGC  >  AE015451/216779‑217216
                                                                                                                                                                                                                       |                                                                                                                                                                                                                              
aaCTCGGCGGTGGCCGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                                                                                     <  1:253812/227‑1 (MQ=255)
     ggcggTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTg                                                                                                                                                                                                        <  1:577826/235‑1 (MQ=255)
                              ttCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTg                                                                                                                                                                                      >  2:679050/1‑228 (MQ=255)
                                    aaCCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCATGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                                                         <  2:426967/235‑1 (MQ=255)
                                      ccTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                                                         <  2:287422/233‑1 (MQ=255)
                                        gcgcgCTGGGCGGTGCCGCGAAGGTGGATGTCCCCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGtca                                                                                                                                                                                                                <  2:254238/188‑1 (MQ=255)
                                                ggCGATGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCa                                                                                                                                                                                                 <  1:802005/199‑1 (MQ=255)
                                                         gcgAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTa                                                                                                                                                   <  1:774959/236‑1 (MQ=255)
                                                            aaGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGtca                                                                                                                                                                                                                <  1:668273/172‑1 (MQ=255)
                                                            aaGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACg                                                                                                                                                 <  1:730742/235‑1 (MQ=255)
                                                                    agTCACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGcc                                                                                                                                                                 <  1:317565/210‑1 (MQ=255)
                                                                            acacCACCGACGAAGTGGTGGCCAAGCTGACAGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACAGTCATCACCGTg                                                                                                                                                                                                        >  2:381365/1‑164 (MQ=255)
                                                                                       gaAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAat                                                                                                                      >  1:535523/1‑235 (MQ=40)
                                                                                         aGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAatcga                                                                                                                   <  1:400944/236‑1 (MQ=40)
                                                                                               ggCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCgg                                                                                                                                                             <  1:679052/188‑1 (MQ=39)
                                                                                                   aaGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGt                                                                                                         >  1:710790/1‑236 (MQ=255)
                                                                                                   aaGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGt                                                                                                         >  1:764391/1‑236 (MQ=39)
                                                                                                          ccGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCCCCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGGGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTc                                                                                                                                               <  2:264569/191‑1 (MQ=255)
                                                                                                            gcgACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAAc                                                                                                            >  1:526202/1‑224 (MQ=255)
                                                                                                                  ccTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGt                                                                                          >  1:231108/1‑236 (MQ=255)
                                                                                                                                     gcgAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTg                                                                       <  1:502642/236‑1 (MQ=255)
                                                                                                                                       gagaTCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTg                                                                       <  2:764389/234‑1 (MQ=255)
                                                                                                                                        agaTCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCaaa                                                                                                                       >  1:427447/1‑185 (MQ=38)
                                                                                                                                               cTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTg                                                                       <  2:784199/226‑1 (MQ=255)
                                                                                                                                                      aTCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGAcccc                                                      <  1:721939/236‑1 (MQ=255)
                                                                                                                                                        cACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACGGCTGACACTGGACAAGAccc                                                       >  1:558869/1‑233 (MQ=255)
                                                                                                                                                          cccTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACag                                                            >  2:426746/1‑225 (MQ=255)
                                                                                                                                                          cccTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGt                                                   >  1:539463/1‑235 (MQ=255)
                                                                                                                                                           ccTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTg                                                                       <  1:70330/214‑1 (MQ=255)
                                                                                                                                                                  aacaaAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTggcgg                                                                                                   >  2:766250/1‑179 (MQ=35)
                                                                                                                                                                       aGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTg                                      >  2:673737/1‑235 (MQ=255)
                                                                                                                                                                       aGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGCCCAACGGCACCGTGGGCACTGCTACCGTCACAGCCCCGGACAACGTCTAc                                                                                                                                                  >  2:219513/1‑127 (MQ=17)
                                                                                                                                                                             ccTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGtt                                                                                         <  2:535920/177‑1 (MQ=25)
                                                                                                                                                                              cTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGtt                                                                                                        >  1:566269/1‑162 (MQ=31)
                                                                                                                                                                                 ccAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGCTGTTGGCACGTTCGaa                                                                             >  2:62392/1‑186 (MQ=255)
                                                                                                                                                                                        acaacCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAACCAGGTTCg                    <  1:749478/236‑1 (MQ=255)
                                                                                                                                                                                                 gTGCATTGACCTTCACGCTGAGTGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGtt                                                                                 <  2:561663/166‑1 (MQ=255)
                                                                                                                                                                                                  tGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTc                                                                                >  2:159817/1‑166 (MQ=40)
                                                                                                                                                                                                  tGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTc                                                                                <  1:159819/166‑1 (MQ=40)
                                                                                                                                                                                                    cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGACCAGCTGACACTGGACAAGACCCCGGt                                                   >  2:210272/1‑193 (MQ=255)
                                                                                                                                                                                                    cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAACCAGGTTCGGGAACTCCAgg         <  1:541110/235‑1 (MQ=255)
                                                                                                                                                                                                     aTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAACCAGGTTCGGGAACTCCAg          >  2:724512/1‑233 (MQ=255)
                                                                                                                                                                                                       tGACCTTCACGCTGAGTGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTgg                                                                <  1:303124/177‑1 (MQ=255)
                                                                                                                                                                                                       tGACCTTCACGCTGAGTGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTgg                                                                >  2:303122/1‑177 (MQ=255)
                                                                                                                                                                                                          ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAACCAGGTTCGGGAACTCCAGGCACCGGc  >  1:30946/1‑236 (MQ=255)
                                                                                                                                                                                                            ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTg                                                                       <  2:741106/165‑1 (MQ=255)
                                                                                                                                                                                                            ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTg                                                                       >  1:741108/1‑165 (MQ=255)
                                                                                                                                                                                                               aCGCTGAGTGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTggcg                                                                                                    >  2:520510/1‑133 (MQ=17)
                                                                                                                                                                                                               aCGCTGAGTGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTggcg                                                                                                    <  1:520512/133‑1 (MQ=17)
                                                                                                                                                                                                               aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGt                                       >  2:460541/1‑194 (MQ=255)
                                                                                                                                                                                                                 gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAAc                                                                            >  2:102407/1‑155 (MQ=255)
                                                                                                                                                                                                                 gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACCGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAg                                                                          >  2:626464/1‑157 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                              
AACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAACCAGGTTCGGGAACTCCAGGCACCGGC  >  AE015451/216779‑217216

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