Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,945,333 A→G 7.5% intergenic (+693/+105) 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,3330AG 7.5% 165.7 / 4.2 53intergenic (+693/+105)PP_5673/flhBputative D‑alanine‑‑D‑alanine ligase/flagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (38/11);  new base G (2/2);  total (40/13)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

TCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAG  >  AE015451/4945138‑4945563
                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
tCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGa                                                                                                                                                                                                                               >  1:951271/1‑205 (MQ=255)
           tCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCCTCGCGAGCAAGCTCGCTCCTCCATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGcac                                                                                                                                                                                                                                       <  1:212841/186‑1 (MQ=255)
                                                           cccgcGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGcaca                                                                                                                                                                                                                                      <  1:508138/139‑1 (MQ=255)
                                                                            gAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGaaa                                                                                                                                                                                                                    <  1:1018201/140‑1 (MQ=255)
                                                                                      gctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGaa                                                                                                                                    >  1:544768/1‑210 (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)
                                                                                                                                                                                                 gCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCa                                                                                             >  2:1019562/1‑142 (MQ=255)
                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
TCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAG  >  AE015451/4945138‑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: