Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,945,329 G→T 7.5% intergenic (+689/+109) PP_5673 → / ← flhB putative D‑alanine‑‑D‑alanine ligase/flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,945,3290GT 7.5% 182.5 / 6.5 53intergenic (+689/+109)PP_5673/flhBputative D‑alanine‑‑D‑alanine ligase/flagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base G (36/13);  new base T (2/2);  total (38/15)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

GGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAG  >  AE015451/4945094‑4945563
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
gggCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCg                                                                                                                                                                                                                                            <  1:1180008/236‑1 (MQ=255)
                           gCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCgg                                                                                                                                                                                                                                           <  1:1162146/210‑1 (MQ=255)
                                            tCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGa                                                                                                                                                                                                                               >  1:951271/1‑205 (MQ=255)
                                                       tCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCCTCGCGAGCAAGCTCGCTCCTCCATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGcac                                                                                                                                                                                                                                       <  1:212841/186‑1 (MQ=255)
                                                                      tGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCgg                                                                                                                                                                                                                                           <  1:731056/167‑1 (MQ=255)
                                                                                        tgAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCg                                                                                                                                                                                                                                            <  1:378929/148‑1 (MQ=255)
                                                                                                       cccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGcaca                                                                                                                                                                                                                                      <  1:508138/139‑1 (MQ=255)
                                                                                                                        gAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGaaa                                                                                                                                                                                                                    <  1:1018201/140‑1 (MQ=255)
                                                                                                                                  gctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGaa                                                                                                                                    >  1:544768/1‑210 (MQ=255)
                                                                                                                                    tcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCg                                                                                                                                                                                                                                            <  2:951268/104‑1 (MQ=255)
                                                                                                                                           tttGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGAt                                                                                                  >  1:527280/1‑235 (MQ=255)
                                                                                                                                               tttCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGc                                                                                              >  1:968531/1‑235 (MQ=255)
                                                                                                                                                 tCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGGTCTCACAGGCAACGCAAACCTATCTCCGAACAGATCAACTGTCGCGGCGCAg                                                                                                              >  1:77221/1‑217 (MQ=255)
                                                                                                                                                 tCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTcc                                                                                                          >  1:772082/1‑221 (MQ=255)
                                                                                                                                                   gggCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGgtc                                                                                         >  2:872756/1‑236 (MQ=255)
                                                                                                                                                       cAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGGCCGGCGGGGGCGGCgggg                                                                                          >  1:106783/1‑230 (MQ=255)
                                                                                                                                                          ttATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCa                                                                                             <  2:527277/225‑1 (MQ=255)
                                                                                                                                                           tATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAAc                                                                                                                           >  1:463168/1‑194 (MQ=255)
                                                                                                                                                                      gATTTACGCCCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCgcggcg                                                                                                                 >  2:522970/1‑193 (MQ=255)
                                                                                                                                                                      ccggcgcgcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCGCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                                                                         >  1:18925/7‑153 (MQ=255)
                                                                                                                                                                       aTTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTg                                                                              >  1:862412/1‑227 (MQ=255)
                                                                                                                                                                       aTTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCCCAGGTCCGGCGGGATCGGCAGgt                                                                                          >  1:9903/1‑215 (MQ=255)
                                                                                                                                                                              cgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCgcggcg                                                                                                                 >  1:1156188/1‑185 (MQ=255)
                                                                                                                                                                                cgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCCCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCATAGGCAACGCAAAACTATCGc                                                                                                                                        >  1:1213459/1‑160 (MQ=255)
                                                                                                                                                                                 gcccGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCaa                                                                                                                                                      <  1:285703/143‑1 (MQ=255)
                                                                                                                                                                                   ccGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCCCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGGTCAACTGTCGCGGCGCAGGTCcggcgg                                                                                                     >  1:575275/1‑192 (MQ=255)
                                                                                                                                                                                                    tCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGaa                                                                                                                                    >  2:163130/1‑144 (MQ=255)
                                                                                                                                                                                                         atTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGACCGGCGGGATCGGCCGGGCGCCCATCacgg                                                                              >  2:252078/1‑190 (MQ=255)
                                                                                                                                                                                                              gcCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAac                               >  1:767969/1‑235 (MQ=255)
                                                                                                                                                                                                                cGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGgagc                                                                                                                                                                                                       >  2:423058/1‑65 (MQ=255)
                                                                                                                                                                                                                cGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCTGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGacc                               <  2:767966/233‑3 (MQ=255)
                                                                                                                                                                                                                cGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGctc                                                                          >  2:809997/1‑190 (MQ=255)
                                                                                                                                                                                                                cGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAacac                             >  1:1233890/1‑235 (MQ=255)
                                                                                                                                                                                                                     aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGga                                                                                                                                                                                                         >  2:342180/1‑58 (MQ=255)
                                                                                                                                                                                                                     aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATgcg                                                                                                                                                                                     >  2:1078697/1‑78 (MQ=255)
                                                                                                                                                                                                                     aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTAtt                                             >  1:711893/1‑214 (MQ=255)
                                                                                                                                                                                                                     aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGGATCTGGAACACATACg                        >  1:562530/1‑235 (MQ=255)
                                                                                                                                                                                                                     aaaaaaGGCGGTCGCGCGCGCCGGGCACAGGAGGAACTGGCCCGAAACAAATGTAGgagc                                                                                                                                                                                                       <  2:136697/60‑1 (MQ=37)
                                                                                                                                                                                                                          aGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgacgcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacct                                                                                                                               >  2:250479‑M2/1‑60 (MQ=255)
                                                                                                                                                                                                                           ggCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGa                                                                                                                                                                                                                      <  1:773176/39‑1 (MQ=255)
                                                                                                                                                                                                                              agtCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaat                                                    <  2:1213167‑M2/196‑142 (MQ=255)
                                                                                                                                                                                                                              agtCGCGCGCGACTGCTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaat                                                    >  1:1213170‑M2/2‑56 (MQ=255)
                                                                                                                                                                                                                               gTCGCGCGCGCCGGGTACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAgg                                                                                                             >  2:583177/1‑140 (MQ=255)
                                                                                                                                                                                                                               gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGa                                                                                                   >  2:518595/1‑150 (MQ=255)
                                                                                                                                                                                                                               gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGACTCTGGAACACATACGc                       >  1:855948/1‑226 (MQ=255)
                                                                                                                                                                                                                               gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTg              >  1:954389/1‑235 (MQ=255)
                                                                                                                                                                                                                               gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGAGCAACTGTCGCGGCGCAGGCCCGGCgggg                                                                                                   >  2:195053/1‑149 (MQ=255)
                                                                                                                                                                                                                               gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCCCCGAACAGATCAACTGTc                                                                                                                       >  2:1131792/1‑130 (MQ=255)
                                                                                                                                                                                                                                 cgcgcgcgCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                                            <  2:630598/91‑1 (MQ=255)
                                                                                                                                                                                                                                  gcgcgcgaCTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaa                                         <  1:176703‑M2/204‑153 (MQ=255)
                                                                                                                                                                                                                                   cgcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGaa                                                                                                                                    <  2:641182/113‑1 (MQ=255)
                                                                                                                                                                                                                                    gcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATgcg                                                                                                                                                                                     >  2:227079/1‑63 (MQ=255)
                                                                                                                                                                                                                                    gcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACAt                           >  2:896539/1‑217 (MQ=255)
                                                                                                                                                                                                                                    gcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCCCCTg              >  2:1244458/1‑230 (MQ=255)
                                                                                                                                                                                                                                      gcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgcc                                                                                                                                                                                   <  1:996533/63‑1 (MQ=255)
                                                                                                                                                                                                                                      gcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGcgcc                                                                                                                                                                                   >  2:996530/1‑63 (MQ=255)
                                                                                                                                                                                                                                       cgcCGGGCACAGGAGGGACTGGCCCGa                                                                                                                                                                                                                      <  2:200473/27‑1 (MQ=255)
                                                                                                                                                                                                                                          cGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAg  <  1:872759/236‑1 (MQ=255)
                                                                                                                                                                                                                                           gggCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGt                                                                                                                        >  2:1255456/1‑117 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
GGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAG  >  AE015451/4945094‑4945563

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