Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,128,535 G→T 8.3% intergenic (+33/+17) 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,5350GT 8.3% 162.3 / 7.0 48intergenic (+33/+17)PP_1887/PP_1888conserved protein of unknown function/Fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base T (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-02

AATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCATCTGAATACCTATGCCACGGGCTGCGGAGTCGTCGGTCA  >  AE015451/2128315‑2128763
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
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                                                                                                                                                                                  <  1:412020/200‑1 (MQ=255)
                                                                         ggTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTCCCTATACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCtt                                                                                                                                                                                  >  2:412019/1‑200 (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                                                                                                                                                                                            <  1:280597/131‑1 (MQ=255)
                                                                                                                                    gCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCa                                                                                                                                                                                            >  2:280596/1‑131 (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                                                                                                                                                       >  2:252786/1‑126 (MQ=255)
                                                                                                                                                                              cGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGtt                                                                                                                                                       <  1:252787/126‑1 (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)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
AATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGATGAGCAGGGGCGGGTACTGACGTATGACGGGCTCAACCGTTTGCTCAGGGTCGAGACCCCTGACGGTGAACGCTGTCGCTACAACTATGACCCAGAGAACATTCTGTCCGGCACGGCCCCTGGCTAAGGGCACCGGGGGGGCGACGAGGCCGGTATAGGGACCGGCCTCGTCGCATTCAGCGGTAGCGAATGGTGAAGTTCAAGGCACCCTTGGCCAAGCCCGTGCGCAGGCTGCTGTTGGGTTCTGTCTGGAAGAAGCGGGCCTGAAAATTCAAAGGCACGTTATCGCTCGTCAGGGGTACCAGGTCGATCTCTTTGTTCAGGGGCATCGATACCCCGTTGTACAGCATCTGAATACCTATGCCACGGGCTGCGGAGTCGTCGGTCA  >  AE015451/2128315‑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: