Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,706 C→T 6.9% N404N (AAC→AAT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,7060CT 6.9% 241.2 / 2.2 83N404N (AAC→AATPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (38/39);  new base T (4/2);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCG  >  AE015451/195525‑195940
                                                                                                                                                                                     |                                                                                                                                                                                                                                          
gTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                      <  1:194325/236‑1 (MQ=255)
   cGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGCTgt                                                                                                                                                                                    >  1:913753/1‑234 (MQ=255)
   cGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACAGTGACCCTGTCCCACGGCCa                                                                                                                                                                                                   >  1:666943/1‑220 (MQ=9)
                           aCGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                                                                                                                                                                                      <  1:488309/161‑1 (MQ=40)
                                 ccGGCGGTAACTTCCAGCCACTGACATCGAATCCGACCCCGGCTCAGACCACGATCAACGACACAGTCGACACCACCACCGCGACCCTGACGGCGACCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAcc                                                                                                                                                                                                                       <  1:906995/170‑1 (MQ=14)
                                  aggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCa                                                                                                                                                                                                                                    >  2:392176/2‑156 (MQ=25)
                                  aggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCa                                                                                                                                                                                                                                    <  1:392179/155‑1 (MQ=25)
                                          aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGacc                                                                                                                                                                                                                 >  2:620516/1‑165 (MQ=2)
                                            ttCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                               >  2:517370/1‑197 (MQ=2)
                                            ttCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                               <  1:517373/199‑3 (MQ=9)
                                              cGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGGTCAACGACACGGTCGACACCACCACCGCGACCCTGACCGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                        <  1:954675/172‑1 (MQ=255)
                                               gAGCAACTGACCCCGACCCCGACCCCGGCTCAGACCACGCTCAACGACGCGGTCGACACCACCACCGCGCCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTaa                                                                                                                                                                                                                         <  1:230839/154‑1 (MQ=255)
                                                 gCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                                                                                                                                                                                      <  1:592506/139‑1 (MQ=40)
                                                              aaTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGt                                                                                                                                                                                                             <  1:842051/151‑1 (MQ=255)
                                                               aTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                        <  2:334152/155‑1 (MQ=255)
                                                                  cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                              <  2:694483/130‑1 (MQ=18)
                                                                  cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                              >  1:694486/1‑130 (MQ=18)
                                                                  cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGc                                                                                                                                                                                                                               <  1:217701/129‑1 (MQ=25)
                                                                             cAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGTCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGt                                                                                                                                                                                                             <  2:709263/136‑1 (MQ=255)
                                                                                 ccACGATGAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                                <  1:167749/145‑1 (MQ=255)
                                                                                       tCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                                    >  1:193678/1‑119 (MQ=39)
                                                                                         accGACACGGTCGACACCACAACCGCGACCCTGACCGCGAGGCCAACCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                        <  1:543927/127‑1 (MQ=11)
                                                                                             aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGt                                                                                                                                                                                                             >  1:984029/1‑120 (MQ=12)
                                                                                                  gtcgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAt                                                                                                                                                        >  1:1140349/8‑168 (MQ=255)
                                                                                                    cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                                >  2:219596/6‑126 (MQ=2)
                                                                                                    cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                                      >  2:536278/6‑120 (MQ=255)
                                                                                                    cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                                      <  1:536281/115‑1 (MQ=255)
                                                                                                    cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                                <  1:219596/121‑1 (MQ=25)
                                                                                                    cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTg                                                                                                                                                                                                              <  1:694146/112‑1 (MQ=255)
                                                                                                    cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTg                                                                                                                                                                                                              >  2:694143/1‑112 (MQ=255)
                                                                                                    cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                      <  2:531115/184‑1 (MQ=255)
                                                                                                    cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                      >  1:531118/1‑184 (MQ=38)
                                                                                                               ccGCGCCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                                    <  2:70329/95‑1 (MQ=255)
                                                                                                               ccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                              >  2:433089/1‑85 (MQ=35)
                                                                                                               ccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                              <  1:433092/85‑1 (MQ=35)
                                                                                                                cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                    >  2:779708/1‑126 (MQ=40)
                                                                                                                cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                    <  1:779711/126‑1 (MQ=40)
                                                                                                                cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAACCCAGGGCAGCGTCGATTTCCAGa                                                                                                                                                 >  1:1161812/1‑159 (MQ=255)
                                                                                                                   gACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGccc                                                                                                                                                                                                                                     <  2:1268198/74‑1 (MQ=34)
                                                                                                                       cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                              <  1:1188983/77‑1 (MQ=34)
                                                                                                                       cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                              >  1:31383/1‑77 (MQ=34)
                                                                                                                       cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                              >  2:1188980/1‑77 (MQ=34)
                                                                                                                       cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCga                                                                                             >  1:1103778/1‑206 (MQ=9)
                                                                                                                          aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                                    >  1:658166/1‑84 (MQ=255)
                                                                                                                          aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                                    <  2:658163/84‑1 (MQ=255)
                                                                                                                                    caaccgTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                               <  2:804354/106‑1 (MQ=255)
                                                                                                                                        ggTGACCGAAAGCGGCGAGAACACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                              <  1:875279/60‑1 (MQ=11)
                                                                                                                                          tCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                  <  1:462178/134‑1 (MQ=21)
                                                                                                                                           cACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                              >  2:226787/1‑57 (MQ=35)
                                                                                                                                              cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTac                                                                                                                             >  1:697662/1‑151 (MQ=38)
                                                                                                                                              cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTac                                                                                                                             <  2:697659/151‑1 (MQ=39)
                                                                                                                                                      gcgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCa                                                                                                                                                                                                                                    <  2:817317/40‑1 (MQ=40)
                                                                                                                                                       cgTGATCTCTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                              <  2:447465/45‑1 (MQ=255)
                                                                                                                                                       cgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGccc                                                                                                                  >  1:226700/1‑151 (MQ=39)
                                                                                                                                                          gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACtt                                                                        >  2:627247/1‑192 (MQ=255)
                                                                                                                                                          gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACtt                                                                        <  1:627250/192‑1 (MQ=40)
                                                                                                                                                              acctACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                  <  2:1145788/111‑1 (MQ=0)
                                                                                                                                                              aCTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTAACCCTGTCCAACGGCCATACGACTACCGTTGAAGCGGGAAAAACCCAgg                                                                                                                                                                   >  1:606335/1‑97 (MQ=17)
                                                                                                                                                                  acacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcc                                                                                                                                         <  2:359822/119‑1 (MQ=17)
                                                                                                                                                                    acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgaga                                                                                           <  2:727776/163‑1 (MQ=255)
                                                                                                                                                                    acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgaga                                                                                           >  1:727779/1‑163 (MQ=255)
                                                                                                                                                                    acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACa                                                                  <  1:235596/188‑1 (MQ=39)
                                                                                                                                                                     ccGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATCACTGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgaga                                                                                           <  2:459335/162‑1 (MQ=255)
                                                                                                                                                                      cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCg                                                                                                                 <  2:946528/139‑1 (MQ=255)
                                                                                                                                                                         gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                                    >  2:808214/1‑37 (MQ=40)
                                                                                                                                                                         gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAggg                                                                                                                                                                  >  2:932941/1‑87 (MQ=11)
                                                                                                                                                                           cccTGAGCAACCCTGCCCAGACGCCGGTAAGCCTGAACCACTCCAAc                                                                                                                                                                                                        >  1:464179/1‑47 (MQ=14)
                                                                                                                                                                           cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                               >  1:476056/1‑72 (MQ=21)
                                                                                                                                                                           cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCAcc                                                                                   >  2:649518/1‑163 (MQ=38)
                                                                                                                                                                           cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgag                                                                                            >  1:236171/1‑155 (MQ=38)
                                                                                                                                                                           cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACaa                                                                 >  1:389921/1‑181 (MQ=39)
                                                                                                                                                                           cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGGCACCac                                                                                                >  2:142776/1‑150 (MQ=255)
                                                                                                                                                                             cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                               >  2:879409/1‑70 (MQ=255)
                                                                                                                                                                             cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                               <  1:879412/70‑1 (MQ=25)
                                                                                                                                                                             cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccac         >  1:788633/1‑236 (MQ=255)
                                                                                                                                                                             cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCg                    >  2:1091521/1‑225 (MQ=255)
                                                                                                                                                                             atGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                               <  2:25499/69‑1 (MQ=255)
                                                                                                                                                                                 gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                        >  1:812475/1‑41 (MQ=255)
                                                                                                                                                                                 gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCa                                                                            >  1:1221827/1‑165 (MQ=40)
                                                                                                                                                                                 gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCCACAACAACGGCTCGCCCGTCAGCGGCACCATCGAGAgcgc                                                                                       >  2:235596/1‑154 (MQ=255)
                                                                                                                                                                                   aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                               <  1:378327/64‑1 (MQ=255)
                                                                                                                                                                                   aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                               >  2:378324/1‑64 (MQ=255)
                                                                                                                                                                                   aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAAccc                                                                                                                                                                      >  1:474211/1‑73 (MQ=255)
                                                                                                                                                                                     cccTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACcgcg  <  1:843525/235‑1 (MQ=255)
                                                                                                                                                                                     |                                                                                                                                                                                                                                          
GTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCG  >  AE015451/195525‑195940

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