Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 388,722 T→A 8.0% L198Q (CTG→CAG)  soxB → sarcosine oxidase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451388,7220TA 8.0% 168.5 / 6.1 50L198Q (CTG→CAG) soxBsarcosine oxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (27/19);  new base A (2/2);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

CCTCTGCCACACCCTGCAGGACATGCGTGACTCCGAGCGTCGGGTCAGCGCCAACCGCCTCAACGGCGTGGATGGCGAGCTGCTGAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGCCGCTGGAATCGCACCCGCTGCAAGCGCTGGTATCCGAGCCGATCAAGCCG  >  AE015451/388489‑388936
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                      
ccTCTGCCACACCCTGCAGGACATGCGTGACTCCGAGCGTCGGGTCAGCGCCAACCGCCTCAACGGCGTGGATGGCGAGCTGCTGAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTgg                                                                                                                                                                                                                      <  1:357048/236‑1 (MQ=255)
 ctctGCCACACCCTGCAGGACATGCGTGACTCCGAGCGTCGGGTCAGCGCCAACCGCCTCAACGGCGTGGATGGCGAGCTGCTGAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTggg                                                                                                                                                                                                                     <  1:82121/236‑1 (MQ=255)
                        gCGTGACTCCGAGCGTCGGGTCAGCGCCATCCGCCTCAACGGCGTGGATGGCGAGCTGCTGAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATc                                                                                                                                                                                                        >  2:487896/1‑226 (MQ=255)
                                                                       aTGGCGAGCTGCTGAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGGCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAgg                                                                                                                                                                         >  1:814981/1‑210 (MQ=255)
                                                                               ctgctgAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAg                                                                                                                                                                                           <  2:418359/184‑1 (MQ=255)
                                                                                   tgAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGt                                                                                                                                                            <  1:281360/211‑1 (MQ=255)
                                                                                      acacTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCtt                                                                                                                                 <  1:10478/235‑1 (MQ=255)
                                                                                             ccAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTcc                                                                                                                                                                               >  2:677086/1‑182 (MQ=255)
                                                                                               aGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGAcc                                                                                                                                                                                                             >  2:595181/1‑150 (MQ=255)
                                                                                               aGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGgc                                                                                                                          >  2:401516/1‑233 (MQ=255)
                                                                                                       gCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTc                                                                                                                                                           <  2:397754/192‑1 (MQ=255)
                                                                                                             aTCCCGTACCTGGACTGCTCGAAGAACGCCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTc                                                                                                                                                                                <  2:313710/165‑1 (MQ=255)
                                                                                                                  gTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTc                                                                                                                                                                                                                 <  1:683143/127‑1 (MQ=255)
                                                                                                                           cTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCa                                                                                                                                                          <  2:178906/173‑1 (MQ=255)
                                                                                                                                cGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTg                                                                                                        >  2:527035/1‑218 (MQ=255)
                                                                                                                                 gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGa                                                                                                                                                                        >  2:711474/1‑153 (MQ=255)
                                                                                                                                 gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCa                                                                                                                                                          >  2:568258/1‑167 (MQ=255)
                                                                                                                                 gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACt                                                                                        >  2:567070/1‑233 (MQ=255)
                                                                                                                                        gcccGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGt                                                                                                                                                            <  2:692547/157‑1 (MQ=255)
                                                                                                                                           cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAGt                                                                                                                                                                                          <  1:872829/125‑1 (MQ=255)
                                                                                                                                           cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAGt                                                                                                                                                                                          >  2:872826/1‑125 (MQ=255)
                                                                                                                                                cccGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCgg                                                                                    <  1:684078/222‑1 (MQ=255)
                                                                                                                                                cccGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCgg                                                                                    >  2:684075/1‑222 (MQ=255)
                                                                                                                                                       cTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTc                                                                                                              <  2:614653/189‑1 (MQ=255)
                                                                                                                                                           gcgcCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccgccg                                                                                              <  2:644233/201‑1 (MQ=255)
                                                                                                                                                              ccACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCg                                                         <  1:903506/235‑1 (MQ=255)
                                                                                                                                                                       aGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGCCGCt                                                <  1:401519/235‑1 (MQ=255)
                                                                                                                                                                            cGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGTGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgtgt                                                                                                                                                   >  2:205684/1‑131 (MQ=255)
                                                                                                                                                                            cGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGc                                                    >  2:453319/1‑226 (MQ=255)
                                                                                                                                                                              tggtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCggg                                                                                   >  2:771897/1‑193 (MQ=255)
                                                                                                                                                                              tggtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGcca                                                                                >  2:28189/1‑194 (MQ=255)
                                                                                                                                                                              tggtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTg                                                                    >  2:828137/1‑208 (MQ=255)
                                                                                                                                                                              tggtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGGGGGCAccgc                                                                                                >  2:229171/1‑180 (MQ=255)
                                                                                                                                                                              tggtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGCCATCGCTGTGGCAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGt                                                                                                         >  2:356425/1‑171 (MQ=255)
                                                                                                                                                                                 tggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGCCGCTg                                               >  2:808002/1‑226 (MQ=255)
                                                                                                                                                                                    cGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGCCGCTGGAATCGCACCCGc                                  <  1:836199/236‑1 (MQ=255)
                                                                                                                                                                                           cGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTaa                                                                                          >  2:905357/1‑173 (MQ=255)
                                                                                                                                                                                           cGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAg                                                                                                                                                                                           <  1:778015/76‑1 (MQ=255)
                                                                                                                                                                                           cGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAg                                                                                                                                                                                           >  2:778012/1‑76 (MQ=255)
                                                                                                                                                                                               aagATGCCGTGGCCTGGGTCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGAGCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAGCGGCGCGGTCATCGGTGTgg                                                                                                                                                 <  2:74703/112‑1 (MQ=255)
                                                                                                                                                                                                cGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGACGTGCTCGGCTTCCGCCAGGAAAACGGCGCGGTCATCGGCGTGGACACcaacc                                                                                                                                        >  2:39321/1‑121 (MQ=255)
                                                                                                                                                                                                          gCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGCCGCTGGAATCGCACCCGCTGCAAGCGCTGGTATCCGagc             <  1:482708/235‑1 (MQ=255)
                                                                                                                                                                                                           ccTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTTGCCGGCTTCCGCCTGCCGCTGGAATCGCACCCg                                   >  2:831030/1‑212 (MQ=255)
                                                                                                                                                                                                           ccTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGCCGCTgg                                              >  2:784026/1‑201 (MQ=255)
                                                                                                                                                                                                                  cTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGCTCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAACCCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGt                                                                                            >  2:12399/1‑148 (MQ=255)
                                                                                                                                                                                                                       cccGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGCCGCTGGAATCGCACCCGCTGCAAGCGCTGGTATCCGagccgatcaagccg  >  2:561927/1‑233 (MQ=255)
                                                                                                                                                                                                                       cccGCGCTGCCGACGCCCTGGCCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCCGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAg                                                                       >  2:840961/1‑164 (MQ=255)
                                                                                                                                                                                                                              ggtcGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCCCATGGCCAAGCTGGCCgg                                                               >  2:549883/4‑165 (MQ=255)
                                                                                                                                                                                                                                       ccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCt                                                                     <  1:493441/150‑1 (MQ=255)
                                                                                                                                                                                                                                       ccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCt                                                                     >  2:493438/1‑150 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                      
CCTCTGCCACACCCTGCAGGACATGCGTGACTCCGAGCGTCGGGTCAGCGCCAACCGCCTCAACGGCGTGGATGGCGAGCTGCTGAACACTGCCCAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGCTGGCCGGCTTCCGCCTGCCGCTGGAATCGCACCCGCTGCAAGCGCTGGTATCCGAGCCGATCAAGCCG  >  AE015451/388489‑388936

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