Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,415,941 G→A 6.0% P80P (CCG→CCA CJ019_04188 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,415,9410GA 6.0% 205.9 / 3.2 66P80P (CCG→CCACJ019_04188hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/39);  new base A (2/2);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GCAACCCACACCGAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGCC  >  ACB122/4415811‑4416075
                                                                                                                                  |                                                                                                                                      
gCAACCCACACCGAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGc                                                                                                                                    <  6:260949/135‑1 (MQ=255)
   aCCCTCACCGAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCag                                                                                                                                 <  6:142544/135‑1 (MQ=255)
    cccTCAGCGAGTGTTCGCCCGCTACCAGGGCGAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCagag                                                                                                                               <  8:369258/136‑1 (MQ=255)
         acCGAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGCGCACACCAAGCGCTTCAACCCGTCGGAAGACTGGGCCGTTGGCGGGCCGAACATGCAGAAGCACAACGACAGCCTGCACAGCCCGCAGCCAGACAGGGa                                                                                                                          >  4:114610/1‑136 (MQ=255)
         acCGAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGa                                                                                                                          >  3:271002/1‑136 (MQ=255)
         acCGAGTGTTCCCCCGCTACCAGGGCCACTCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATCCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCACAGTgcg                                                                                                                           >  6:57916/1‑133 (MQ=255)
          ccGAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGa                                                                                                                          <  1:191172/135‑1 (MQ=255)
          ccGAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGa                                                                                                                          <  6:280884/135‑1 (MQ=255)
            gAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGacta                                                                                                                       >  7:126395/1‑136 (MQ=255)
            gAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGact                                                                                                                        >  6:249606/1‑135 (MQ=255)
             aGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGactac                                                                                                                      <  8:126395/136‑1 (MQ=255)
               tgtTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGactact                                                                                                                     <  1:56498/135‑1 (MQ=255)
               tgtTCGCCCGCGACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGactact                                                                                                                     <  6:201107/135‑1 (MQ=255)
                gtTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGCCTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTg                                                                                                                    >  2:164519/1‑135 (MQ=255)
                  tCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTgg                                                                                                                   <  1:298328/134‑1 (MQ=255)
                  tCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGc                                                                                                                 >  3:331686/1‑136 (MQ=255)
                    gCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGc                                                                                                                 <  5:231550/134‑1 (MQ=255)
                       cGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCt                                                                                                            <  7:356548/136‑1 (MQ=255)
                        gCTCCCAGGGCCATGCCATCGAGCAGACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACTACGTCAGCCTGCACTGGCCGCTGCCAGAGTGGGACTACTGGGCAGCCt                                                                                                            <  1:282905/135‑1 (MQ=255)
                            ccAGGGCCAGGCCATCGAGCACACCACGCGCGTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGa                                                                                                         >  2:379216/1‑134 (MQ=255)
                                ggccaggccaTCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAAcc                                                                                                    <  6:18082/135‑1 (MQ=255)
                                       ccaTCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCa                                                                                                                                  >  5:216331/1‑98 (MQ=255)
                                       ccaTCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCa                                                                                                                                  <  6:216331/98‑1 (MQ=255)
                                       ccaTCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACgg                                                                                            >  8:174815/1‑136 (MQ=255)
                                       ccaTCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACg                                                                                             >  1:129168/1‑135 (MQ=255)
                                       ccaTCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTCGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAAcc                                                                                             >  2:282905/1‑134 (MQ=255)
                                        caTCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACgg                                                                                            <  6:105597/135‑1 (MQ=255)
                                            gAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAgg                                                                                       <  7:221083/136‑1 (MQ=255)
                                              gCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGAc                                                                                     >  3:168022/1‑136 (MQ=255)
                                              gCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGAc                                                                                     >  4:89221/1‑136 (MQ=255)
                                                  acCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCg                                                                                 >  4:1519/1‑136 (MQ=255)
                                                   ccAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGc                                                                                <  3:201436/136‑1 (MQ=255)
                                                   ccAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGc                                                                                <  7:174815/136‑1 (MQ=255)
                                                       gcgcTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAACCTGCACTGCCCGCAGCCAGAGTCGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCa                                                                             >  6:1359/1‑135 (MQ=255)
                                                        cgcTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAgg                                                                           <  3:114610/136‑1 (MQ=255)
                                                          cTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGc                                                                                                                                       <  3:266750/74‑1 (MQ=255)
                                                          cTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGc                                                                                                                                       >  4:266750/1‑74 (MQ=255)
                                                             cAACCCGTGGGAAGACTGGGCGGGTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCt                                                                      <  8:58638/136‑1 (MQ=255)
                                                               aCCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGa                                                                    >  3:41009/1‑136 (MQ=255)
                                                               aCCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCACAGCCAGAGTGGGACTACTGGGc                                                                                                                 <  2:216255/91‑1 (MQ=255)
                                                               aCCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCACAGCCAGAGTGGGACTACTGGGc                                                                                                                 <  2:209782/91‑1 (MQ=255)
                                                               aCCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCACAGCCAGAGTGGGACTACTGGGc                                                                                                                 >  1:209782/1‑91 (MQ=255)
                                                               aCCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCACAGCCAGAGTGGGACTACTGGGc                                                                                                                 >  1:216255/1‑91 (MQ=255)
                                                                   gTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCt                                                                 <  5:143044/135‑1 (MQ=255)
                                                                        aaGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACg                                                                                             >  7:213797/1‑102 (MQ=255)
                                                                        aaGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACg                                                                                             <  8:213797/102‑1 (MQ=255)
                                                                                 cGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCggcggc                                                  <  4:331686/136‑1 (MQ=255)
                                                                                  ggTTGGCGGTCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGgcggcg                                                 <  8:325982/136‑1 (MQ=255)
                                                                                  ggTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGgcggcg                                                 <  8:242947/136‑1 (MQ=255)
                                                                                   gatGGCGGGCCGATCATGCAGAAGCACACCGTCAGCCTGCACTGGCCGCAGCCAGAGGGGGACTACTGGGCAGCCTGTATAAGCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGgcggcg                                                 <  5:57916/133‑1 (MQ=255)
                                                                                   gTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGgcggcg                                                 <  1:370596/135‑1 (MQ=255)
                                                                                   gTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGgcggcg                                                 <  5:227819/135‑1 (MQ=255)
                                                                                   gTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGgcggcg                                                 <  2:129168/135‑1 (MQ=255)
                                                                                   gTTGGCGGGCCGATCATGCAGAAGCACAACGGCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGgcggcg                                                 <  6:269026/135‑1 (MQ=255)
                                                                                             cGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCAt                                      <  8:7059/136‑1 (MQ=255)
                                                                                                tCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAg                                    <  5:15882/135‑1 (MQ=255)
                                                                                                  aTGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTc                                  <  7:41657/135‑1 (MQ=255)
                                                                                                       gAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCccac                            <  3:380823/136‑1 (MQ=255)
                                                                                                       gAAGCACAACGGCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCccac                            <  7:211923/136‑1 (MQ=255)
                                                                                                        aaGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCccac                            <  1:379216/135‑1 (MQ=255)
                                                                                                             caACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTTGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCt                       >  6:195402/1‑135 (MQ=255)
                                                                                                                cGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATACCCGATAACGGCAAGGACGTCGCCCAGGGCGCTGAGCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCgg                    <  1:187373/135‑1 (MQ=255)
                                                                                                                cGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCgg                    >  5:388017/1‑135 (MQ=255)
                                                                                                                             tGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTccctg      >  3:288132/1‑133 (MQ=255)
                                                                                                                                 cGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGc                                                                                <  3:271724/58‑1 (MQ=255)
                                                                                                                                 cGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGc                                                                                >  4:271724/1‑58 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCccac                            <  7:148863/109‑1 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCccac                            >  8:148863/1‑109 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCccac                            <  2:229607/109‑1 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCccac                            >  1:229607/1‑109 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATCCCGTCCAGGTggc  >  1:178458/1‑133 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGcc  >  3:93298/1‑135 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGcc  >  4:301911/1‑135 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGcc  >  2:347331/1‑135 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGcc  >  8:219463/1‑135 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGcc  >  1:378662/1‑135 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGcc  >  4:92778/1‑135 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGc   >  2:305241/1‑134 (MQ=255)
                                                                                                                                  gCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGc   >  6:357217/1‑134 (MQ=255)
                                                                                                                                  |                                                                                                                                      
GCAACCCACACCGAGTGTTCGCCCGCTACCAGGGCCAGGCCATCGAGCACACCAAGCGCTTCAACCCGTGGGAAGACTGGGCGGTTGGCGGGCCGATCATGCAGAAGCACAACGTCAGCCTGCACTGCCCGCAGCCAGAGTGGGACTACTGGGCAGCCTGGATAACCGATAACGGCAAGGACGTCGCCCAGGGCGCTGATCTGCCGTTGCCGGCGGCGTGCCGGGCCATAGTCGCCCACCAGCTCGGCGATACCGTCCAGGTGCC  >  ACB122/4415811‑4416075

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