Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 376,966 A→C 9.2% K781T (AAA→ACA)  srfAA → surfactin non‑ribosomal peptide synthetase SrfAA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP034943376,9660AC 9.2% 809.6 / 14.6 357K781T (AAA→ACA) srfAAsurfactin non‑ribosomal peptide synthetase SrfAA
Reads supporting (aligned to +/- strand):  ref base A (145/163);  new base C (46/2);  total (191/166)
Fisher's exact test for biased strand distribution p-value = 8.50e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

GTTGAAGCTATTTAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGACGGCTGAG  >  NZ_CP034943/376829‑377102
                                                                                                                                         |                                                                                                                                        
   gAAGCTATTTAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAgcgc                                                                                                                                     <  1:3077179/140‑1 (MQ=255)
   gAAGCTATTTAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAgcgc                                                                                                                                     <  2:1024901/140‑1 (MQ=255)
   gAAGCTATTTAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAgcgc                                                                                                                                     <  1:3348523/140‑1 (MQ=255)
     aGCTATTTAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTg                                                                                                                                   <  1:314695/140‑1 (MQ=255)
     aGCTATTTAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTg                                                                                                                                   <  2:1400563/140‑1 (MQ=255)
         aTTTAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCAGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATAAGCGTATGTGCCGATCGGGAAAGCGCTGGGGa                                                                                                                               >  2:328545/1‑140 (MQ=255)
          tttCAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGaa                                                                                                                              <  1:2195869/140‑1 (MQ=255)
           ttAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGact                                                                                                                             >  1:3538124/1‑138 (MQ=255)
           ttAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAAt                                                                                                                             >  1:1979840/1‑140 (MQ=255)
           ttAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAAt                                                                                                                             >  2:1349698/1‑140 (MQ=255)
            tAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATa                                                                                                                            <  2:695786/140‑1 (MQ=255)
             aaaGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATaa                                                                                                                           <  1:2122446/140‑1 (MQ=255)
                    cAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAACGCGCTGGGGAATAACCGTTta                                                                                                                    >  2:4250549/1‑140 (MQ=255)
                    cAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTta                                                                                                                    >  2:2937817/1‑140 (MQ=255)
                     aGGCACATAGCAGCGTTTGACCAATTTGTACGGGCCGACTGAAACGTGTGTTGAAGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTAAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTtat                                                                                                                   <  1:434786/140‑1 (MQ=255)
                     aGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTtat                                                                                                                   <  1:3762090/140‑1 (MQ=255)
                             cgccgcGTTTGACCAATGTGTACGGGCCGACTGACACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGAATAACCGTTTATATATTTTg                                                                                                           >  2:3344249/1‑140 (MQ=255)
                             cgccgcGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGACGAACCGTTTATATATTTTg                                                                                                           >  1:3718366/1‑140 (MQ=255)
                             cgccgcGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTg                                                                                                           >  2:3322522/1‑140 (MQ=255)
                             cgccgcGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTg                                                                                                           >  2:434189/1‑140 (MQ=255)
                              gccgcGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTgg                                                                                                          <  2:3006284/140‑1 (MQ=255)
                               ccgcGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCCTCAGTTTATCCCAGAGAGTGCCTTTCAATCAGCGTATGTGCCGATCGTGACAGCGCTGGGGAATCACCGTTTATATATTTGGGa                                                                                                         >  1:2132969/1‑140 (MQ=255)
                               ccgcGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGa                                                                                                         >  1:3697062/1‑140 (MQ=255)
                                  cgTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGACTAACCGTTTATATATTTTGGATCa                                                                                                      >  1:1020245/1‑140 (MQ=255)
                                       gACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGACTAACCGTTTATATATTTTGGATCAAagag                                                                                                 >  1:477377/1‑140 (MQ=255)
                                       gACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCACATAGAGTGCCGTTCAATCATCGTATGTGCAGATCTGGAACGCGCTGGGGAATAACCGTTTATATATTTTGTATCAAagag                                                                                                 >  1:505594/1‑140 (MQ=255)
                                        aCCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCGGGGGAATAACCGTTTATATATTTTGGATCAAAGagg                                                                                                <  2:177092/140‑1 (MQ=255)
                                            aTGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAACGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGaggcagg                                                                                            >  2:526451/1‑140 (MQ=255)
                                               tgtACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTg                                                                                         <  1:1600610/140‑1 (MQ=255)
                                                    gggCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGcc                                                                                    <  1:550327/140‑1 (MQ=255)
                                                    gggCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGcc                                                                                    <  2:1550297/140‑1 (MQ=255)
                                                    gggCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGcc                                                                                    <  1:1567331/140‑1 (MQ=255)
                                                    gggCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTGATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGcc                                                                                    <  1:808685/140‑1 (MQ=255)
                                                      gCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCgg                                                                                  >  1:2544606/1‑140 (MQ=255)
                                                      gCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCgg                                                                                  >  1:3432238/1‑140 (MQ=255)
                                                         gACTGACCCGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATCACCGGTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAg                                                                               >  2:863883/1‑140 (MQ=255)
                                                         gACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCCGCGGATGTGCCGATCGGGACAGCGCTGGGGAATAACCGTTTATATCTTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAg                                                                               >  2:3632281/1‑140 (MQ=255)
                                                         gACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGAATAACCGTGTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAg                                                                               >  2:1703650/1‑140 (MQ=255)
                                                         gACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAg                                                                               >  2:2408937/1‑140 (MQ=255)
                                                                    gtgtTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTggcgggc                                                                    <  2:345898/140‑1 (MQ=255)
                                                                    gtgtTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTggcgggc                                                                    <  1:2925094/140‑1 (MQ=255)
                                                                    gtgtTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTggcgggc                                                                    <  2:249312/140‑1 (MQ=255)
                                                                     tgtTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTgtcgggcg                                                                   >  1:1650154/1‑140 (MQ=255)
                                                                       ttGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAg                                                                 <  1:2658872/140‑1 (MQ=255)
                                                                       ttGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAg                                                                 <  1:3264408/140‑1 (MQ=255)
                                                                          aCGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCtt                                                              >  1:2166814/1‑140 (MQ=255)
                                                                                 cGTTCATGAGCTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCg                                                       <  1:2812614/140‑1 (MQ=255)
                                                                                 cGTTCATCAGGTTATCCCAGAGCGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTGTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCg                                                       <  2:1517891/140‑1 (MQ=255)
                                                                                 cGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTggcgggc                                                                    >  2:289760/1‑127 (MQ=255)
                                                                                 cGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTggcgggc                                                                    <  1:289760/127‑1 (MQ=255)
                                                                                        ccggTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  1:1366003/138‑1 (MQ=255)
                                                                                        caGTTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGCTCGCGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCACAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  2:266310/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGCGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGTGACTAACCGTTTATATATGTTGGATCAAAGAGGCAGTCTGCAGCCGGAAGGTGTGCCGGGCGAGCTTTATATCGCAGGCGa                                                <  1:2381480/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  2:3489726/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  1:1938729/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  1:2408937/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  1:3642983/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  2:1143662/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  2:4214012/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  1:1141695/140‑1 (MQ=255)
                                                                                        caGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  1:3014565/140‑1 (MQ=255)
                                                                                        caGGTGATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                <  1:4257/140‑1 (MQ=255)
                                                                                         aGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGAc                                               <  1:3969141/140‑1 (MQ=255)
                                                                                              aTCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAg                                                                 >  2:4069091/1‑117 (MQ=255)
                                                                                              aTCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAg                                                                 <  1:4069091/117‑1 (MQ=255)
                                                                                                cccAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGACCGCGCTGGGGACTAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTg                                        >  2:1196135/1‑140 (MQ=255)
                                                                                                    gagagTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCg                                    <  2:2132969/140‑1 (MQ=255)
                                                                                                    gagagTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCg                                    <  2:1943993/140‑1 (MQ=255)
                                                                                                     agagTGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGa                                   >  2:3215241/1‑140 (MQ=255)
                                                                                                     agagTGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGGGTTGGCCGa                                   >  2:2212279/1‑140 (MQ=255)
                                                                                                     agagTGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGAATAACCGGTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGa                                   >  2:2772280/1‑140 (MQ=255)
                                                                                                     agagTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGa                                   >  1:1836128/1‑140 (MQ=255)
                                                                                                       agTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAgg                                 <  1:4116585/140‑1 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  2:616286/1‑140 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGACAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  2:3692574/1‑140 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGAACGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  2:4147111/1‑140 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  2:3666846/1‑140 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  1:1208759/1‑140 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  2:3761504/1‑140 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  1:3863157/1‑140 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  2:2992285/1‑140 (MQ=255)
                                                                                                         tGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGGTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCt                               >  2:1654111/1‑140 (MQ=255)
                                                                                                          gCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTggcgggc                                                                    <  1:1677446/102‑1 (MQ=255)
                                                                                                           ccGTTCAATCAGCGTCGGTGCCGATCGGGACCGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTGATATCGCAGGCGACGGAGTTGGCCGAGGCTAt                             >  2:2062170/1‑140 (MQ=255)
                                                                                                              ttCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGa                                                                                 >  1:1864317/1‑85 (MQ=255)
                                                                                                              ttCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGa                                                                                 <  2:1864317/85‑1 (MQ=255)
                                                                                                                 aaTCCGCGTATGTGGCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGTGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGGGTTGGCCGAGGCTATTTACAt                       >  2:1272625/1‑140 (MQ=255)
                                                                                                                 aaTCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGa                                                >  2:1059176/1‑115 (MQ=255)
                                                                                                                   tCAGCGTATGTGCCGATCGGGACAGCGCTGGGGACTAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACAttt                     >  1:2993227/1‑140 (MQ=255)
                                                                                                                     aGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGc                   <  1:4228079/140‑1 (MQ=255)
                                                                                                                     aGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGc                   <  2:3984103/140‑1 (MQ=255)
                                                                                                                       cGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCt                 <  2:227618/140‑1 (MQ=255)
                                                                                                                            gtgCCGATCGGGAAAGCGCTGGGGACTACCCGTTTATATATTTTGGCTCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAAtt            <  2:2821776/140‑1 (MQ=255)
                                                                                                                              gCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGa          <  1:4178386/140‑1 (MQ=255)
                                                                                                                              gCCGATCGGCAAAGCGCGGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGa          <  2:3397129/140‑1 (MQ=255)
                                                                                                                                 gATCGGGAACGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGACgg       >  1:4322954/1‑140 (MQ=255)
                                                                                                                                 gATCGGGAACGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGACgg       >  2:792020/1‑140 (MQ=255)
                                                                                                                                 gATCGGGAAAGCGCTGGGGACTAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGACgg       >  1:2301162/1‑140 (MQ=255)
                                                                                                                                 gATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGACgg       >  1:658305/1‑140 (MQ=255)
                                                                                                                                      ggAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGACGGCTgag  <  1:1748357/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GTTGAAGCTATTTAAAGAAGCAGGCACAGCGCCGCGTTTGACCAATGTGTACGGGCCGACTGAAACGTGTGTTGACGCGTCCGTTCATCAGGTTATCCCAGAGAGTGCCGTTCAATCAGCGTATGTGCCGATCGGGAAAGCGCTGGGGAATAACCGTTTATATATTTTGGATCAAAGAGGCAGGCTGCAGCCGGAAGGTGTGGCGGGCGAGCTTTATATCGCAGGCGACGGTGTTGGCCGAGGCTATTTACATTTGCCTGAATTGACGGCTGAG  >  NZ_CP034943/376829‑377102

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: