Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,022,307 G→C 8.8% intergenic (‑2/+666) PP_3547 ← / ← emrB putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/multidrug efflux transport system ‑ permease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,022,3070GC 8.8% 225.1 / 13.5 68intergenic (‑2/+666)PP_3547/emrBputative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/multidrug efflux transport system ‑ permease subunit
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base C (4/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCCGGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAA  >  AE015451/4022084‑4022527
                                                                                                                                                                                                                               |                                                                                                                                                                                                                            
gccgGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCaa                                                                                                                                                                                                                   >  1:731607/1‑235 (MQ=255)
   ggTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGT                                                                                                                                                                                                                      <  1:1237211‑M1/229‑1 (MQ=255)
     tgatgaCAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcgc                                                                                                                                                                                                               <  2:979549/234‑1 (MQ=255)
               gACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGA                                                                                                                                                                                                                             >  1:934612‑M1/1‑210 (MQ=255)
                       ccATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGACCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGGCCGTCAAGCGATGAGGCCCCa                                                                                                                                                                                                     >  2:685033/1‑226 (MQ=255)
                           nggTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                                                                                                       <  2:517225/235‑1 (MQ=255)
                           gggTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                                                                                                       <  1:965515/236‑1 (MQ=255)
                             gTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                                                                                            <  2:522536/229‑1 (MQ=255)
                             gTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                                                                                            >  1:522539/1‑229 (MQ=255)
                             gTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                                                                                                       <  2:1078172/234‑1 (MQ=255)
                               ttGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTGCAgcgc                                                                                                                                                                                   >  1:834560/1‑236 (MQ=255)
                                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                                                                                                 <  2:548474/202‑1 (MQ=255)
                                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                                                                                                 >  1:548477/1‑202 (MQ=255)
                                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTc                                                                                                                                                                               >  1:457668/1‑236 (MQ=255)
                                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAT                                                                                                                                                                                                                            >  2:1194674‑M1/1‑191 (MQ=255)
                                            ggggAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                                                                                            <  1:374292/214‑1 (MQ=255)
                                            ggggAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCgg                                                                                                                                                                        <  1:658708/234‑1 (MQ=255)
                                                     aaCATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCg                                                                                                                                                             >  1:634388/1‑236 (MQ=255)
                                                     aaCATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGC                                                                                                                                                                                                                          >  2:956321‑M1/1‑175 (MQ=255)
                                                       cATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGc                                                                                                                                                               <  1:92425/232‑1 (MQ=255)
                                                               aTAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                                                                                                  >  1:590760/1‑205 (MQ=255)
                                                               aTAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                                                                                                  <  2:590757/205‑1 (MQ=255)
                                                                  ggCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCCTGGATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagccgcaatccgcg                                                                                                                                                                    >  2:787201‑M1/1‑168 (MQ=255)
                                                                      aaaTCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCggg                                                                                                                                                                       >  2:826717/1‑209 (MQ=255)
                                                                       aaTCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCTGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGa                                                                                                                                           <  1:566227/236‑1 (MQ=255)
                                                                        aTCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                                                                                                                                                                                             <  1:479082/169‑1 (MQ=255)
                                                                          cgcgCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                                                                                                  >  1:791534/1‑194 (MQ=255)
                                                                          cgcgCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                                                                                                  <  2:791531/194‑1 (MQ=255)
                                                                              cTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgca                                                                                                                                      <  2:1057358/234‑1 (MQ=255)
                                                                               tGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCGCGGagcgcagc                                                                                                                                    >  1:786173/1‑235 (MQ=255)
                                                                                gCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcagcg                                                                                                                                   <  1:789570/235‑1 (MQ=255)
                                                                                   ttGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAgg                                                                                                                               <  1:62250/236‑1 (MQ=255)
                                                                                    tGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGa                                                                                                                                                           >  2:632744/1‑207 (MQ=255)
                                                                                    tGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAgcgcagcgca                                                                                                                                 >  2:1201040/1‑233 (MQ=255)
                                                                                              aaCATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcagcgc                                                                                                                                  >  2:578285/1‑222 (MQ=255)
                                                                                               aCATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgag                                                                                                                   <  1:956324‑M1/236‑98 (MQ=255)
                                                                                                   cAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                                                                                                       <  1:968763/164‑1 (MQ=255)
                                                                                                   cAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                                                                                                       >  2:968760/1‑164 (MQ=255)
                                                                                                   cAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAgg                                                                                                                               <  2:871891/220‑1 (MQ=255)
                                                                                                            cAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGt                                                                                                       >  2:557943/1‑235 (MQ=255)
                                                                                                                 ctctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgaga                                                                                                                  >  1:858014‑M1/1‑121 (MQ=255)
                                                                                                                   ctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgaga                                                                                                <  1:733656‑M1/235‑117 (MQ=255)
                                                                                                                     ggttttgatctggcttttgatctggcTTTTGATCTTGATCTTGGTCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGa                                                                                                                                                                      <  1:271503/161‑1 (MQ=255)
                                                                                                                        tttgatctggcttttgatctggctTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgagagcga                                                                                            >  2:1005485‑M1/1‑114 (MQ=255)
                                                                                                                              ctggcttttgatctggcttttgatCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAT                                                                                                                                                                                                                            >  2:1139094‑M1/1‑100 (MQ=255)
                                                                                                                              ctggcttttgatctggcttttgatCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAT                                                                                                                                                                                                                            <  1:1139097‑M1/100‑1 (MQ=255)
                                                                                                                                   ttttgatctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAACGCAGGGGCCGGATGATGGGAGCGCAGCg                                                                                                        >  1:521546/1‑211 (MQ=39)
                                                                                                                                            ggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTg                                                                             <  2:834557/228‑1 (MQ=11)
                                                                                                                                              cttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCCAGCGTTTTTTGGTTACTTTTTGTCGCGGTTGACaaa                                                                        >  1:859570/1‑230 (MQ=11)
                                                                                                                                                        tgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCg                                                           <  1:588996/235‑1 (MQ=31)
                                                                                                                                                            cttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAgg                                                      >  1:1082078/1‑236 (MQ=31)
                                                                                                                                                              tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGGGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAg                                                                                                          <  1:604526/182‑1 (MQ=14)
                                                                                                                                                              tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCCcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgagagc                                                                                              <  1:1006520‑M1/194‑119 (MQ=255)
                                                                                                                                                                    tgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGaaa                                                <  2:634385/234‑1 (MQ=33)
                                                                                                                                                                              gACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGc                                    >  1:1163335/1‑236 (MQ=33)
                                                                                                                                                                                    aggaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgagagcgagagcgagagcatcgggatcttgccttgatttta                                                           <  1:1005488‑M1/207‑154 (MQ=255)
                                                                                                                                                                                      gaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgagagcgagagcgagagcatcgggatcttgccttgattttagctcaatgccccgcc                                            >  2:908129‑M1/1‑52 (MQ=255)
                                                                                                                                                                                              aCGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgagagcgagagcgagagcatcgggatcttgc                                                                     <  1:1061127‑M1/187‑144 (MQ=255)
                                                                                                                                                                                                gCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAacgg                   >  1:1204358/1‑235 (MQ=38)
                                                                                                                                                                                                gCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgagagcgagagcgagagcatcgggatcttgccttgattttagctcaatgccccgccgatgccgggcccgccttagcggta                    >  2:1129470‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                 cAGGCCTTGCGGAGGCAGGGGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAacg                    <  2:878681/233‑1 (MQ=37)
                                                                                                                                                                                                     ccTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTCTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGtt             <  1:176480/236‑1 (MQ=38)
                                                                                                                                                                                                          cGGAGGCAGGTGACGGGCATCCATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgagagcgagagcgagagcatcgggatcttgccttgattttagctcaatgccccgccgatgccgggcccgccttagcggtaaagggcggcttg        >  1:56832‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                           ggagggagGTGACGGGTATGGATGCCCGTCAAGAGCTGAGGCCCCAGGTTGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAAcc                                                                                                                                                <  2:1131021/99‑1 (MQ=11)
                                                                                                                                                                                                           ggAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAACGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGc                                                             <  2:521543/182‑1 (MQ=14)
                                                                                                                                                                                                               gCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCGTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCTGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCg                                                              >  1:442367/1‑177 (MQ=14)
                                                                                                                                                                                                                 gcttgacGGGCATGGATGCCCGTCAcctcccttcgcaagccctgcgctcggcctcctgaagtcgcaatccgcggcgcctgggctatcgcgcgcttagaagcaaaatcaaaagccgaagcgagagcgagagcgagagcgagagcgagagcgagagcatcgggatcttgccttgattttagctcaatgccccgccgatgccgggcccgccttagcggcaaagggcggcttggcgaac  >  1:845376‑M1/4‑25 (MQ=255)
                                                                                                                                                                                                                            aTGGATGCCCGTCAAGCGCTGAGGCCCCAGGAGGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGt                                                                                                       >  2:178609/1‑123 (MQ=17)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                            
GCCGGTGATGACAATGACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAA  >  AE015451/4022084‑4022527

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