Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,945,370 C→T 24.2% intergenic (+730/+68) 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,3700CT24.2% 117.7 / 34.3 62intergenic (+730/+68)PP_5673/flhBputative D‑alanine‑‑D‑alanine ligase/flagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (37/10);  new base T (6/9);  total (43/19)
Fisher's exact test for biased strand distribution p-value = 8.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGAT  >  AE015451/4945137‑4945581
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                     
aTCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCt                                                                                                                                                                                                                      <  1:538452/235‑1 (MQ=255)
                                                                                       gctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTc                                                                                                                                                                                      >  1:12338/1‑178 (MQ=255)
                                                                                       gctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGa                                                                                                                                                                                    >  2:696418/1‑180 (MQ=255)
                                                                                       gctcCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCgcggc                                                                                                                                    >  1:912334/1‑228 (MQ=255)
                                                                                                      tCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCcgg                                                                                                                          >  1:983576/1‑223 (MQ=255)
                                                                                                      tCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCACGAAACAAATGGAGGAGCGAGCTTGCt                                                                                                                                                                                                                  >  2:241629/1‑137 (MQ=255)
                                                                                                      tCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCCGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTcg                                                                                                                                                                                                                >  2:760429/1‑139 (MQ=255)
                                                                                                                 aTTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTc                                                                                                                                                                                      >  2:178373/1‑152 (MQ=255)
                                                                                                                       gCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTg                                                                                                                                           <  2:383165/189‑1 (MQ=255)
                                                                                                                       gCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTg                                                                                                                                           >  1:383168/1‑189 (MQ=255)
                                                                                                                           gATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactccaaattcaa                                                                                                                   <  1:62589‑M2/210‑99 (MQ=255)
                                                                                                                           gATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccacc                                                                                           <  1:19803‑M2/234‑123 (MQ=255)
                                                                                                                           ccggcgcgcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTcgc                                                                                                                                       >  1:440725/7‑189 (MQ=255)
                                                                                                                            cggcgcgcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGAt                                                                                                                                                 >  2:1029228/6‑178 (MQ=255)
                                                                                                                            aTTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGCGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGc                                                                                                                                                                          >  1:7212/1‑153 (MQ=255)
                                                                                                                            aTTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCg                                                                                                                  >  2:862622/1‑209 (MQ=255)
                                                                                                                                gcgcgcgACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCcggc                                                                                                                         >  2:584707/2‑198 (MQ=255)
                                                                                                                                      gACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgccctt                                                                             >  2:1009553/1‑236 (MQ=255)
                                                                                                                                          gCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGggcggc                                                                                                        >  1:269280/1‑203 (MQ=255)
                                                                                                                                                        cTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCcggcgg                                                                                                                       >  2:231915/1‑176 (MQ=255)
                                                                                                                                                        cTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTAtt                                                               <  1:178375/232‑1 (MQ=255)
                                                                                                                                                        cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcggagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtcccctcaaaa                                                                                                                        <  2:125031‑M2/176‑94 (MQ=255)
                                                                                                                                                         tCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGa                                                          >  1:433894/1‑236 (MQ=255)
                                                                                                                                                           gATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAAt                                                        >  1:748517/1‑236 (MQ=255)
                                                                                                                                                             tatTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCt                                                      <  1:231918/236‑1 (MQ=255)
                                                                                                                                                               ttGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGaa                                                                                                                    >  2:1073901/1‑171 (MQ=255)
                                                                                                                                                                  cgtcgtACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccc                                                        >  1:125033‑M2/1‑73 (MQ=255)
                                                                                                                                                                     cGTACAAACAAGGCGGTCGCGCGCGCCGGGCCCAGGAGGGACTGGCCCGAAACACATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGcca                                                                                                                                                                        >  2:34739/1‑112 (MQ=255)
                                                                                                                                                                          aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAgg                                                                                                             <  2:957746/168‑1 (MQ=255)
                                                                                                                                                                          aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAgg                                                                                                             >  1:957749/1‑168 (MQ=255)
                                                                                                                                                                          aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGc                                         >  1:859023/1‑236 (MQ=255)
                                                                                                                                                                          aaacaaGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCCCGATCTCACAGGCAACGCAAACCTATCGCCGCACAGATc                                                                                                                                                >  2:168936/1‑133 (MQ=255)
                                                                                                                                                                                    gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGCACAGATCAACTGTCGCGGCGCAgggc                                                                                                                             >  2:123138/1‑140 (MQ=255)
                                                                                                                                                                                    gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGgtcgtc                                                                                                        >  1:334892/1‑163 (MQ=255)
                                                                                                                                                                                    gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCtgc                                                                                                    >  2:460077/1‑165 (MQ=255)
                                                                                                                                                                                    gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCg                                                           >  2:700373/1‑208 (MQ=255)
                                                                                                                                                                                    gTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTgc                               >  1:741178/1‑236 (MQ=255)
                                                                                                                                                                                        cgcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCCGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTgcg                              >  2:671869/1‑233 (MQ=255)
                                                                                                                                                                                         gcgcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTg                                                     >  2:673408/1‑209 (MQ=255)
                                                                                                                                                                                           gcgcCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCAcc                         <  1:696421/235‑1 (MQ=255)
                                                                                                                                                                                             gactggtgcaggagggacTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcctg                        >  1:415166‑M2/9‑46 (MQ=255)
                                                                                                                                                                                              ccGGGCACAGGAGGGGCTGGCCCGAAACAAATGTAGGAGCGCGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGAt                                                                                                                                                                                   <  2:981405/78‑1 (MQ=255)
                                                                                                                                                                                                         aggGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGccc           >  1:1037099/1‑235 (MQ=255)
                                                                                                                                                                                                         aggGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGGCCTTCAGTGGct                                                                                             >  1:167600/1‑153 (MQ=255)
                                                                                                                                                                                                         aggGACTGGCCGAAACACACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaa                                           >  1:962472‑M2/1‑34 (MQ=255)
                                                                                                                                                                                                          gggACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttact                               <  2:415163‑M2/215‑183 (MQ=255)
                                                                                                                                                                                                          gggACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaat                                                                        >  2:679278‑M2/1‑33 (MQ=255)
                                                                                                                                                                                                          gggACTGGCCGAAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaat                                                                        <  1:679281‑M2/174‑142 (MQ=255)
                                                                                                                                                                                                            gACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCg                                                                                                                                                        <  2:338673/91‑1 (MQ=255)
                                                                                                                                                                                                            gACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCg                                                                                                                                                        >  1:338676/1‑91 (MQ=255)
                                                                                                                                                                                                              cTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGcttgccc                                                                               >  2:891259/1‑162 (MQ=255)
                                                                                                                                                                                                               tGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCt                                                      <  2:334889/186‑1 (MQ=255)
                                                                                                                                                                                                                ggCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGc                                                                                                                >  1:473790/1‑127 (MQ=255)
                                                                                                                                                                                                                 ggccgaAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaccacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaa                                                                                                    <  2:51421‑M2/138‑114 (MQ=255)
                                                                                                                                                                                                                 ggccgaAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcacaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactc                           <  2:401835‑M2/211‑187 (MQ=255)
                                                                                                                                                                                                                 ggccgaAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactc                                    <  2:377346‑M2/202‑178 (MQ=255)
                                                                                                                                                                                                                 ggccgaAACAAACGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactc                                    >  1:377349‑M2/3‑27 (MQ=255)
                                                                                                                                                                                                                   ccGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAgggc                                                                                                           >  2:731823/1‑127 (MQ=255)
                                                                                                                                                                                                                     gAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGa                                                          >  1:1005687/1‑176 (MQ=255)
                                                                                                                                                                                                                      aaacaaacGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatcc                                                                      <  1:570717‑M2/165‑144 (MQ=255)
                                                                                                                                                                                                                      aaacaaacGTAGGAGCGGGTTTacccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatcc                                                                      >  2:570714‑M2/1‑22 (MQ=255)
                                                                                                                                                                                                                             tGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGAt  <  1:903546/224‑1 (MQ=255)
                                                                                                                                                                                                                                    gcgagcTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGAt  >  2:872094/1‑217 (MQ=255)
                                                                                                                                                                                                                                     cgagcTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGgtcg                                                                                                          <  2:108554/112‑1 (MQ=255)
                                                                                                                                                                                                                                        gcTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGc                                         >  2:534887/1‑174 (MQ=255)
                                                                                                                                                                                                                                         cTTGCTCGCGATGCGCCGCGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTg                                >  2:394306/1‑184 (MQ=255)
                                                                                                                                                                                                                                         cTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGGTCGGCAGgt                                                                                                            >  2:243327/1‑106 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                     
ATCTGCTGTGAATCGGGGTGTTCGCCGTGAGGTTTTTGTCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCATCGCGAGCAAGCTCGCTCCTACATTTGTTTCGGGCCAGTTATTCCTGCGGGATTTACGCGCGACCGCCTTCGCGTACACCTCGATATTGCGCCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGACTGGCCCGAAACAAATGTAGGAGCGAGCTTGCTCGCGATGCGC‑‑CGCGCGGGCGGCGCTCGATCTCACAGGCAACGCAAACCTATCGCCGAACAGATCAACTGTCGCGGCGCAGGTCCGGCGGGATCGGCAGGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGAT  >  AE015451/4945137‑4945581

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