Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 200,338 G→C 7.0% V1948V (GTG→GTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451200,3380GC 7.0% 216.3 / 6.9 68V1948V (GTG→GTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (26/37);  new base C (3/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

GCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGT  >  AE015451/200106‑200547
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                 
gCAACAGTGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCGTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCCACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCAcc                                                                                                                                                                                                                <  1:1048288/236‑1 (MQ=14)
          cgcgccctcaacggttgcgatcgacatcacgacagggtcgttagcgccgacgtagacgttgtccggagcagtcacggtggcagtgcccacggtgccgttggccggcacagtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggttgttgatcggcagaccgtctttgttggtcagcgtgatggtgtaagtgatctcaccgCCTTCGGTCACCGAAGGCGGTg                                                                                                                                                                                                      >  1:1080262‑M1/215‑236 (MQ=255)
               tCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATc                                                                                                                                                                                                 <  1:168383/236‑1 (MQ=255)
                 aCCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                                                                                                                                               <  1:621934/236‑1 (MQ=255)
                  ccATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                                                                                                                                               <  1:578039/235‑1 (MQ=255)
                            tGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATc                                                                                                                                                                                     <  1:940820/235‑1 (MQ=255)
                                gAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCAccc                                                                                                                                                                                 >  1:407884/1‑235 (MQ=255)
                                               tACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGTCAGTAGACGCCACTGGTGATGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                                                                                                                                               <  1:1059220/206‑1 (MQ=255)
                                                 cacCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                               <  1:779197/236‑1 (MQ=37)
                                                    ccACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGt                                                                                                                                                                                                       <  2:441481/193‑1 (MQ=255)
                                                            cgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaa                                                                                                                                                     >  1:594235/1‑235 (MQ=40)
                                                             gcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCAccc                                                                                                                                                                                 <  1:438985/206‑1 (MQ=40)
                                                                  gCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACg                                                                                                                                                                  <  1:209249/216‑1 (MQ=39)
                                                                  gCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCAc                                                                                                                                              <  1:605921/236‑1 (MQ=255)
                                                                       gACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCg                                                                                                                                                          <  1:350084/219‑1 (MQ=39)
                                                                             tATCCGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCt                                                                                                                                                              <  1:61396/209‑1 (MQ=255)
                                                                               tCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACACTCGACGAAGTAGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCAACAACATCAc                                                                                                                                                                                   >  1:14346/1‑186 (MQ=255)
                                                                                 aGGACGCTGGCCAGATCAGCCTGGGCATCAGCTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTc                                                                                                                               <  1:1083032/236‑1 (MQ=255)
                                                                                    cggtggcagtgcccacggtgccgttggccggcacagtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggttgttgatcggcagaccgtctttgttggtcagcgtgatggtgtaggtgatctcgccgCCTTCGGTGACCGAAgg                                                                                                                                                                                                           <  2:593264‑M1/17‑1 (MQ=255)
                                                                                        tGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaacaa                                                                                                                                                  >  2:725069/1‑210 (MQ=255)
                                                                                          gCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaa                                                                                                                                                     <  1:896280/205‑1 (MQ=255)
                                                                                          gCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaa                                                                                                                                                     >  2:896277/1‑205 (MQ=40)
                                                                                           ccAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGAcc                                                                                                                      >  1:78004/1‑235 (MQ=40)
                                                                                             aGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACg                                                                                                                                                                  <  2:866430/189‑1 (MQ=255)
                                                                                             aGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGaa                                                                                                                                                        >  2:1083029/1‑198 (MQ=255)
                                                                                                tCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACgg                                                                                                                 <  1:982254/235‑1 (MQ=40)
                                                                                                       gggCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                >  2:801388/1‑197 (MQ=14)
                                                                                                       gggCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc                                                                                                                                                <  1:801391/197‑1 (MQ=14)
                                                                                                           aTCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgtcgt                                                                                                     >  1:367797/1‑236 (MQ=255)
                                                                                                                cTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTCACCGAAGGCGGTg                                                                                                                                                                                                      >  2:14537/1‑134 (MQ=255)
                                                                                                                 tCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCt                                                                                                                                      >  1:229785/1‑197 (MQ=40)
                                                                                                                     cggTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGGCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCAc                                                                                                                                                                                               <  1:1024124/134‑1 (MQ=255)
                                                                                                                      aGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGt                                                                                                                                    <  1:935576/194‑1 (MQ=39)
                                                                                                                      aGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGcc                                                                                           >  1:635535/1‑235 (MQ=40)
                                                                                                                                    gTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                <  1:688126/168‑1 (MQ=39)
                                                                                                                                     tgctgcATTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATc                                                                                                                                                                                     <  1:887554/130‑1 (MQ=17)
                                                                                                                                            ttCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTCACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACTCGGAGc                                                                                                                                       >  2:3669/1‑168 (MQ=255)
                                                                                                                                               gagaACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAg                                                                                                                            >  2:775181/1‑177 (MQ=35)
                                                                                                                                                 gaACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGtt                                                                           >  2:372765/1‑224 (MQ=39)
                                                                                                                                                        gAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCAcc                                                        <  1:906162/236‑1 (MQ=255)
                                                                                                                                                           cTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCa                                                       <  2:635532/234‑1 (MQ=14)
                                                                                                                                                           cTGGGCGGTGCTGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCAcc                                                        <  2:426070/233‑1 (MQ=25)
                                                                                                                                                            tGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCa                                                       <  2:594232/233‑1 (MQ=255)
                                                                                                                                                              gcactgtggttgttgatcggcagaccatctttgttggtcagggtgatggtgtaggtgatctcaccaCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGCTGGTCTGCCGATCAACAACCAc                                                                                                                                              <  1:294872‑M1/78‑1 (MQ=255)
                                                                                                                                                                     ggttattgatcggcaagccgtctttgttggtcagggtaatggtgtaagtgatctcaccgCCTTCGGTGACCGAAGGCGGTGAGATCACTTATACCATTACCCTGACCAACAAAGACg                                                                                                                                                                  <  2:912280‑M1/58‑1 (MQ=255)
                                                                                                                                                                      cGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTCCTTCAAGCTGACCGACGGCACCACtgtcgt                                                                                                     >  2:174132/1‑177 (MQ=2)
                                                                                                                                                                      cGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATg                                          >  1:927612/1‑236 (MQ=37)
                                                                                                                                                                        cTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGt                                                                                                                                    <  2:888745/144‑1 (MQ=32)
                                                                                                                                                                          tCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGt                                         <  2:407881/233‑1 (MQ=38)
                                                                                                                                                                                cAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGCCTAt                                           >  1:107323/1‑225 (MQ=37)
                                                                                                                                                                                      aCCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACgcgc                           <  1:658737/235‑1 (MQ=38)
                                                                                                                                                                                           cacCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGt                                                                            <  1:509190/181‑1 (MQ=37)
                                                                                                                                                                                           cacCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGt                                                                            >  2:509187/1‑181 (MQ=35)
                                                                                                                                                                                             ggttggtcagggtaatggtgtaagtgatctcaccgCCTTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACTACTGTCGTCGTGGCAGCCAACAGca                                                                                    >  1:451587‑M1/36‑171 (MQ=255)
                                                                                                                                                                                              ccTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCggt                       <  2:927609/231‑1 (MQ=38)
                                                                                                                                                                                                 cgacgaAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGCTCGGCCACTGCAACCGCACCAg                                                      >  2:1067893/1‑197 (MQ=14)
                                                                                                                                                                                                 cgacgaAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACGGGtt                                                                           >  2:798824/1‑176 (MQ=17)
                                                                                                                                                                                                      aaGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCAcc                                                        <  2:78002/190‑1 (MQ=39)
                                                                                                                                                                                                           gtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGc              >  2:820581/1‑227 (MQ=39)
                                                                                                                                                                                                           gtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCCATGTCGGCACCAACGCGCCg                         >  2:532125/1‑216 (MQ=39)
                                                                                                                                                                                                           gtggCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCAGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGCTCTGCCCATCAACAACCACTc                                                                                                                                            >  2:64252/1‑101 (MQ=255)
                                                                                                                                                                                                             ggCCAAGCTCACTCCCACCCCTTCGGTGACCGAAGGCGGTGAGCTCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCg         <  1:124266/230‑1 (MQ=255)
                                                                                                                                                                                                                  aggtgatctcaccgCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGCCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATc          <  1:454162‑M1/210‑1 (MQ=255)
                                                                                                                                                                                                                             aCCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGt  <  2:405995/221‑1 (MQ=25)
                                                                                                                                                                                                                                CCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGAt                                                                                                                                                        >  2:1042113‑M1/1‑68 (MQ=255)
                                                                                                                                                                                                                                CCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGAt                                                                                                                                                        <  1:1042116‑M1/68‑1 (MQ=255)
                                                                                                                                                                                                                                  TTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATg               >  1:728956‑M1/1‑203 (MQ=255)
                                                                                                                                                                                                                                    CGGTGACCGAAGGCGGTGAGATCACCTACAGCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGa                                                                                                                                         >  2:31616‑M1/1‑79 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                 
GCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGATGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTAGACGCCACTGGTGCTGCGTTCGAGAACCTTGAGCTGGGCGGTTCCGCTTCGGTTCAAGTCACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCTTCGGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGT  >  AE015451/200106‑200547

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