Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 213,874 T→G 11.6% T6460T (ACT→ACG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451213,8740TG11.6% 153.8 / 3.1 51T6460T (ACT→ACGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (16/27);  new base G (3/4);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAT  >  AE015451/213651‑214099
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACcgac                                                                                                                                                                                                                       <  1:159654/236‑1 (MQ=255)
            acaacCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTAcc                                                                                                                                                                                                           <  1:440519/236‑1 (MQ=255)
                                ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGAccc                                                                                                                                                                                       <  1:416720/236‑1 (MQ=255)
                                     gCTGGGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGActtggtc                                                                                                                                                                                   <  1:315879/235‑7 (MQ=255)
                                            gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGAcct                                                                                                                                                                                       >  1:28816/1‑223 (MQ=255)
                                                        gTCATCACGCTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                                                                                                                                                            <  2:776678/207‑1 (MQ=255)
                                                                                 ccGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCAcc                                                                                                                                                                                                                          <  2:358813/152‑1 (MQ=255)
                                                                                       gTACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAcc                                                                                                                                 <  1:399919/235‑1 (MQ=255)
                                                                                        tACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAcc                                                                                                                                 <  2:416122/234‑1 (MQ=255)
                                                                                              cACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGCCCAGACCTCGGTGGCCGCGAACGTCACACCGACATTCa                                                                                                                         >  2:176798/1‑236 (MQ=255)
                                                                                                 cGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACt                                                                                                                       >  1:772126/1‑235 (MQ=255)
                                                                                                  gTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCgg                                                                                                                                                                                     >  2:664998/1‑172 (MQ=255)
                                                                                                    gACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGCTCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACCCTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCCCCGCCGACCAGACCTCGGTGGCCGAGAACGTCaa                                                                                                                                    <  1:206734/219‑1 (MQ=255)
                                                                                                     actccgccGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGCCCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGc                                                                                                                                                <  1:269005/200‑1 (MQ=255)
                                                                                                       tGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCAc                                                                                                                 >  1:698167/1‑235 (MQ=255)
                                                                                                                      cTACGTCGGTACCAACGACCCTGTCGTGCTGTCGATCGCTACCGTTGAGGGTGCGGCCGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTc                                                                                                                                                                                   >  2:508978/1‑154 (MQ=255)
                                                                                                                      cTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGcc                                                                                                                                               <  1:624452/190‑1 (MQ=255)
                                                                                                                      cTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAc                                                                                                                                  >  2:556376/1‑203 (MQ=255)
                                                                                                                      cTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGa                                                                                                                               >  2:690695/1‑206 (MQ=255)
                                                                                                                             ggcaCCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGa                                                                                          <  1:409707/233‑1 (MQ=255)
                                                                                                                                 ccAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCgg                                                                                                                                                    >  1:633438/1‑174 (MQ=12)
                                                                                                                                 ccAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCgg                                                                                                                                                    <  2:633436/174‑1 (MQ=2)
                                                                                                                                  cAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCATGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGAcc                                                                                                                                                        <  2:451906/169‑1 (MQ=255)
                                                                                                                                    aCGACCCTGTCGTGATGTCGATCGCTACCGTTGAGGGTGCGGAAGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTc                                                                                                                                      <  1:443689/185‑1 (MQ=255)
                                                                                                                                        cccTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACcgac                                                                                                                                                                                                                       >  1:595434/1‑100 (MQ=255)
                                                                                                                                        cccTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACcgac                                                                                                                                                                                                                       <  2:595432/100‑1 (MQ=255)
                                                                                                                                               gTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                                                        >  1:280188/1‑236 (MQ=255)
                                                                                                                                                          cGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAgg                                                                                                                                                                                             >  1:235375/1‑108 (MQ=255)
                                                                                                                                                                            ggACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCCCCATCAAgg                                              >  2:423786/1‑233 (MQ=255)
                                                                                                                                                                               cGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACcgac                                                                                                                                                                                                                       <  2:786970/61‑1 (MQ=255)
                                                                                                                                                                               cGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACcgac                                                                                                                                                                                                                       >  1:786972/1‑61 (MQ=255)
                                                                                                                                                                                      aaGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGc                                                                                                       <  2:470820/166‑1 (MQ=255)
                                                                                                                                                                                      aaGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGc                                                                                                       >  1:470822/1‑166 (MQ=255)
                                                                                                                                                                                      aaGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                  <  1:469195/235‑1 (MQ=255)
                                                                                                                                                                                           cGAGCACCGGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGa                                                                                          <  2:153473/174‑1 (MQ=255)
                                                                                                                                                                                           cGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATc                                                  >  1:641672/1‑214 (MQ=255)
                                                                                                                                                                                            gAGCAACTGACACTCGACAAGACCCCGGTCAGCACCTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGcc                           <  1:706104/236‑1 (MQ=255)
                                                                                                                                                                                             aGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGt                                                                                                                                                                                                              >  2:109987/1‑56 (MQ=35)
                                                                                                                                                                                                aaCTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCGGGTACGCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                                                                                <  1:702644/179‑1 (MQ=255)
                                                                                                                                                                                                   tGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                           <  2:756181/197‑1 (MQ=255)
                                                                                                                                                                                                     acacTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGCCCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCa                                  <  1:192720/220‑1 (MQ=255)
                                                                                                                                                                                                         tCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCCCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCacgac              <  1:430782/236‑1 (MQ=255)
                                                                                                                                                                                                           gACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCCTCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                     <  2:442193/163‑1 (MQ=255)
                                                                                                                                                                                                            aCAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCacgac              <  1:697320/233‑1 (MQ=255)
                                                                                                                                                                                                               agcccccGGCCAGCCCTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTc                                                                                                                                      <  2:202855/107‑1 (MQ=255)
                                                                                                                                                                                                                   cccGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGc     <  1:737898/235‑1 (MQ=255)
                                                                                                                                                                                                                   cccGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCg    <  1:690697/236‑1 (MQ=255)
                                                                                                                                                                                                                        tCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGCGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCCCACGCCCTGGTCGTGACCCTGAGCCACAACGc                                                               >  2:206073/1‑172 (MQ=255)
                                                                                                                                                                                                                          aGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCCTCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAt  >  1:42565/1‑231 (MQ=255)
                                                                                                                                                                                                                           gCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCaacaac                                                                 <  2:336712/167‑1 (MQ=255)
                                                                                                                                                                                                                             aCTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCacgac              <  2:581601/216‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
GCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAGACCCCGGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAT  >  AE015451/213651‑214099

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