Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,128,532 A→C 8.5% intergenic (+30/+20) PP_1887 → / ← PP_1888 conserved protein of unknown function/Fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,128,5320AC 8.5% 154.4 / 6.4 47intergenic (+30/+20)PP_1887/PP_1888conserved protein of unknown function/Fimbrial protein
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base C (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

TTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCATCTGAATACCTATGCCACGGGCTGCGGAGTCGTCGGTCA  >  AE015451/2128302‑2128763
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       
ttCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAgg                                                                                                                                                                                                                                       >  1:352230/1‑233 (MQ=255)
             aaTACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCtc                                                                                                                                                                                                                             >  1:13130/1‑230 (MQ=255)
                       gCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTa                                                                                                                                                                                                             <  1:457395/236‑1 (MQ=255)
                                        tGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCa                                                                                                                                                                                            <  1:328417/236‑1 (MQ=255)
                                             tACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCa                                                                                                                                                                                       <  1:610835/236‑1 (MQ=255)
                                                                         atgaGCAGGGGCGGGTACTGACGGATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCa                                                                                                                                                                                            <  2:13130/203‑1 (MQ=255)
                                                                               aGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGTACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCtcgtcg                                                                                                                                                                                                                         >  1:351476/1‑168 (MQ=255)
                                                                                      ggTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTCCCTATACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCtt                                                                                                                                                                                  >  2:412019/1‑200 (MQ=255)
                                                                                      ggTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTCCCTATACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCtt                                                                                                                                                                                  <  1:412020/200‑1 (MQ=255)
                                                                                         aCTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTCCCTATACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTgc                                                                                                                                                                    >  2:457614/1‑211 (MQ=255)
                                                                                         aCTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTCCCTATACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTgc                                                                                                                                                                    <  1:457615/211‑1 (MQ=255)
                                                                                         aCTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGt                                                                                                                                                                      >  2:742474/1‑209 (MQ=255)
                                                                                             aCGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCa                                                                                                                                                                                       <  2:352229/188‑1 (MQ=255)
                                                                                                   gACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCtcgtc                                                                                                                                                                                                                          >  1:705032/1‑147 (MQ=255)
                                                                                                              ccGTTTGCGCAGGGTCGAGACCCCTGTCGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGcc                                                                                                                                                                              <  2:210754/180‑1 (MQ=255)
                                                                                                                                 aCCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTgg                                                                                                                                                     <  2:46515/186‑1 (MQ=255)
                                                                                                                                  ccccTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTc                                                                                                   >  1:290875/1‑235 (MQ=255)
                                                                                                                                         cGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCa                                                                                                                                                                                       <  1:645843/144‑1 (MQ=255)
                                                                                                                                                 gCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCa                                                                                                                                                                                            >  2:280596/1‑131 (MQ=255)
                                                                                                                                                 gCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCa                                                                                                                                                                                            <  1:280597/131‑1 (MQ=255)
                                                                                                                                                          aCAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAgg                                                                                                                                                               <  2:466347/151‑1 (MQ=255)
                                                                                                                                                             aCTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTg                                                                        <  1:742475/235‑1 (MQ=255)
                                                                                                                                                                aTGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCa                                                                                                                                                                                                                   >  2:383686/1‑93 (MQ=255)
                                                                                                                                                                    cccAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTc                                                                                                   >  1:207125/1‑201 (MQ=255)
                                                                                                                                                                      cAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCGTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGcc                                                                                                                                                                         >  2:455520/1‑129 (MQ=255)
                                                                                                                                                                      cAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCGGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCa                                                              >  1:220935/1‑236 (MQ=255)
                                                                                                                                                                       agagAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTa                                                                                                         <  2:310523/192‑1 (MQ=255)
                                                                                                                                                                          gaACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGt                                                                                                                                                                      <  2:145654/128‑1 (MQ=255)
                                                                                                                                                                                     ccGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGTCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTgg                                                                                                                                                     <  2:321857/134‑1 (MQ=255)
                                                                                                                                                                                           cGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGtt                                                                                                                                                       <  1:252787/126‑1 (MQ=255)
                                                                                                                                                                                           cGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGtt                                                                                                                                                       >  2:252786/1‑126 (MQ=255)
                                                                                                                                                                                                 cTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTg                                                <  1:641198/223‑1 (MQ=255)
                                                                                                                                                                                                 cTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTg                                                >  2:641197/1‑223 (MQ=255)
                                                                                                                                                                                                 cTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCATCTGaa                                   >  1:654195/1‑236 (MQ=255)
                                                                                                                                                                                                    gtaaaGGGAACGGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAg                                           <  1:273413/222‑1 (MQ=255)
                                                                                                                                                                                                           gCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCATCTGAATACCTTTGc                          >  1:599901/1‑235 (MQ=255)
                                                                                                                                                                                                             aCCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGt                                                                                                 >  2:314781/1‑162 (MQ=255)
                                                                                                                                                                                                             aCCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGt                                                                                                 <  1:314782/162‑1 (MQ=255)
                                                                                                                                                                                                                 ggggggCGACGAGGCCGGTATAGGGACCGGCCtcgt                                                                                                                                                                                                                           >  2:486937/1‑36 (MQ=255)
                                                                                                                                                                                                                  gggggCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGt                                                                                                           >  1:45543/1‑147 (MQ=255)
                                                                                                                                                                                                                  gggggCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGt                                                                                                           <  2:45543/147‑1 (MQ=255)
                                                                                                                                                                                                                       cgacgaGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCATCTGAATACCTATGCCACGGGCt                  >  1:747625/1‑231 (MQ=255)
                                                                                                                                                                                                                        gacgaGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGAtc                                                                              <  1:513805/170‑1 (MQ=255)
                                                                                                                                                                                                                        gacgaGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGAtc                                                                              >  2:513804/1‑170 (MQ=255)
                                                                                                                                                                                                                         acgaGGCCGGTATAGGGACCGGCCtcg                                                                                                                                                                                                                            <  2:210556/27‑1 (MQ=255)
                                                                                                                                                                                                                            aGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGc                                                               >  2:633911/1‑181 (MQ=255)
                                                                                                                                                                                                                                  ggataGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCAt                                        <  2:536038/196‑1 (MQ=255)
                                                                                                                                                                                                                                  gTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCATCTGAATACCTATGCCACGGGCTGCGGAGTCGTCGGTCa  >  1:477435/1‑236 (MQ=255)
                                                                                                                                                                                                                                      aGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGTCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGGTGTACAGCATCTGAATACCTATGCCACGGGc                   >  1:235385/1‑215 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       
TTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCATCTGAATACCTATGCCACGGGCTGCGGAGTCGTCGGTCA  >  AE015451/2128302‑2128763

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