Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,772,546 T→G 45.0% intergenic (‑140/+21) grsB ← / ← lgrB_1 Gramicidin S synthase 2/Linear gramicidin synthase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,772,5460TG45.0% ‑3.5 / 34.7 86intergenic (‑140/+21)grsB/lgrB_1Gramicidin S synthase 2/Linear gramicidin synthase subunit B
Reads supporting (aligned to +/- strand):  ref base T (3/42);  new base G (31/3);  total (36/50)
Fisher's exact test for biased strand distribution p-value = 3.37e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

CATCTCCTACCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACAC  >  ACB143/4772412‑4772668
                                                                                                                                        |                                                                                                                          
cATCTCCTACCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGt                                                                                                                            >  2:97454/1‑135 (MQ=255)
cATCTCCTACCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGt                                                                                                                            >  3:242729/1‑135 (MQ=255)
      cTACCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGGCCGGGCATagcc                                                                                                                     >  5:125432/1‑133 (MQ=255)
      cTACCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGGCCCGGCGTACCg                                                                                                                     >  4:341766/1‑136 (MQ=255)
      cTACCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGGCATaacc                                                                                                                     >  7:299817/1‑133 (MQ=255)
      cTACCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGGCATAACg                                                                                                                     >  4:331883/1‑136 (MQ=255)
      cTACCAGGGCACCGGTATTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGTGCCGGGCATaacc                                                                                                                     >  7:276163/1‑133 (MQ=255)
        aCCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGCAAACCcgg                                                                                                                   >  1:135259/1‑136 (MQ=255)
           gggccACCGGTCTTTGCCTGCGG‑AACGTGCCCGCCCGCGACAGGCTCATCTTGCTGCCTAAAAAATGAACGCGCCGGACGGCTGCCGGCTCTGGGTGTACTGACGTTGACGACCC‑GGGCCCGGCCATACCGGGccc                                                                                                                <  8:315622/132‑1 (MQ=255)
           cgggCACCGGGCTTTGCCTGCGG‑ACGGTCCCAGCCCGCGACAGGCTCTGGATCCTTCCCAAAAAATCAACGAGCCGGACGGCGGCCGGCTCTGGGTGTACTCAAGTTGACGAGCC‑GGGCCCGGTAATACCGGGccc                                                                                                                <  2:305632/135‑1 (MQ=255)
           aGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGGCCGGGCCTACCGGGccc                                                                                                                >  5:157364/1‑136 (MQ=255)
           aGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGccc                                                                                                                <  3:331883/136‑1 (MQ=255)
              gccccGGTCTTTCCCTGCGG‑ACCGTGCCAGCCCGCGACAGGCTCATGTTGCTGCCGAACAAATGAAAGACCCGGACGGCTGCCGGCTGTGGGTGTACTCAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCa                                                                                                             <  6:177281/133‑1 (MQ=255)
              gccccGGTCTTTCCCTCCGC‑AACGTCCCAGCCCGCGACAGGCTCATGATCCTCCCCAACAAATCAAAGCCCCGCACGGCGGCCGGCTCTGGGTGTACTCAAGTTGACGAGCC‑GGGCCCGGTAATACCGGGCCCGCa                                                                                                             <  6:129206/133‑1 (MQ=255)
              gCACCGGGTTTTGCCTGCGA‑ACCGTGCCCGCCCGCAACAGGCTCATGATCCTCCCGAACAAAGGAACCAGCGGGCCGGCTGCCGTCTCTGTGTGTAGTGAAGTTGAGGAGCC‑GGGCCCGGTCATACCGGGCCCGCa                                                                                                             <  8:177616/136‑1 (MQ=255)
              gCACCGGGCTTTGCCGGCGG‑AACGGGCCAGCCCGCGACAGGCGCATGAGGCTGCCGACCAACTGACCGCGCCGGACGGCTGCCGGCTCTGGGGGTCCGGACGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCa                                                                                                             <  6:56544/136‑1 (MQ=255)
                aacGGTCTTTGCCTGCGG‑AAAGTTCCAGCCCGCGACAGGATCATGATGCTGCCGAACAAATGAAAGAGACGGACGGCTTCCGTCTATGGGTGTAATGAAGTTGACGAGCC‑GGGGCCGGTGAGAACGGGacccg                                                                                                              >  7:287746/3‑131 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGGCCCGGGCTCACCGGGccccccg                                                                                                            >  6:20421/1‑132 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGGCCCGGGCGGACCGGcccccc                                                                                                              >  2:17901/1‑132 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGGCCCGGCAATACC‑GGCCCGCcga                                                                                                           >  3:156912/1‑136 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGGCCGGGGATACCCGGccccccca                                                                                                           >  5:160281/1‑131 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGGCCCGGGATAACCGGccccccgc                                                                                                           >  7:42225/1‑131 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGGCCCGGCATACCGGGCCCGccgc                                                                                                           >  5:364272/1‑133 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGGCTACCGGGccccccga                                                                                                           >  2:243425/1‑131 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGGCTACCCGGCCCGccgc                                                                                                           >  2:278982/1‑133 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGGATACCGGGCCCGCagc                                                                                                           >  7:101339/1‑135 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGCCATCCGGGcccccc                                                                                                             >  7:110406/1‑131 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGCATAACCGGCCCGCa                                                                                                             >  7:68434/1‑134 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGGGCTACCGGGCCCGccgc                                                                                                           >  6:156738/1‑133 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGGGATACCCGGcccccccc                                                                                                           >  1:96402/1‑131 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGGCCAACCCGGcccccccc                                                                                                           >  8:217237/1‑131 (MQ=255)
                aCCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCAGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGCATAGCGGGCCCCCagc                                                                                                           >  2:302615/1‑135 (MQ=255)
                   gggcTTTGCCTGCGC‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGCGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAgcg                                                                                                        <  7:243693/133‑1 (MQ=255)
                   ggTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGGATACCGGGCCCCCCGAgcg                                                                                                        >  2:15768/1‑136 (MQ=255)
                   ggTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGGAGAACGGGGCCGCcgcgcg                                                                                                        >  1:13751/1‑136 (MQ=255)
                   ggTCTTTGCCGGCGC‑ACCGGGCCAGCCCGCGACAGGCTCATGATGCTGCCTAACAAATCAAAGCGCCGGCCGGCTGCGGGCTGTGGGTGTACTCACGTTCACGAGCC‑GGGCCCGGTCAAACCGGGCCCGCAGAgcg                                                                                                        <  8:108524/136‑1 (MQ=255)
                   ggTCTTTGCC‑GCGGCACCGGCCCCGCCCGCGACAGGCTCATGATGCTGCCGAACAAATCAAAGCGCCGCCCGGCGGCCGGCTGTGGGTGTACTGAAGTTGACGAGCC‑GGGCCGGGTAATCCCGGGCCCGCAGAgcg                                                                                                        <  8:76752/136‑1 (MQ=255)
                   ggTCTTGCCCTGCGC‑AACGTCCCCCCCCGCGACAGGCCCAGGAGCCGGCCAAACAAATGAACCCCCCGCACGGCTGCCGGCTCTGGGGTTACTGACGTTGCCGCGCC‑GGGCCCGGTAATACCGGGCCCGCAGAgcg                                                                                                        <  6:269121/136‑1 (MQ=255)
                    gTCTTTCCCGGCGG‑ACCGTGCCAGCCGGCGACAGGCTCTTCATGCTGCCCAACAAATGACAGCGCCGGCCGGCTGCCGGCTCTGGGGGTCCTTAAGTAGACGCGCC‑GGGCCCGGTCAAACCGGGCCCGCAGAgcgc                                                                                                       <  7:67804/136‑1 (MQ=255)
                      cattGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGCCAGAGCGGGCGCCCAGAccgc                                                                                                       >  3:94565/3‑134 (MQ=255)
                       tttGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGGCAAGCCCGGCCCGCcgcgcgggt                                                                                                     >  6:273465/1‑132 (MQ=255)
                            cTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGGCCCGGGATAACCGGCCCCCCCAGCGCGACggggg                                                                                               >  1:11546/1‑135 (MQ=255)
                                 gtaccGTGCCCGCCGGCGACTGGCTCTTGATGCTCCCAAACAAAGGAACGCGCGGGCGGGCGGCGGGCTCTGGGTGTACTGACGTTGACGACCG‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAAc                                                                                           <  6:186366/132‑1 (MQ=255)
                                   aaCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGGGCGGGGGTACAGGGCCCCCCGCGCGCGGCGGGGTCGAc                                                                                           >  7:349190/1‑134 (MQ=255)
                                   aaCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATTACAGGCCCGCAGAGCCCTAATGGGGGAACtc                                                                                         >  8:128890/1‑136 (MQ=255)
                                    aCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACtcg                                                                                        <  7:226236/136‑1 (MQ=255)
                                      gggCCAGCCCGCTACAGGTTCATTATGCTGCCTAAAAAATGAAAGAGCGGGACGGGTGCCGGCTCTGGGTTTACTGAAGTTCACGAGCC‑GGGCCGGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACtcgtc                                                                                      <  3:344996/134‑1 (MQ=255)
                                      gTGCCCCCCCGCGACAGGCTCATGAGGATGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGCCGCGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACtcgtc                                                                                      <  6:224576/136‑1 (MQ=255)
                                      gTGCCACCCCGCCACCGGCTCTTTATGCTCCCCAACAAATCAAAGCGCCGGCGGGCTGCCGGCTGGGGGTGTCCTCACGTTCACGCCCC‑GGGCCCGGGCATCCCGGGCCCGCAGAGCGCTACTGGGTCAACtcgtc                                                                                      <  8:160900/136‑1 (MQ=255)
                                       gcccggCCCGCGACAGGCTCATGATGCGGCCGAACAAATGAAAGAGCCGGACGGCGGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGCATACCGGGCCCGCAGAGCGCTACTGGGTCAACtcgtc                                                                                      <  7:294183/130‑1 (MQ=255)
                                        cccAGCCCGCGACAGGCTCTTGATGCTGCGGAACAAATGAAAGACCCGGACGGCTGCCGGCTCTGGGTGTACTGACGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACtcgtc                                                                                      <  8:79791/133‑1 (MQ=255)
                                          cAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGGATACCGGGCCCCCCGCGCGCGACTGGGGGAACTCGTCGAtc                                                                                  >  1:243019/1‑136 (MQ=255)
                                           cgCCGGCGACAGGCTCTTGATGCTCCCTAAAAATTCAAAGCGCCGCAGGGGTGCGGGCTGTGGGTTTCCTTATGTGGACTCGCC‑GGGCCCGGAAATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAtct                                                                                 <  2:16589/135‑1 (MQ=255)
                                           cgCCCGCGACAGGCTCATGATGCTGCCGAACAAATCAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGAtct                                                                                 <  1:69658/135‑1 (MQ=255)
                                                      ggCTCATGTTGCGGCCTAAAAATTCAACGCGCGGGAGGGGTGCCGGGTGTGGGTTTACTGACGATGCCGCCCC‑GGGCCCGTTCATACGGGCCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGt                                                                       <  4:163258/135‑1 (MQ=255)
                                                                     ctaaCAAATGAACGCGCCGGACGGCTGCGGGCTCTGGGTGTACTTAAGTTGACGCCCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcg                                                       <  6:366786/134‑1 (MQ=255)
                                                                     cGAACAAATTAACGCCCCGGACGGCGGCGGGGTCTGGGTGTACTGAAGTTGACCAGCC‑GGGCCGGGTCATCCCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcg                                                       <  7:108803/136‑1 (MQ=255)
                                                                     cGAACAAATGAACGAGCCGGCCGGGGGCGGGCTCTGGGGGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcg                                                       <  8:30858/136‑1 (MQ=255)
                                                                     cGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGGATACCGGGCCCCCCGAGCGCGACTTGGGGAACTCGGAGATGTCCGCTGCAGTTAGAAGGTGCAAggcg                                                       >  3:156689/1‑136 (MQ=255)
                                                                      gAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTGAACTAGGAGATATACTCTGCAGTTATACGTTTAAAtgcgg                                                      >  6:115106/1‑136 (MQ=255)
                                                                      gAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGGCATACCGGGCCCGCCGACCGCGACGGGGTGAACTACTGGCTCTACGCTGGAGGAAGACGTTTGAAggcgg                                                      >  4:151098/1‑136 (MQ=255)
                                                                      gAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGGCAGAACGGGCCCGCCGAGCGCGACAGGGGGACCTCGTCGAGCGCCACTGGACTAAGACGGTTGAAggcgg                                                      >  2:126833/1‑136 (MQ=255)
                                                                        aaaattGAACGCGCCGGAGGGCTGCGGGCTCTGGGGTTACTCACGTGCCCGAGCG‑GGGCCCGGTAATACGGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcgg                                                      <  4:130920/129‑1 (MQ=255)
                                                                        aCAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCCGGGCTACCCGGCCCCCACCCCCCGACTGGGGGCACGCGGCGCTCGCCGCGGCCGTAAGCCGGTTAAAGGCGGcc                                                    >  7:125397/1‑136 (MQ=255)
                                                                                    gCCGGAGGGGGGCCGGCTCTGGGTTCACTTAAGTTCCCGCGCG‑GGGCCCGGTCATCCCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGcc                                        <  2:135214/136‑1 (MQ=255)
                                                                                             cTGCCGGCTCTGGGTGTACTGAAGTGGACGCGCC‑GGGCCCGGTAATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTcc                               <  1:63899/136‑1 (MQ=255)
                                                                                                  ggCTCTGGGTGTACTGAAGTTGACGCGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAAt                          <  1:3778/136‑1 (MQ=255)
                                                                                                  ggCTCTGGGTGTACTGAAGTTGACGCGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAAt                          <  4:93543/136‑1 (MQ=255)
                                                                                                  ggCGCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAAt                          <  7:123907/136‑1 (MQ=255)
                                                                                                   gcgcggggTGAACTTAAGTGGACGACCC‑GGGCCCGGAAATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATc                         <  4:166478/131‑1 (MQ=255)
                                                                                                   gCTCTGGGTGTACTGAAGTTTACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGAAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATc                         <  6:289641/136‑1 (MQ=255)
                                                                                                        ggggGTACTGAAGTTGCGGCGCC‑GGGCCGGGCAATCCGGGCCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgc                    <  8:158194/136‑1 (MQ=255)
                                                                                                        gggTGTACGCAAGTTGACGCCCC‑GGGCCCGGTAATACCGGCCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgc                    <  6:357578/136‑1 (MQ=255)
                                                                                                         ggtcaaCTGAAGTTGACGACCC‑GGCCCCGTTAATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgc                    <  7:288584/130‑1 (MQ=255)
                                                                                                         ggTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGAcgc                    <  5:156738/135‑1 (MQ=255)
                                                                                                          ttgtACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc                  <  8:207441/135‑1 (MQ=255)
                                                                                                          gttcaCTGAAGTTGCCGCGCC‑GGGCCCGGCCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc                  <  8:36110/132‑1 (MQ=255)
                                                                                                          gtgtAGTGAAGTTGACGAGCC‑GGGCCCGGTAATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc                  <  7:353134/136‑1 (MQ=255)
                                                                                                          gtgtACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc                  <  7:300780/136‑1 (MQ=255)
                                                                                                          gtgtACTCAAGTTGCCGAGCC‑GGGCCCGGGCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc                  <  8:212245/136‑1 (MQ=255)
                                                                                                            ttACTGAAGTGGCCGAGCC‑GGGCCGGGTAATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccc                <  1:128251/135‑1 (MQ=255)
                                                                                                            gcagtGAAGTTGACGCGCC‑GGGCCCGGTCATACGGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccc                <  7:252483/132‑1 (MQ=255)
                                                                                                             tACTTACGTTGACGAGCC‑GGGCCCGGTCATACCGGGCACGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccc                <  7:101538/135‑1 (MQ=255)
                                                                                                             tACTGACGTGGCCGAGCC‑GGGCCGGGTAATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccc                <  7:11881/135‑1 (MQ=255)
                                                                                                             tACTGAAGTGCACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccct               <  3:322852/136‑1 (MQ=255)
                                                                                                                  acgTTGCCGCGCC‑GGGCCGGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccctgccc           <  7:217237/133‑1 (MQ=255)
                                                                                                                  aaGTTCACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccctgccc           <  8:34068/135‑1 (MQ=255)
                                                                                                                          gcccc‑gggccgggTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAacac  <  1:17901/133‑1 (MQ=255)
                                                                                                                                        |                                                                                                                          
CATCTCCTACCAGGGCACCGGTCTTTGCCTGCGG‑AACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCC‑GGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACAC  >  ACB143/4772412‑4772668

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