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.4% 1021.6 / ‑5.5 312pseudogene (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 (118/191);  minor base G (2/0);  total (120/192)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 313 total aligned reads displayed.

GCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATA  >  NZ_CP025268/4512924‑4513197
                                                                                                                                         |                                                                                                                                        
cccAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCACAGTGACAGCCa                                                                                                                                        <  2:4458456/139‑1 (MQ=255)
   aGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAggg                                                                                                                                     >  1:3524155/1‑140 (MQ=255)
     aaaTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGAt                                                                                                                                   >  1:1320620/1‑140 (MQ=255)
     aaaTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGAt                                                                                                                                   <  1:4639588/140‑1 (MQ=255)
       aTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTa                                                                                                                                 <  1:1661285/140‑1 (MQ=255)
          cgcgcTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGa                                                                                                                              >  2:2345479/1‑140 (MQ=255)
           gcgctggcgctgCCCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGAc                                                                                                                             <  2:1694034/140‑1 (MQ=255)
           gcgctggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGGTCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGAc                                                                                                                             <  2:1671884/140‑1 (MQ=255)
           gcgctggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGAc                                                                                                                             <  1:3131079/140‑1 (MQ=255)
           gcgctggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGAc                                                                                                                             <  2:4070503/140‑1 (MQ=255)
            cgctggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGAcc                                                                                                                            >  1:300877/1‑140 (MQ=255)
                    ctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGaaa                                                                                                                    <  1:3676714/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:75868/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:4886217/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:888198/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:431479/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  1:1827134/140‑1 (MQ=255)
                      gACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAcc                                                                                                                  <  2:896630/140‑1 (MQ=255)
                      gACCGCAAGAAAGCTCCGTTAACCGACGAACGGCAAGATGATTGCTCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTAACAGCCAGCGATTACGACCATACGAAccc                                                                                                                  <  2:967919/140‑1 (MQ=255)
                          gccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCt                                                                                                              <  1:3288146/140‑1 (MQ=255)
                           ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTc                                                                                                             <  1:214631/140‑1 (MQ=255)
                           ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTc                                                                                                             <  2:1778138/140‑1 (MQ=255)
                            cgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGCTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGGGACAGCCAGGGATTACGACCATACGAAACCAGCTcc                                                                                                            >  2:150387/1‑140 (MQ=255)
                             ggAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCt                                                                                                           <  1:2410524/140‑1 (MQ=255)
                              gAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTg                                                                                                          <  1:380111/140‑1 (MQ=255)
                                aaGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCAGGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                                        <  1:4214350/140‑1 (MQ=255)
                                  gCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACt                                                                                                      <  2:4518960/140‑1 (MQ=255)
                                   cTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACtg                                                                                                     >  1:3748499/1‑140 (MQ=255)
                                   cTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGAGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGTCtg                                                                                                     >  1:3728531/1‑140 (MQ=255)
                                     cGGTTAACCGACGAACGGCAAAATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGTCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACtgtg                                                                                                   <  2:3581229/140‑1 (MQ=255)
                                      ggTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGa                                                                                                  >  1:4387081/1‑140 (MQ=255)
                                      ggTTAACCGACGAACGGCAAGATGATTGATCCGGGGCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGa                                                                                                  >  1:2260736/1‑140 (MQ=255)
                                       gTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTACCTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGAt                                                                                                 <  1:3250946/140‑1 (MQ=255)
                                       gTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGAt                                                                                                 <  2:1822682/140‑1 (MQ=255)
                                       gTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGAt                                                                                                 <  1:1385278/140‑1 (MQ=255)
                                       gTTAACCCACCAACGGCAAGATGATTGATCCGGGTCTGGTTCCTGGCGACATTAACCCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGAGTAGGACCATACGCAACCAGCTCCTGGACTGTGAt                                                                                                 <  2:329144/140‑1 (MQ=255)
                                             cgacgaACGGCAAGATGCTTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctc                                                                                           <  2:4901575/140‑1 (MQ=255)
                                             cgacgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctc                                                                                           <  1:4000688/140‑1 (MQ=255)
                                             cgacgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctc                                                                                           <  2:346094/140‑1 (MQ=255)
                                              gacgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctct                                                                                          <  1:5012130/140‑1 (MQ=255)
                                               acgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTg                                                                                         <  2:3914667/140‑1 (MQ=255)
                                                cgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTgg                                                                                        >  2:4211253/1‑140 (MQ=255)
                                                  aaCGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCg                                                                                      <  1:3850659/140‑1 (MQ=255)
                                                  aaCGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGt                                                                                     >  2:2391269/1‑140 (MQ=255)
                                                       cAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                                 <  1:4591360/140‑1 (MQ=255)
                                                               ttgatCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTg                                                                         >  2:2456785/1‑140 (MQ=255)
                                                                     cGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgcag                                                                   <  2:3241885/140‑1 (MQ=255)
                                                                     cGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTAATCCCCTCTGGCGTTAAATGGTGCTGAGgcag                                                                   <  1:2299985/140‑1 (MQ=255)
                                                                        ggcTGGTTCAGGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGGCGCCCCCAAAGTGACAGCCAGGGATGGCGACGGTACGAAACCAGCTCCTGGACTGTGAGCCCCTCGGGCGTGAAATGGTGCGGAGGcagcag                                                                >  2:1044563/3‑140 (MQ=255)
                                                                        gTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGcagcag                                                                >  1:129435/1‑140 (MQ=255)
                                                                            ggTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGaa                                                            <  1:1541763/140‑1 (MQ=255)
                                                                                  tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgca                                                                    <  2:3987444/126‑1 (MQ=255)
                                                                                  tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTg                                                      >  1:2736583/1‑140 (MQ=255)
                                                                                  tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTg                                                      >  1:1429573/1‑140 (MQ=255)
                                                                                  tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGGGACAGCCAGGGGTTACGACCATACGAAACCAGCGCCTGGACTGGGAGCCCCCCTGGCGGTAAATGGGGCTGACGCAGCAGCGCACGCCTg                                                      >  2:2881318/1‑140 (MQ=255)
                                                                                   ggCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGc                                                     >  2:4161903/1‑140 (MQ=255)
                                                                                   ggCGACATTAACTCGTGCATTGGCTTCAGCGGAGAGACGCCCCCAGAGTTGCAGCCAGGGATTACGACCATACGAAACCAGCTCCAGGACTGTGATGCCCTCTGGCGTGAAATGGTGCTGAGGCAACAACGAAAGACTGc                                                     >  2:1970897/1‑140 (MQ=255)
                                                                                       aCATTAACTCGTGCATTGTCTTCGGCAGAGTGACGCCCCCAGAGTGACAGCGAGGGATTACGACCATACGAAGCCGGCTCCTGGACTGTGCTCACCTGTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGc                                                 >  1:3059992/1‑140 (MQ=255)
                                                                                       aCATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGc                                                 >  2:2958955/1‑140 (MQ=255)
                                                                                       aCATGAACTAGTGACTTGTCTGCAGCCGAGAGACGCCCCCAACGTTACAGGCACGGATTCCGCCGATACGAAAACACCCCCCAAAAGGCAATCCCCTCCCCCGTCAACGGGCGCTGCGGCAGCAGCGAAAGCCTGcacgc                                                 >  2:2031689/1‑136 (MQ=255)
                                                                                        cATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGcc                                                >  1:510854/1‑140 (MQ=255)
                                                                                                  tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGcgcg                                      >  1:4063184/1‑140 (MQ=255)
                                                                                                  tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGcgcg                                      <  2:305443/140‑1 (MQ=255)
                                                                                                  tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTCACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGcgcg                                      >  1:3196278/1‑140 (MQ=255)
                                                                                                   gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGa                                     <  1:2598290/140‑1 (MQ=255)
                                                                                                   gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGCGATCGCCTCTGGCGGTAAATGGTGCTGAGGCAGCCGCAAAAGCCTGCGGGCCAACTGGcgcgc                                     >  2:3377049/1‑139 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  1:2525288/140‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  2:2534572/140‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  1:4749752/140‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  1:4430421/140‑1 (MQ=255)
                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCCCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                    <  2:2625426/140‑1 (MQ=255)
                                                                                                     aTTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATg                                   <  1:2138194/140‑1 (MQ=255)
                                                                                                      ttGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGAGCTGAGGCAGCATCGAAAGCCTGCGGGCCAACTGGCGCGATga                                  <  2:4665999/140‑1 (MQ=255)
                                                                                                      ttGTCTGCAGCGGAGAGACGCCCCCAGCGTGACAGCCAGGGATTACGCCCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATga                                  <  2:674528/140‑1 (MQ=255)
                                                                                                        gTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATgaga                                <  1:1911152/140‑1 (MQ=255)
                                                                                                         tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                               <  2:2073739/140‑1 (MQ=255)
                                                                                                          cTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                              <  1:348279/140‑1 (MQ=255)
                                                                                                          cTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                              <  1:2785923/140‑1 (MQ=255)
                                                                                                           ttCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCa                             <  1:4799628/140‑1 (MQ=255)
                                                                                                           ttCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCa                             >  1:2217598/1‑140 (MQ=255)
                                                                                                           tgcAGCGGACAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCa                             <  2:4836095/138‑1 (MQ=255)
                                                                                                            tCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCat                            <  1:1258130/140‑1 (MQ=255)
                                                                                                            tCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCat                            <  2:1255314/140‑1 (MQ=255)
                                                                                                             cAGCGGCGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGCCACCCGCCCCTGCACTGTGATCCCCCCGGGCGTTAACCGGCGCGGACGCAGCAGCGAAACCCTGCGGGCCAGCCGGGGCCGTGAGAGCatc                           >  2:4320808/1‑139 (MQ=255)
                                                                                                                cGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATAttt                        >  1:2903575/1‑140 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATg                     >  1:1991223/1‑140 (MQ=255)
                                                                                                                   agagaCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATg                     >  1:659688/1‑140 (MQ=255)
                                                                                                                            cccAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCg            <  1:3333591/140‑1 (MQ=255)
                                                                                                                              cAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGcc          <  2:2294534/140‑1 (MQ=255)
                                                                                                                              cAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGcc          <  1:4144278/140‑1 (MQ=255)
                                                                                                                               cgagTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCCGGACTGTGATCCCCCCTGGCGGTAAATGGTGCAGAGGCAGCAGCGAAAGCCTGCGGGCCAACGGGCGCGATGAGAGCATATTTATGAGGTCAACGCCg         >  2:2661279/2‑140 (MQ=255)
                                                                                                                               agagTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCg         >  1:3619856/1‑140 (MQ=255)
                                                                                                                                    gACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGaaa    <  1:1203300/140‑1 (MQ=255)
                                                                                                                                     aCAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAAt   >  2:3703197/1‑140 (MQ=255)
                                                                                                                                     aCAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAAt   >  1:2801881/1‑140 (MQ=255)
                                                                                                                                     aCAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAAt   >  2:5065225/1‑140 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  1:2064974/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  1:970459/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  2:1127807/140‑1 (MQ=255)
                                                                                                                                      cAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCACCCTGGCGTTACATGGTGATGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCCATGAGAGCATATTTATGGGATTATCGCCGAGAAATa  <  2:2139622/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: