Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,128,828 T→G 7.0% K93Q (AAG→CAG)  PP_4515 ← Transcriptional regulator, MarR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,128,8280TG 7.0% 194.1 / 7.3 57K93Q (AAG→CAG) PP_4515Transcriptional regulator, MarR family
Reads supporting (aligned to +/- strand):  ref base T (35/18);  new base G (2/2);  total (37/20)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

ACCTTGGGCACACCTTCTGAGGCTTGCACTGCACCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACAGGTTGAGCCCGGCCAGGCGGCGGTCCAGCTCGGCACGCCAACCTCGGGACATCTGG  >  AE015451/5128615‑5129060
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
aCCTTGGGCACACCTTCTGAGGCTTGCACTGCACCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTTAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCAcc                                                                                                                                                                                                                    <  1:637577/236‑1 (MQ=255)
    tGGGCACACCTTCTGAGGCTTGCACTGCACCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCt                                                                                                                                                                                                                >  1:519642/1‑236 (MQ=255)
        cacaCCTTCTGAGGCTTGCACTGCACCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCg                                                                                                                                                                                                            >  1:1172623/1‑236 (MQ=255)
        cacaCCTTCTGAGGCTTGCACTGCACCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGCACAATCTTCTTCGCCCGCCGATCCTCCAGCAccgcccgc                                                                                                                                                                                                              >  1:157227/1‑234 (MQ=255)
                   aGGCTTGCACTGCACCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGc                                                                                                                                                                                                  <  1:556398/235‑1 (MQ=255)
                        tgcactgcacCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTg                                                                                                                                                                                             >  1:800808/1‑235 (MQ=255)
                           actgcacCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGAGCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTcgcccgc                                                                                                                                                                                                                                   >  1:67461/1‑194 (MQ=255)
                                  cTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGagc                                                                                                                                                                                                                              >  1:249072/1‑190 (MQ=255)
                                        cGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGa                                                                                                                                                                            <  1:795841/236‑1 (MQ=255)
                                          aggggTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCg                                                                                                                                                                                                            <  1:599265/200‑1 (MQ=255)
                                            gggTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCCCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGACCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTg                                                                                                                                                                                                               <  1:1084382/197‑1 (MQ=255)
                                            gggTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATCCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTg                                                                                                                                                                                                               >  2:1084379/1‑197 (MQ=255)
                                             ggTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAg                                                                                                                                                                        >  1:887130/1‑235 (MQ=255)
                                                  gCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCTGCGGGCACAGTACAATCTTCTTCGCCCGACGATCCTCCAGCACCGCCTGCCGGCGCACCAAcc                                                                                                                                                                                                 >  2:857357/1‑205 (MQ=255)
                                                  gCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCg                                                                                                                                                                             >  1:944977/1‑225 (MQ=255)
                                                  gCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACgcgc                                                                                                                                                                  >  1:811765/1‑236 (MQ=255)
                                                                gAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGgcg                                                                                                                                                                                                         >  2:1178926/1‑183 (MQ=255)
                                                                gAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGccc                                                                                                                                                                                                >  1:18630/1‑192 (MQ=255)
                                                                gAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCgg                                                                                                                                                         >  1:909654/1‑231 (MQ=255)
                                                                                      tGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGc                                                                                                                                                                                                              <  2:471472/156‑1 (MQ=255)
                                                                                         aCGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTc                                                                                                                                                                              >  1:511591/1‑185 (MQ=255)
                                                                                            cgcATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTgcg                                                                                                                                                                                           <  1:457417/169‑1 (MQ=255)
                                                                                                  cAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTc                                                                                                                           <  1:1022283/227‑1 (MQ=255)
                                                                                                        tGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGCCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCgggg                                                                                                             >  1:806077/1‑235 (MQ=255)
                                                                                                          cAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTg                                                                                                                                                                                >  2:635214/1‑166 (MQ=255)
                                                                                                          cAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGAACCTCCAGccc                                                                                                                                                                                                                     >  1:289046/1‑127 (MQ=255)
                                                                                                            aGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGAccc                                                                                                                                              >  2:1132753/1‑196 (MQ=255)
                                                                                                            aGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCCCCCGCCGGCGCACCCGCCCtt                                                                                                                                                                                              >  2:113818/1‑150 (MQ=255)
                                                                                                                         cACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGcc                                                                                            <  1:516146/235‑1 (MQ=255)
                                                                                                                          aCGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAg                                                                                          >  1:545857/1‑236 (MQ=255)
                                                                                                                               ggTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGc                                                                                                                             >  2:787890/1‑196 (MQ=255)
                                                                                                                                gTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACg                                                                                                                                               <  2:402147/177‑1 (MQ=255)
                                                                                                                                    aCAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCGAAGAAGATTGTACTGTGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCgggg                                                                                                             >  1:1121170/1‑207 (MQ=255)
                                                                                                                                          ctACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGcttctt                                                                                                       >  2:1225676/1‑207 (MQ=255)
                                                                                                                                           tACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGctct                                                                                                                                     <  1:796469/176‑1 (MQ=255)
                                                                                                                                           tACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTg                                                                                                                                    >  1:396738/1‑177 (MQ=255)
                                                                                                                                           tACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACAccagc                                                                          >  1:588853/1‑235 (MQ=255)
                                                                                                                                            aCCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACAccagcc                                                                         <  1:635217/235‑1 (MQ=255)
                                                                                                                                             cccGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCCGGCTGTCGAGCAGACGCGCCAGGGGCGGGCCTTCGACGCCTACGCTCTGGGCCCGCTCCCGCTGGGTCggg                                                                                                              >  1:215643/1‑197 (MQ=255)
                                                                                                                                             cccGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGCgccgca                                                                                  >  1:1026014/1‑225 (MQ=255)
                                                                                                                                              ccGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACAccagccag                                                                       >  1:1201787/1‑235 (MQ=255)
                                                                                                                                                gCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCGAAGAAGATTGTACTGTGCCCGCCGATCCTCCAGCACCGcc                                                                                                                                                                                                                 <  1:1258356/95‑1 (MQ=21)
                                                                                                                                                gCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCGAAGAAGATTGTACTGTGCCCGCCGATCCTCCAGCACCGcc                                                                                                                                                                                                                 >  2:1258353/1‑95 (MQ=16)
                                                                                                                                                gCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAgcg                                                                     >  1:1243689/1‑235 (MQ=255)
                                                                                                                                                     tGGTTGGCGTTGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTCCGATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCtt                                                                                                                                                                                        <  2:902259/115‑1 (MQ=255)
                                                                                                                                                     tGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCGAAGAAGATTGTACTGTGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCa                                                                                                                          <  2:1121167/177‑1 (MQ=255)
                                                                                                                                                      ggTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAgcg                                                                     <  1:723712/229‑1 (MQ=255)
                                                                                                                                                          ggCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGGCCAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCCCCgg                                                                                                                                                                                                                  >  2:382153/1‑83 (MQ=255)
                                                                                                                                                          ggCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGAc                                                                                                                                                                      >  2:243044/1‑128 (MQ=255)
                                                                                                                                                                    tcGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACa                                                          <  1:704205/226‑1 (MQ=255)
                                                                                                                                                                    tcGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACAGGTTGAGcc                                                 >  1:779032/1‑235 (MQ=255)
                                                                                                                                                                     cGATCTGACCGATCAACGGCTTGGCCGGCGGGGACAGTACAATCTGCTTCGCcagccg                                                                                                                                                                                                                                 <  2:75724/58‑1 (MQ=255)
                                                                                                                                                                               gatcAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACAGGTTGAGCCCGGCCAggcggc                                     >  1:1067693/1‑236 (MQ=255)
                                                                                                                                                                                        gaggatcggcggGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCCGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCCCACGCCGTGTCgg                                                                                                               >  2:1010922/7‑153 (MQ=255)
                                                                                                                                                                                                gcggGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTggg                                                             >  2:1098376/1‑195 (MQ=255)
                                                                                                                                                                                                          tACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTg                                                                                                                                                                                                               <  2:1127172/39‑1 (MQ=255)
                                                                                                                                                                                                            cAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACAGGTTGAGCCCGGCCAGGCGGCGGTCCAGCTCGGCACGCCAACCTc             <  1:1202598/231‑1 (MQ=255)
                                                                                                                                                                                                                     tcGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTg                                                                                       <  1:844396/148‑1 (MQ=255)
                                                                                                                                                                                                                     tcGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACAGGtt                                                      >  2:89788/1‑181 (MQ=255)
                                                                                                                                                                                                                     tcGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACAGGTTGAGCCCGGCCAGGCGGCGGTCCAGCTCGGCACGCCAACCTCGGGACATCTgg  >  2:589984/1‑233 (MQ=255)
                                                                                                                                                                                                                     tcGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCCCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACAGGTTGAGCCCGGCCAGGCGGCGGTCCAGCTCGGCACGCCAACCTCGGGACATCTgg  <  1:895393/233‑1 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
ACCTTGGGCACACCTTCTGAGGCTTGCACTGCACCTGCATCGATGGGTGAGCTCAGGACTTTTCGAGGTTGGCCAGGATCTTGGCATGCACGCGCATACAGACCTGCAAGTCCGCCTCGTCCACGCCGGTAAACAGCTCTACCCGCAGCTGGTTGGCGATGGTCTCGATCTGATCGATCAACGGCTTGGCCGGCGGGCACAGTACAATCTTCTTCGCCCGCCGATCCTCCAGCACCGCCTGCCGGCGCACCAGCCCTTGCGCTTCCAGGCTGTCGAGCAGACGCGCCAGGGTCGGGCCTTCGACGCCTACGCTCTGGGCCAGCTCACGCTGTGTCGGGGCTTCTTCGAAGCGGGCCAGGTGCAGCAACACCAGCCAGCGCGCCTGGGACAGGTTGAGCCCGGCCAGGCGGCGGTCCAGCTCGGCACGCCAACCTCGGGACATCTGG  >  AE015451/5128615‑5129060

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