Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 247,291 A→G *250* (TAA→TAG yafL → C40 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP025268247,2910AG99.6% 907.7 / ‑5.9 276*250* (TAA→TAGyafLC40 family peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (120/155);  minor base T (0/1);  total (120/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 276 total aligned reads displayed.

GTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGCC  >  NZ_CP025268/247154‑247423
                                                                                                                                         |                                                                                                                                      
gtgtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGa                                                                                                                                      <  1:1151592/140‑1 (MQ=255)
gtgtATTTGGGCGATGGGAAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGa                                                                                                                                      <  1:4473630/140‑1 (MQ=255)
 tgtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGAc                                                                                                                                     >  1:4593151/1‑140 (MQ=255)
 tgtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGAc                                                                                                                                     >  2:819093/1‑140 (MQ=255)
 tgtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGAc                                                                                                                                     >  1:5112174/1‑140 (MQ=255)
 tgtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGAc                                                                                                                                     >  1:2178872/1‑140 (MQ=255)
  gtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACg                                                                                                                                   >  1:5014702/1‑140 (MQ=255)
  gtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACg                                                                                                                                   >  1:3514903/1‑140 (MQ=255)
  gtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACg                                                                                                                                   >  1:2034624/1‑140 (MQ=255)
   tATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGCTTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACgg                                                                                                                                  >  2:1807007/1‑139 (MQ=255)
   tATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACgg                                                                                                                                  >  1:371428/1‑139 (MQ=255)
       tGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACggataa                                                                                                                              <  1:2327594/140‑6 (MQ=255)
        gggCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATCAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAg                                                                                                                             <  1:2941525/140‑1 (MQ=255)
        gggCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGAGTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGATACGTTTTTGTAGGACTGATAAg                                                                                                                             >  2:4540261/1‑140 (MQ=255)
         ggCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAgg                                                                                                                            >  2:4662523/1‑140 (MQ=255)
         ggCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAgg                                                                                                                            >  2:3674880/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCTACGTACCGGCGAAACCATTCGGATAAGCCGGTTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  2:2664305/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  2:3846066/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  1:2601441/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  1:264672/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  2:1164558/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  2:1211557/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  2:1743680/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        <  2:2561219/140‑1 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  2:603194/1‑140 (MQ=255)
             aTGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGtt                                                                                                                        >  1:374912/1‑140 (MQ=255)
                ggCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTAc                                                                                                                     >  2:3345783/1‑140 (MQ=255)
                ggCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTAc                                                                                                                     >  1:3676533/1‑140 (MQ=255)
                    aTTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTAcgccg                                                                                                                 <  1:1807670/140‑1 (MQ=255)
                         tCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATcc                                                                                                            <  1:5313653/140‑1 (MQ=255)
                           gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCgg                                                                                                          <  1:952949/140‑1 (MQ=255)
                           gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCgg                                                                                                          <  2:72583/140‑1 (MQ=255)
                           gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCgg                                                                                                          <  1:441816/140‑1 (MQ=255)
                           gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCgg                                                                                                          <  1:2918864/140‑1 (MQ=255)
                           gAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCgg                                                                                                          >  1:1291221/1‑140 (MQ=255)
                              tCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAg                                                                                                       <  1:112297/140‑1 (MQ=255)
                                  cACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGt                                                                                                   <  2:2962373/140‑1 (MQ=255)
                                   aCGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTa                                                                                                  <  1:1376327/140‑1 (MQ=255)
                                   aCGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTa                                                                                                  <  1:2869977/140‑1 (MQ=255)
                                   aCGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTa                                                                                                  <  1:2069424/140‑1 (MQ=255)
                                     gTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTCTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACg                                                                                                <  2:4899265/140‑1 (MQ=255)
                                       aCCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCa                                                                                              <  1:807444/140‑1 (MQ=255)
                                          ggCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGt                                                                                           <  1:3881570/140‑1 (MQ=255)
                                           gCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTGTGGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTg                                                                                          >  2:291116/1‑140 (MQ=38)
                                            cGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGc                                                                                         >  1:1440116/1‑140 (MQ=255)
                                              aaaCCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCt                                                                                       <  1:1870419/140‑1 (MQ=255)
                                              aaaCCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCt                                                                                       <  1:685016/140‑1 (MQ=255)
                                                         aTAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCt                                                                            >  1:812621/1‑140 (MQ=255)
                                                           aaGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTgg                                                                          <  1:4532173/140‑1 (MQ=255)
                                                            aGCCGATTAGCCTAACCTTTCTGGCAGGACCTTTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACTGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGTAGGTGCCTGATGCGACGCTGgc                                                                         <  2:4419436/140‑1 (MQ=38)
                                                            aGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         <  1:3621014/140‑1 (MQ=255)
                                                            aGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         <  1:3407056/140‑1 (MQ=255)
                                                            aGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         <  1:4295168/140‑1 (MQ=255)
                                                            aGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         <  1:4229024/140‑1 (MQ=255)
                                                            aGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         <  1:4689675/140‑1 (MQ=255)
                                                            aGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         <  2:1556602/140‑1 (MQ=255)
                                                            aGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGGTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                                         <  2:3910561/140‑1 (MQ=255)
                                                             gCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgcg                                                                        <  1:3054031/140‑1 (MQ=255)
                                                              ccGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgcgc                                                                       >  1:4148567/1‑140 (MQ=255)
                                                                gATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGt                                                                     <  2:1399858/140‑1 (MQ=255)
                                                                gATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGt                                                                     <  2:5154602/140‑1 (MQ=255)
                                                                 aTTAGCCGAACGTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAATAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTc                                                                    >  2:3808465/1‑140 (MQ=37)
                                                                 aTTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTc                                                                    >  2:1045055/1‑140 (MQ=255)
                                                                  ttAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCt                                                                   >  1:4798029/1‑140 (MQ=255)
                                                                    aGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCtga                                                                 <  1:2499753/140‑3 (MQ=255)
                                                                    aGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTa                                                                 <  1:1012656/140‑1 (MQ=255)
                                                                      ccGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATc                                                               <  2:4616137/140‑1 (MQ=255)
                                                                      ccGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATc                                                               <  2:3137321/140‑1 (MQ=255)
                                                                        gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCTGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                             <  2:2286239/140‑1 (MQ=255)
                                                                        gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                             <  2:2633157/140‑1 (MQ=255)
                                                                        gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                             <  1:3279948/140‑1 (MQ=255)
                                                                              ttCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTAc                                                       >  1:5228146/1‑140 (MQ=255)
                                                                               tCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACg                                                      >  1:2421554/1‑140 (MQ=255)
                                                                                cTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGa                                                     <  2:426665/140‑1 (MQ=255)
                                                                                 tGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAg                                                    <  2:2200255/140‑1 (MQ=255)
                                                                                 tGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAg                                                    <  2:4323502/140‑1 (MQ=255)
                                                                                    cAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGccc                                                 <  1:2213771/140‑1 (MQ=255)
                                                                                      ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCa                                                                                              <  2:4944117/93‑1 (MQ=37)
                                                                                      ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCcgc                                               <  2:2996779/140‑1 (MQ=255)
                                                                                      ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCcgc                                               >  1:3060054/1‑140 (MQ=255)
                                                                                        aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGa                                             >  2:2802498/1‑140 (MQ=255)
                                                                                                            aagATTTTGGGGGACGGGGCGATTTTGTATGGGGAATAAGGCGTTGACGCAGCACCCGACAGTTGTACGCAGGTGCCTGATGCGGCGCAGTCGAGTCTAATCATGCCTACGAGCCCGCGCATATTTGCCAGCGGCtttgc                         <  2:580699/138‑3 (MQ=255)
                                                                                                            aGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                         <  2:3452524/140‑1 (MQ=255)
                                                                                                            aGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                         <  2:90991/140‑1 (MQ=255)
                                                                                                            aGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                         <  1:4392286/140‑1 (MQ=255)
                                                                                                            aGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc                         <  1:4210353/140‑1 (MQ=255)
                                                                                                             ggATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTccc                        <  1:1953097/140‑1 (MQ=255)
                                                                                                               cttttGACGGAAGAGACGATTTTGTAGGACGGAGACGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGa                      <  2:2029703/139‑1 (MQ=25)
                                                                                                                ttttGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGat                     <  2:4806485/140‑1 (MQ=255)
                                                                                                                ttttGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGat                     <  2:3983403/140‑1 (MQ=255)
                                                                                                                ttttGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGat                     <  2:3541283/140‑1 (MQ=255)
                                                                                                                  ttGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGatat                   <  2:1487219/140‑1 (MQ=255)
                                                                                                                   tGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTAcg                                                                                                                    <  2:5273204/42‑1 (MQ=18)
                                                                                                                       ggAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAaac              <  2:3059869/140‑1 (MQ=255)
                                                                                                                        gAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAaaca             <  2:328232/140‑1 (MQ=255)
                                                                                                                        gAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAaaca             >  1:488858/1‑140 (MQ=255)
                                                                                                                         aaGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAaacaa            >  2:4559006/1‑140 (MQ=255)
                                                                                                                           gacaCGATTTTGTAGGACGGATAAGGCGTTTACGCCGCCTCCGGCAGTTCTCCGAAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAccc          <  2:1134207/137‑1 (MQ=21)
                                                                                                                                   tttGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGcc  <  1:122455/140‑1 (MQ=255)
                                                                                                                                   tttGTAGGACGGATAAGCCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCGTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGcc  <  1:4059523/140‑1 (MQ=25)
                                                                                                                                         |                                                                                                                                      
GTGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGATATAAAACAACCTCATTGCC  >  NZ_CP025268/247154‑247423

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