Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,424 G→A 10.4% V310V (GTG→GTA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,4240GA10.4% 193.1 / 7.5 76V310V (GTG→GTAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base A (3/5);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGC  >  AE015451/195195‑195655
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
gtCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccgtg                                                                                                                                                                                                                                   <  1:1408181/235‑1 (MQ=255)
             cacacCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                       <  1:739508/235‑1 (MQ=255)
               cacCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGc                                                                                                                                                                                                                    <  1:817010/236‑1 (MQ=255)
                    tACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGt                                                                                                                                                                                                               <  1:1346606/236‑1 (MQ=255)
                         cGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTCCCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtga                                                                                                                                                                                                                                        <  1:464980/206‑1 (MQ=255)
                              gTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGtt                                                                                                                                                                                                     >  1:1063846/1‑236 (MQ=255)
                               tGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCa                                                                                                                                                                                                                  >  1:907749/1‑222 (MQ=255)
                               tGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTg                                                                                                                                                                                                    >  1:978667/1‑236 (MQ=255)
                                     gATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGt                                                                                                                                                                                                               <  2:853320/219‑1 (MQ=255)
                                     gATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                                              <  1:993894/236‑1 (MQ=255)
                                                                cGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCaaa                                                                                                                                                                                                                >  1:1047213/1‑189 (MQ=255)
                                                                     gtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgacc                                                                                                                                                                                                                                      >  1:897159/1‑164 (MQ=255)
                                                                     gtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                                              >  1:1176642/1‑236 (MQ=255)
                                                                     gtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                                              >  1:1070043/1‑236 (MQ=255)
                                                                     gtggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                                              >  1:885565/1‑236 (MQ=255)
                                                                       ggCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGgtgaccgtgac                                                                                                                                                                                                                                 >  2:64729/1‑167 (MQ=255)
                                                                          gggCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACTATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                                   >  1:814219/1‑178 (MQ=255)
                                                                          gggCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCGAGACGCCGGTGACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                                   >  1:797564/1‑178 (MQ=255)
                                                                          gggCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                                                                         >  1:1493648/1‑236 (MQ=255)
                                                                           ggCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCCTGActcttat                                                                                                                                                                                                                           >  2:555797/1‑166 (MQ=255)
                                                                                   gTGACATTGCCGCAAACGCCAGCGGTTCCCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                                                                                                                                                              <  1:231223/190‑1 (MQ=255)
                                                                                                  aCGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACggc                                                                                                                                   >  1:948808/1‑233 (MQ=255)
                                                                                                         cGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACg                                                                                                                                                                                                                      <  2:1202149/144‑1 (MQ=255)
                                                                                                         cGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                                 <  2:1243966/197‑1 (MQ=255)
                                                                                                          ggTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACg                                                                                                                                     <  2:1105393/224‑1 (MQ=255)
                                                                                                                 cAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGCTCATCACCGTTGAAGCCGGCAAAACCCAGGGc                                                                                                                                                                                >  1:99130/1‑174 (MQ=255)
                                                                                                                   attcaacGATACGGTCGACACCACAACCGCGACCCGGACGGCCTCGCCAAGGGTGACCGAAGGTGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                                                                                                                       <  2:114911/127‑1 (MQ=11)
                                                                                                                                aTCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGt                                                                                                                                                                                                               <  2:485772/128‑1 (MQ=255)
                                                                                                                                  cGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                       <  2:1070040/182‑1 (MQ=255)
                                                                                                                                  cGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                       <  2:808161/182‑1 (MQ=255)
                                                                                                                                    acacTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                                   >  1:568108/1‑120 (MQ=255)
                                                                                                                                     cacTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACcgg                                                                                                <  2:1176639/234‑1 (MQ=255)
                                                                                                                                         aCAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGACCAATGACGTCTAcaacaa                                                                                                                                       >  1:565761/1‑191 (MQ=255)
                                                                                                                                            aCCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                       <  2:410444/108‑1 (MQ=255)
                                                                                                                                              cccgACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                       <  2:218402/104‑1 (MQ=255)
                                                                                                                                                  aCCCTGCCCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                                                  <  2:348557/171‑1 (MQ=255)
                                                                                                                                                  aCCCTGACCGCCTCGCCAAGCGTGACCGAAGGTAGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACCGGCGGCAACTTCGaac                                                                                  >  1:1149853/1‑233 (MQ=255)
                                                                                                                                                   cccTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGAttt                                                                                                                                                                     <  2:835827/151‑1 (MQ=255)
                                                                                                                                                    ccTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAAc                                                                                                                                                                                                                       <  1:482360/100‑1 (MQ=255)
                                                                                                                                                     cTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCaa                                                                                                                                                      <  2:26225/164‑1 (MQ=255)
                                                                                                                                                     cTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCaa                                                                                                                                                      >  1:26225/1‑164 (MQ=255)
                                                                                                                                                      tGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTg                                                                             <  1:355322/236‑1 (MQ=255)
                                                                                                                                                        ccctccTCGCCAAGCGTCACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACgg                                                                                                                                    <  2:424572/175‑1 (MQ=255)
                                                                                                                                                           gCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGt                                                                                                                           >  1:1017197/1‑185 (MQ=255)
                                                                                                                                                                gCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTTACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCa                                                                                                                                                                               >  1:1511213/1‑128 (MQ=255)
                                                                                                                                                                gCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATc                                                                   >  1:1037408/1‑236 (MQ=255)
                                                                                                                                                                     gcgcaaCCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCACCCCTGACCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACGCAGGGCAGCGTCGATTTCCAGa                                                                                                                                                                <  2:187901/133‑1 (MQ=14)
                                                                                                                                                                           ccGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCg                                                                       <  2:885563/221‑1 (MQ=255)
                                                                                                                                                                             gAAGGTGGCGTGATCACCTACACCGTGACCCTGCGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGt                                                                                                                                   <  2:474436/159‑1 (MQ=18)
                                                                                                                                                                                     cGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCaaac                                                                                                                                                                                                               >  1:272923/1‑72 (MQ=0)
                                                                                                                                                                                        gATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                                                       >  1:402870/1‑64 (MQ=14)
                                                                                                                                                                                        gATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                                                       <  2:402868/64‑1 (MQ=14)
                                                                                                                                                                                         aTCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAAc                                                                               >  1:329887/1‑199 (MQ=39)
                                                                                                                                                                                              cTACCCCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTg                                                                             <  2:194770/196‑1 (MQ=255)
                                                                                                                                                                                                   ccGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                       >  1:869137/1‑117 (MQ=35)
                                                                                                                                                                                                   ccGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAGCGATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                       <  2:869135/117‑1 (MQ=255)
                                                                                                                                                                                                       gACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGCCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACcggcggc                                                                                            >  2:652216/1‑171 (MQ=17)
                                                                                                                                                                                                         cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTc                                                                                                                                                                                                                           >  1:235909/1‑43 (MQ=255)
                                                                                                                                                                                                         cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCa                                                                                                                                                                  >  1:226562/1‑100 (MQ=37)
                                                                                                                                                                                                         cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGCCCAc                                                >  1:1009989/1‑214 (MQ=255)
                                                                                                                                                                                                         cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCAAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACt                                     <  1:905149/225‑1 (MQ=255)
                                                                                                                                                                                                         cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCa                                                                                                                                                                  <  2:226561/100‑1 (MQ=37)
                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTGACTGTGACCCTGt                                                                                                                                                                                                                            >  1:237779/1‑40 (MQ=255)
                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccac                        >  1:1173793/1‑236 (MQ=255)
                                                                                                                                                                                                           cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccac                        >  1:1439904/1‑236 (MQ=255)
                                                                                                                                                                                                               gCAACCCTGCCCAGACGCCGGTGACCGTGAccc                                                                                                                                                                                                                               >  2:117960/1‑33 (MQ=255)
                                                                                                                                                                                                                ccaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCATGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGcc                                                                                                     <  2:1184573/152‑1 (MQ=18)
                                                                                                                                                                                                                cAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                           <  2:897157/196‑1 (MQ=255)
                                                                                                                                                                                                                   tccTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                                       <  2:1183219/100‑1 (MQ=14)
                                                                                                                                                                                                                   cccTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                 <  2:1493645/187‑1 (MQ=255)
                                                                                                                                                                                                                    ccTGCCCAGACGCCGGTGACCGGTACCCGGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaa                                                                                                                                          <  2:65796/113‑1 (MQ=255)
                                                                                                                                                                                                                      tGCCCAGACGCCGGTGACCGTGACCCTGTCCAACgg                                                                                                                                                                                                                     >  1:338923/1‑36 (MQ=40)
                                                                                                                                                                                                                       gCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAACTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACg                                               >  2:712372/1‑201 (MQ=25)
                                                                                                                                                                                                                           aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                                                                                                                  >  1:1444161/1‑98 (MQ=34)
                                                                                                                                                                                                                           aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACCGTTGAAGCGGGTAAAACCCCGGGctcc                                                                                                                                                                             >  2:683347/1‑68 (MQ=11)
                                                                                                                                                                                                                            gACGCCGGTAACTGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGCCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCa                                                                                 <  2:823440/162‑1 (MQ=255)
                                                                                                                                                                                                                             ccGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                           <  2:13440/182‑1 (MQ=255)
                                                                                                                                                                                                                               gCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                 <  1:1187655/175‑1 (MQ=255)
                                                                                                                                                                                                                                ccGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAAccc                                                                  >  2:916407/1‑173 (MQ=38)
                                                                                                                                                                                                                                ccGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGCCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTAcaacaa                                                                                                                                       >  1:152721/1‑104 (MQ=17)
                                                                                                                                                                                                                                    tgaccgtgaccCTGTCCAACGGCCAGGTCATCACTGtt                                                                                                                                                                                                     >  1:76757/1‑38 (MQ=37)
                                                                                                                                                                                                                                     gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGGCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGc  >  2:996753/1‑232 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
TTCCGGCAAACGACACACCTTACATCGACGGTGGCCAGATTTCGGCCACCGTGACAGGTAGTACCGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGC  >  AE015451/195195‑195655

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: