Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,844,971 T→G 26.1% intergenic (+47/+72) apt → / ← recR Adenine phosphoribosyltransferase/Recombination protein RecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,844,9710TG26.1% 140.8 / 17.7 151intergenic (+47/+72)apt/recRAdenine phosphoribosyltransferase/Recombination protein RecR
Reads supporting (aligned to +/- strand):  ref base T (22/77);  new base G (28/10);  total (59/92)
Fisher's exact test for biased strand distribution p-value = 6.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

CCGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAT  >  ACB143/4844841‑4845102
                                                                                                                                         |                                                                                                                                         
ccGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCt                                                                                                                                        <  4:315143/134‑1 (MQ=255)
ccGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTAc                                                                                                                                      <  3:358089/136‑1 (MQ=255)
 cGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACg                                                                                                                                     <  5:70768/136‑1 (MQ=255)
 cGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACg                                                                                                                                     <  2:121890/136‑1 (MQ=255)
  ggAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCCCGGCTACg                                                                                                                                     >  2:216047/1‑135 (MQ=255)
  ggAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑AGCCGGTTAAACCCCCGCCTAc                                                                                                                                      >  2:159218/1‑135 (MQ=255)
   gAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGaa                                                                                                                                   <  7:374152/136‑1 (MQ=255)
         gACGGGTCGCGTCGGTTGGAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCcgc                                                                                                                           <  2:357320/136‑1 (MQ=255)
         gACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCcgc                                                                                                                           <  8:367774/136‑1 (MQ=255)
           cGGGTCGCGTGGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGt                                                                                                                        <  6:337382/136‑1 (MQ=255)
           cGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGt                                                                                                                        <  6:139926/136‑1 (MQ=255)
             gggcgggTCGGTGGCACGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTAt                                                                                                                      <  1:387440/130‑1 (MQ=255)
             ggTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAAACCGCGCCGACGGAAGCCGCGTAt                                                                                                                      >  7:253223/1‑136 (MQ=255)
                 gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTTGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAAACCCCTACATAGAAAAGCCGGGTGTGtt                                                                                                                  >  3:214100/1‑136 (MQ=255)
                 gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCAGCGAAGAAAGCCGCCTATGttt                                                                                                                 >  6:156764/1‑136 (MQ=39)
                 gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAAACCGCACCTACGAAAGCCGCGGATGtgt                                                                                                                 >  2:119748/1‑134 (MQ=255)
                    tCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAACCCGCGCCAACGCAAGCCGCGGATGtttttt                                                                                                              >  4:304887/1‑134 (MQ=39)
                     cgttggCAAGCGGCCGGGTTGCCTACCTTCTGCCTGACCGAGTTATCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTa                                                                                                             <  8:303748/131‑1 (MQ=255)
                     cGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTa                                                                                                             <  8:188273/136‑1 (MQ=255)
                      ggTTCCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTa                                                                                                             <  7:93118/135‑1 (MQ=255)
                       gTTGCAAGCGGCCGGGGGGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAgg                                                                                                           <  6:169223/136‑1 (MQ=255)
                           aaaGCGGACGGGGTGCCTACATTCTGCCTTACCGAGTTTTCCCTCAGCGCGTATTGATTGGGCATTTGGGGCTGAATTTCAGCACAAT‑‑CGCCGGTAAAACCCCGCCGAAGAAAACCGCGGATTTTTTTTATgcgcc                                                                                                      >  2:263309/2‑135 (MQ=255)
                                ggCCGGGGTGCCTACCTTCTGCCTGCCCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAACCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGttt                                                                                                 <  3:54352/136‑1 (MQ=255)
                                ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCTGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAAACCCCCGCACAGGAAACCCCCGTTTTTTTTTTAGGGGGGGGtt                                                                                                  >  6:135253/1‑136 (MQ=255)
                                ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTAGCGGGGGTAAACCCCTCGAACAAAAGCCGCGAGTTTTTTTTTGGCggggtgt                                                                                                  >  4:15870/1‑133 (MQ=255)
                                ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAAACCCCGCCGACCCAAGCCGCGGATGTTTTTTTGGGGCGGGgtt                                                                                                 >  4:141549/1‑136 (MQ=255)
                                ggCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑AGCGGGGTAACCCGCGCCCAACGAAGCCGCGGTTGTGTTTTTGGGGCGGGgtt                                                                                                 >  7:21234/1‑136 (MQ=255)
                                    ggggTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAACCCGCGCCGACGGAAGCCGCGGTTGTGTTGTAGGCGCGGGGTTAccc                                                                                             >  5:266157/1‑136 (MQ=38)
                                                      ctcaCCGAGTTTCCCCTCAGCGAGTATTGCGTGGGCTTTTGGGGCTGCTTTCCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTCCGAACGCCGCGTATGTTTTGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                           <  4:84109/133‑1 (MQ=38)
                                                       tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCCCTCCTACACAAACCGCCGATGTTTTTTATGCGCGGGTTTTACCCCCAAAGGgccgcccaccg                                                                          >  6:96373/1‑131 (MQ=25)
                                                        gACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAAACCGCTGCTACGAAAGCCGCGGATGTTTTGTTTGCGCGGGGTTTCCCGCGCAAGGGCCGCCaacagg                                                                         >  2:109702/1‑132 (MQ=37)
                                                        gACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAAACCCCGCCAACCAAATCCGCCGATGTTTTTTAAGCGCGGGTTTAACCGCCAAATGGCCGCCAACCgg                                                                         >  3:44475/1‑136 (MQ=255)
                                                          ccGAGTTTTCCCTCTCCGAGTATTGAGTGGGCTTGTGGGGCGGCTTGGCAGCCAATT‑‑CGCGGGTAAACCCGCTCCTACGAAACCCGCGTATGTTTTGTAGGCCCGGGTTTACCCGCGAATGGGCCGCCAAGCGGcc                                                                       <  2:279331/136‑1 (MQ=37)
                                                           cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGGATGTTTTGTAGGCGCCGGGTTACCCGCGAATGGGCCGCCCAGCGGccc                                                                      >  1:252567/1‑136 (MQ=255)
                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAACCCGCTCCGACGCAACCCGCGGATGGTTTTTAGGCGCGGGGTTTCCCGCGAATGGGGCCCCCAGCGGcccc                                                                     >  5:392701/1‑136 (MQ=37)
                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAACCCGCGCCAAACCAAGCCCCCGATGTTTTTTTTGCGCGGGTTTTCCCCCCAAATGGCCCCCCACCGGgccc                                                                     >  5:302619/1‑136 (MQ=16)
                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAACCCCCTCCGACGAAAGCCGCGTTTGTTTTTTAAGCGCGGGGTTAACCGCGAAAGGGCCGCCCAGCGGccc                                                                      >  1:48927/1‑135 (MQ=37)
                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAACCCCCGCCAACGAAAACCGCGGATTTTTTTTATGCGCGGGGTTTACCGCGAATTGGCCGCCCAGCGGgcc                                                                      >  5:363973/1‑135 (MQ=255)
                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAACCCCCCCCTAACAAAGCCCCGTATTTTTTTTTGGCGCGGGTTTACCCCCGAATGGGCCGCCAAACCGcccc                                                                     >  6:267263/1‑136 (MQ=37)
                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAAACCCCTCCGAACCAAGCCCCGGTTGTTTTTTTTGCCCCGGGTTTACCCCGAAAGGGGCCCCCACCGGcccc                                                                     >  3:165137/1‑136 (MQ=255)
                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAAACCGCTCCAACCAAACCCGCGGATGTTTTTTAGGCGCGGGTTTACCCCCGAATGGGCCGCCCACCCGcccc                                                                     >  8:188238/1‑136 (MQ=37)
                                                            gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCACATT‑‑CGCGGGTTAACCCGCTCCGACGAAAGCCGCGGATGTTTTGTAGGCGCGGGTTTTACCGCGCATGGGGCGCCCAGCCGcccc                                                                     >  5:91497/1‑136 (MQ=255)
                                                                 ttCCCTCAGCGAGTATGGAGTGGGCTTTTGGGGCGGCTTTGCAGCCCATT‑‑CGCGGGAAAACCCGCTCCTACGAAAGCCGCGTATGTTTGGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATg                                                                 <  5:103297/135‑1 (MQ=38)
                                                                          gtaGTTTTGAGTGGGCTTTGGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAACCCCGCTCCTACGAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCg                                                      <  6:264861/134‑1 (MQ=38)
                                                                          cGATTATTGCGGGGGCTTTTGGGGGTGCTTTGCAGCCCATT‑‑CGCGGGTAACCCCGAGCCTCCGAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCg                                                      <  8:253223/136‑1 (MQ=37)
                                                                           gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAAACCGCGGCTAAGAAAGCCGCGTATTTTTTTTATGCCCGGGGTTTCCCGCGGATTGGCCCCCAAGCGGCCCCCATGAATCAGAGCg                                                      >  8:80066/1‑135 (MQ=255)
                                                                           gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTTAAACCCCTGCGTCGCAAGCCGCGGTTGTTTTGTAGGCGCGGGGTTTCCCGCGAATGGGCCCCCCACCGGCCCCAATGAATAAGcgcg                                                      >  7:35531/1‑135 (MQ=255)
                                                                           gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGGATTTTTTGTTGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATAAGAGCg                                                      >  5:183453/1‑135 (MQ=255)
                                                                           gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTAACTAAGGCGCGTATGTTTTGTAGGAGCCGGTTTAACCGCCAATTGGGCGCCAAACCGCCCCAAAGAATCAGAGCg                                                      >  7:378865/1‑135 (MQ=25)
                                                                               aTTGAGTGCGGTTTTGGGGCGTCTTCCCACCAATTGCCCGGTAAAACCCGCTCAAACTAAAGCGGCTTTTGTTTTGTCGCCGGGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATg                                                 <  4:395567/136‑1 (MQ=255)
                                                                                    gTGGGCTTTTGGGGCGGCTGTGCCGCCCATT‑‑CGCGGGTAACCCCGCTCCTACGAAAGCCGCGTATGTTTTGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                            <  7:167221/136‑1 (MQ=38)
                                                                                     tGGGGTTTTGGGGCGGCTTTGCTGCCCATT‑‑GGGGGGTAAACCCGCGCCTACGAACGCGGCGTATGTTTTGATGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                            <  5:181469/135‑1 (MQ=37)
                                                                                     tGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCCGGGTAAACCCCGCCATCGCAAGCCCGGGTTTGTTTTTTGGCGCCGGGTTTGCCCCCCAAGGGGCCCCCACCGGGCCCAATTAATACGCGCGCAATTGGCCgg                                           >  8:156446/1‑136 (MQ=255)
                                                                                     tGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                            <  8:106378/135‑1 (MQ=255)
                                                                                     gtgggttttGGGGGGGCATTTGCCCCCAATC‑‑CGCGGTAAACCCCCCCCAACCCAACCCGGCGTTTTTTTTGAGGGCGGGGTTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                            <  5:32174/131‑1 (MQ=255)
                                                                                     gggggTTTGGGGGCTGGTTTCCCGCCCTTG‑‑CGCGGGTAACCCCGCTCCTACCAACGCCGCTTATTTTTTGTAGGGGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                            <  7:233511/134‑1 (MQ=25)
                                                                                     ggggCTTTGGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                            <  2:263323/134‑1 (MQ=255)
                                                                                      tgggttttGGGGCGTTCTTTGCTGCCCTTC‑‑CGCGGGTAACCCCGCTCTTCCGAACGCCGCTTTTGTTTTGTAGGGGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCgg                                           <  8:227611/132‑1 (MQ=25)
                                                                                       ggCTTTTGGGGCTGCTTTCCAGCCAATT‑‑CGCGGGAAACCCCGCTCCTACCAACCCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGgc                                          <  4:340065/135‑1 (MQ=38)
                                                                                           tttGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTGCTACGAAAGCCGCGGATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCt                                     >  3:293872/1‑136 (MQ=255)
                                                                                            ttGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGTAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGGATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCt                                     <  7:7712/135‑1 (MQ=255)
                                                                                                ggCTGCTTTGCAGCCCATT‑‑CGCGGGTTAACCCGCTGCTACGAAAGCCGCGGTTGTTTTGTAGGCGCGGGGTTACCCGCGCATGGGCCGCCCAGCGGCCCCAATGGATCAGAGCGATATGGGCCGGCGCCCGGCa                                  >  3:304530/1‑134 (MQ=37)
                                                                                                ggCTGCTTTGCAGCCCATT‑‑CGCGGGTTAAACCCCGCCGACGGAAGGCGCGGTTTGTTTTTTTGGGCGGGGTTTACCGCGGATTGGGCCCCCACCCGGCCCACAGAATGAGCGCGGAAATGGGCCGCGccccggaa                                 >  4:270630/1‑130 (MQ=255)
                                                                                                 gctgctTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAg                               <  1:309734/136‑1 (MQ=255)
                                                                                                 gctgctTTGCAGCCCATT‑‑CGCGGGAAAACCCGCGCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAg                               <  7:279688/136‑1 (MQ=255)
                                                                                                 gctgctTTCCAGCCCTTT‑‑CGCGGGTAACCCCGCTCCTCCGAAAGCCGCTTTTGTTTTGTCGGGGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAg                               <  8:36434/136‑1 (MQ=38)
                                                                                                 gctgctTGGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAg                               <  2:61475/136‑1 (MQ=255)
                                                                                                  ggggctTTGCAGCCCATT‑‑CGCGGGAAAACCCGCTCATACGAAAGCCGCGTTTGTTTTGAAGGCGCGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAg                               <  4:353100/133‑1 (MQ=255)
                                                                                                  gcggctTTGCAGCCCATT‑‑CGCGGCTAACCCCGATCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAg                               <  1:101953/133‑1 (MQ=255)
                                                                                                  ctgctTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACTAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAgg                              <  1:55620/136‑1 (MQ=255)
                                                                                                   gcggatttGCAGCCCATT‑‑CGCGGTAAACCCCGCTCCTCCGAACGCCGCGTATGTTTTGTCGGCGGGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAg                               <  7:154047/131‑1 (MQ=255)
                                                                                                      ctTTGCAGCCCATT‑‑CGCGGGTTAACCCGCGGCGACGGAAGCCGGGGTTGTTTTGTTGGCGCGGGGTTTCCCGCGAATGGGGCGGCCCGCGGGCCCCATGGATGAGGGGGGAATGGGGCCGCGCCCGGGAAAGGGAt                           >  8:263892/1‑136 (MQ=255)
                                                                                                        tttGCAGCCCATT‑‑CGCGGGTAAACCCGCGGCTACGAAAGGCGCGGTTGGTTTGTTGGCGCGGGGTTTCCCGCGAATGGGGCGCCCAGCGGGCCCCATGGATCAGAGCGCAATGGGCCGGCGCCCGGGCAAGGGATg                          >  6:374516/1‑136 (MQ=255)
                                                                                                        tttGCAGCCCATT‑‑CGCGGGAAACCCCGCTCCAACGAACGCCGCGTATGTTTTGTCGGCGCGGGTTTCCCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  8:121688/136‑1 (MQ=255)
                                                                                                        tttGCAGCCCATT‑‑CGCGGGAAAACCCGCTCCTACGAACGCCGCGTAGGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  7:314865/136‑1 (MQ=255)
                                                                                                        tttCCAGCCCATT‑‑GGCGGGTAAACCCGCTCCTACTAACGCCGCGTATGTTTTGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  5:386537/136‑1 (MQ=39)
                                                                                                        ttggCAGCCCATT‑‑CGCGGGTAAACCCGCGCCTACTAAAGCCGCGTATGTTTTTTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCAGCAAAGGCATg                          <  4:332060/133‑1 (MQ=38)
                                                                                                        gtcgctgCCCTTT‑‑GGCGGGTACACCCGCTCCAACGAAAGCCGCGTATTTTTTGTCGGCGGGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  4:210564/130‑1 (MQ=255)
                                                                                                        gcggCACCCCATG‑‑CGCGGTACAACCCGCGCATCCCAATGCCCCGTATGTTGTGACGGCGCGGTTATACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  8:333960/133‑1 (MQ=255)
                                                                                                         ttggcccccATT‑‑GGGGGGTAAACCCGGCTCCACCGAACGCCGCTTTTGTTTTGTAGGGGGGGTTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  3:197687/130‑1 (MQ=255)
                                                                                                         ttGCAGCCCATT‑‑CGCGGGGTAAACCGCGGCGTCGGAAGGCGGGGTTGTTTTGTTGGCGGGGGGTTTCCCGCGGATGGGGCGCCCAGCGGGCCCCATGGATGAGAGCGCAATGGGGCGGCGCCCGGGAAAgggaggg                         >  6:193721/1‑131 (MQ=255)
                                                                                                          tGCAGCCCATT‑‑CGCGGGTAACCCCGCTCCTACTAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  8:232609/134‑1 (MQ=255)
                                                                                                          tGCAGCCCATT‑‑CGCGGGTAACCCCGCGCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  1:136977/134‑1 (MQ=255)
                                                                                                          tGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  8:128026/134‑1 (MQ=255)
                                                                                                          tGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGAAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  7:18316/134‑1 (MQ=255)
                                                                                                          tGCAGCCCATT‑‑CGCGGGTAAACCCGCGCCTACGAAAGCCGCGTTTGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  2:8185/134‑1 (MQ=255)
                                                                                                          tGCAGCCCATT‑‑CGCGGGAAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATTGGCCGGCGCCCTGCAAAGGCATg                          <  2:256620/134‑1 (MQ=255)
                                                                                                          tGCACCCCATT‑‑CGCGGGCAAACCCGCTCCTACCAAAGCCGCTTATGTTTTGACGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  5:206888/134‑1 (MQ=39)
                                                                                                          ggCAGCCCATT‑‑CGCGGGAACACCCGCGCCTACGAAACCCGCGTATGTTTTGTCGGCGGGGTTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                          <  1:384515/133‑1 (MQ=255)
                                                                                                           gCAGCCCATT‑‑CGCGGGAAAACCCGCTCCTACGAACGCCGCGTATGTTTTGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATggg                        <  1:49092/135‑1 (MQ=255)
                                                                                                             aGCCCATT‑‑CGCGGGTAAACCGGCTCCTACGAACGCCGCGTATGTTTTGTCGGCGCGGGTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCa                     <  2:313617/136‑1 (MQ=255)
                                                                                                                    t‑‑CGGGGCTAAACCCGCTCCTACGAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAc               <  6:178916/135‑1 (MQ=255)
                                                                                                                       tcccgGGTACACCGGCTCCACCGACCGCCGGGTTTGTGTGGTCGGCGCGGGTATCCCGGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTa                <  4:102381/131‑1 (MQ=255)
                                                                                                                        gcgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccgc            <  2:287538/136‑1 (MQ=255)
                                                                                                                           gggTAAACCCGCTCCTACTAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccgccg          <  7:267054/136‑1 (MQ=255)
                                                                                                                            ggTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt         <  2:250873/136‑1 (MQ=255)
                                                                                                                             gttAAACCCGCTCCTACGAACGCCGCGTATGTTTTGAAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGt         <  1:220502/134‑1 (MQ=255)
                                                                                                                               tAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc       <  6:357226/136‑1 (MQ=255)
                                                                                                                                aaaCCCGCGCCTGCGAAAGGCGCGGATGTTTTGTAGGCGCGGGGTTTGCCGCGCATGGGGCGCCAAGCGGCCCCAATGAATCAGAGGGAAATGGGCCGGCGCCCTGCAAAGGCAGGGGCCCGGGTGCCGCCCTccc      >  2:68144/1‑135 (MQ=255)
                                                                                                                                    ccGCTCCTACGAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAt  <  1:345034/136‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
CCGGAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATT‑‑CGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCAT  >  ACB143/4844841‑4845102

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