Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 4,513,061 T→C pseudogene (312/768 nt) fecE ← Fe(3+) dicitrate ABC transporter ATP‑binding protein FecE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252684,513,0610TC99.7% 1055.3 / NA 317pseudogene (312/768 nt)fecEFe(3+) dicitrate ABC transporter ATP‑binding protein FecE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (142/174);  minor base G (1/0);  total (143/174)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

GCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATA  >  NZ_CP025268/4512924‑4513197
                                                                                                                                         |                                                                                                                                        
gCCAGAAATGCGCGCTGGCGCTGACCGCAGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTACCTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCa                                                                                                                                        <  1:1310689/140‑1 (MQ=255)
   aGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGGTAACCGACGAACGGCAAGATGACGGATCCGGGTCGGGCTCATGGCGACATCAACGCGTGCATTGTCTTCAGCGGAGAGACGGCCCCACAGGGACAGCCAgcg                                                                                                                                     >  2:281091/1‑138 (MQ=255)
     aaaTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGAt                                                                                                                                   <  2:2427918/140‑1 (MQ=255)
      aaTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGAtt                                                                                                                                  <  2:5066016/140‑1 (MQ=255)
           gcgctggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGAc                                                                                                                             <  1:828485/140‑1 (MQ=255)
            cgctggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGAcc                                                                                                                            >  1:2265054/1‑140 (MQ=255)
               tggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATa                                                                                                                         <  2:2604094/140‑1 (MQ=255)
                      gacagacgGAAAGCTCGGTTAACCGCCCAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTACCTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  2:1026054/134‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:2862247/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:746601/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:64187/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:4789049/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  2:3948204/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  2:4298897/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  2:4664097/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  2:4992067/140‑1 (MQ=255)
                       aCCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCa                                                                                                                 >  1:1700468/1‑140 (MQ=255)
                        ccgccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAg                                                                                                                <  1:3232244/140‑1 (MQ=255)
                           ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACCTTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTc                                                                                                             >  2:4459416/1‑140 (MQ=255)
                           ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTc                                                                                                             >  1:4342358/1‑140 (MQ=255)
                                aaGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                                        <  2:1220008/140‑1 (MQ=255)
                                aaGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                                        <  2:1585861/140‑1 (MQ=255)
                                aaGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                                        <  2:3998264/140‑1 (MQ=255)
                                aaGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                                        <  2:1979775/140‑1 (MQ=255)
                                       gTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGTCCATACGATACCAGCTCCTGGACTGTGAt                                                                                                 <  2:3654438/140‑1 (MQ=255)
                                            ccGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                            <  2:2334336/140‑1 (MQ=255)
                                             cgacgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctc                                                                                           >  2:2696590/1‑140 (MQ=255)
                                             cgacgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctc                                                                                           <  1:701873/140‑1 (MQ=255)
                                             cgacgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctc                                                                                           <  2:4801792/140‑1 (MQ=255)
                                               acgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTg                                                                                         <  1:173995/140‑1 (MQ=255)
                                               acgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTg                                                                                         <  1:1484521/140‑1 (MQ=255)
                                                  aaCGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAGATCGTGCATTGTCTTCAGCGGCGCGACGCCCCCATAGTGACAGCCACGGATTACGACCCTACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCg                                                                                      <  2:882003/140‑1 (MQ=255)
                                                  aaCGGCAAGATGATCGGTCCGGGTCTGGTTCATGGCGACATCACCTCGTACATTGTCTCCAGCGGAGAGACCCCCCCAGAGCGACCGCCGAGGATTACGACCATACGAAACCAGCTCCCGGACCGTGATCCCCTCAGGcc                                                                                      >  2:1483192/1‑139 (MQ=255)
                                                    cGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGtt                                                                                    <  1:2510882/140‑1 (MQ=255)
                                                       cAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                                 <  2:341811/140‑1 (MQ=255)
                                                       cAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                                 <  1:2356719/140‑1 (MQ=255)
                                                       cAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                                 <  1:1273771/140‑1 (MQ=255)
                                                                   tCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACTGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                                     <  2:4707964/140‑1 (MQ=255)
                                                                     cGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTTGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGATgcag                                                                   <  2:1899417/140‑1 (MQ=255)
                                                                     cGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgcag                                                                   >  1:2433358/1‑140 (MQ=255)
                                                                     cGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgcag                                                                   >  1:3697527/1‑140 (MQ=255)
                                                                        gTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGcagcag                                                                <  2:4937422/140‑1 (MQ=255)
                                                                            ggTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGAGTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTCAATGGTGCGGAGGCAGCAGCGaa                                                            >  2:784918/1‑140 (MQ=255)
                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  2:300185/85‑1 (MQ=255)
                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGaaa                                                           >  1:1106852/1‑140 (MQ=255)
                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGGGACATCCAGGGATGACGATCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGTCGTGAAATGGTGCGGAGGCGGCAGCGaaa                                                           >  1:1356891/1‑140 (MQ=255)
                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAAACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGaaa                                                           >  2:1568616/1‑140 (MQ=255)
                                                                             gTTCATGGCGACATTAACTCGTGCATTGTATTCAGCGGAGAGACGCCCCCCGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTCGCGTTAAATGGTGCTGAGGCAGCAGCGaaa                                                           >  2:4210106/1‑140 (MQ=255)
                                                                                  tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTg                                                      >  1:1076303/1‑140 (MQ=255)
                                                                                     cGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCgg                                                   <  1:2161724/140‑1 (MQ=255)
                                                                                      gACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCggg                                                  <  1:4507602/140‑1 (MQ=255)
                                                                                          ttAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCaa                                              <  2:81783/140‑1 (MQ=255)
                                                                                          ttAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCaa                                              <  2:3720728/140‑1 (MQ=255)
                                                                                              cTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTgg                                          <  1:3121282/140‑1 (MQ=255)
                                                                                                  tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGcgcg                                      >  1:1568759/1‑140 (MQ=255)
                                                                                                  tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGcgcg                                      <  1:1643217/140‑1 (MQ=255)
                                                                                                  tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGcgcg                                      >  1:1790517/1‑140 (MQ=255)
                                                                                                  tGCATTGTCTGCAGCGGAGGGAGGCGCCCAGAGGGACAGCCAGGGGTTACGACCATACGACACCAGGCCCTGGACTGTGATCCCCTCTGGCGGGAAAGGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGcgcg                                      >  2:654184/1‑140 (MQ=255)
                                                                                                   gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGCCAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCGTGCGGGCCAACTGGCGCGa                                     >  2:3167773/1‑140 (MQ=255)
                                                                                                    ctttatCTTCAGCGAGGAGACGCCCCCACGGCGACAGCGAAGGATTACGACCATACGAACCCAGCCCTTCGACTGTGAACCCCCCTGGCGTGAAATGGTGCGGAGGCAGCCGCGAAAGCCTGCGGGACAACTGGCGCGAt                                    <  2:4054725/135‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  2:4988081/140‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  2:3298229/140‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  1:2083329/140‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  2:1712979/140‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  2:1880757/140‑1 (MQ=255)
                                                                                                      ttGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATga                                  <  2:2577183/140‑1 (MQ=255)
                                                                                                       tGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATgag                                 >  1:4895996/1‑140 (MQ=255)
                                                                                                         tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                               <  1:3723401/140‑1 (MQ=255)
                                                                                                         tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                               >  2:2516065/1‑140 (MQ=255)
                                                                                                         tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                               <  2:3121091/140‑1 (MQ=255)
                                                                                                         tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                               >  1:1194707/1‑140 (MQ=255)
                                                                                                          cTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                              <  2:2179605/140‑1 (MQ=255)
                                                                                                              aGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCatat                          >  2:2860780/1‑140 (MQ=255)
                                                                                                              aGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCatat                          >  1:2757024/1‑140 (MQ=255)
                                                                                                              aGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCatat                          >  2:4344682/1‑140 (MQ=255)
                                                                                                              aGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGCTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACGGGCGCCAGGCGAGCacgt                          >  2:1161329/1‑137 (MQ=255)
                                                                                                                  gagagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTAt                      <  2:2265812/140‑1 (MQ=255)
                                                                                                                   atagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCATGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCAAATTTacg                     >  2:2478869/3‑138 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAGCCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAGCTGGCGCGATGAGAGCATATTTATg                     >  2:3929074/1‑140 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATg                     >  1:4471097/1‑140 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATg                     >  2:3594721/1‑140 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATg                     >  1:1953959/1‑140 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATg                     >  1:2026387/1‑140 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATg                     >  1:252168/1‑140 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATg                     >  1:4155724/1‑140 (MQ=255)
                                                                                                                        cGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGAtt                >  2:3836895/1‑140 (MQ=255)
                                                                                                                          cccccAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGTGATTAt              >  2:4948720/1‑140 (MQ=255)
                                                                                                                          cccccAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGCGCTGAGGCAGCAGCGAACGCCTGCGGGCCAACTGGCGCGATGAGCGCATATTTATGGGCTTAt              >  2:537468/1‑140 (MQ=255)
                                                                                                                              cAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGcc          <  1:3343009/140‑1 (MQ=255)
                                                                                                                               agagTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCg         <  1:2917038/140‑1 (MQ=255)
                                                                                                                                 agTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCgag       <  2:3016410/140‑1 (MQ=255)
                                                                                                                                     aCAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAAt   >  1:4100150/1‑140 (MQ=255)
                                                                                                                                     aCAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAAt   >  2:3246349/1‑140 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATCCGAAACCAGCTCCTGGACTGTGATCGCCTCTGGCGTTAAATGGTGCGGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  1:4604393/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  2:4471097/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  1:4856042/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  1:3516561/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  1:4930735/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCAGATTTATGGGATTATCGCCGAGAAATa  <  1:3969122/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCGGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  1:510851/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATA  >  NZ_CP025268/4512924‑4513197

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