Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 3,971,997 A→G pseudogene (1048/1048 nt) wzzE → ECA polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252683,971,9970AG99.6% 914.0 / NA 275pseudogene (1048/1048 nt)wzzEECA polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (156/117);  total (156/119)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

CCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCC  >  NZ_CP025268/3971860‑3972132
                                                                                                                                         |                                                                                                                                       
ccTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATagcaa                                                                                                                                      <  2:355316/140‑1 (MQ=255)
ccTACCGCTATTTGCGTACGCCGGAAGAACCGGTAAATCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCAGTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATagcaa                                                                                                                                      <  1:4479678/140‑1 (MQ=255)
  tcccGCTATTTGCGCTCGCCGGAAGAACCGGTACACCGCGATAGCCTACGTCGTGCCTGCCTGATGAGTATGCGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACGCGCCGTTGCTCGAAATAgcaaca                                                                                                                                    <  2:4438639/137‑1 (MQ=255)
    tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACAct                                                                                                                                  >  1:1455095/1‑140 (MQ=255)
    tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACAct                                                                                                                                  >  2:1576049/1‑140 (MQ=255)
    tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACAct                                                                                                                                  >  1:786735/1‑140 (MQ=255)
       cTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACActgct                                                                                                                               <  1:1278488/140‑1 (MQ=255)
         aTTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACtgctgc                                                                                                                             >  2:2783573/1‑140 (MQ=255)
            tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGt                                                                                                                          <  1:1310625/140‑1 (MQ=255)
            tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGt                                                                                                                          <  1:2527315/140‑1 (MQ=255)
            tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGt                                                                                                                          <  2:4789909/140‑1 (MQ=255)
            tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGt                                                                                                                          <  1:2359785/140‑1 (MQ=255)
              cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGa                                                                                                                        >  1:744547/1‑140 (MQ=255)
              cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGa                                                                                                                        >  2:3419957/1‑140 (MQ=255)
               gTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAg                                                                                                                       <  2:3668500/140‑1 (MQ=255)
                        aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGgc                                                                                                              >  1:4153996/1‑140 (MQ=255)
                           aaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                                                                           <  2:2647686/140‑1 (MQ=255)
                              cgctAAAGCGCGATAGCCCACGTCGTGCCTTCCTGATGAATATGTGGGGCATTGTCGGGGGGCTGAGCGGGGCTGGTGTCGAATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGTTcgc                                                                                                        <  2:1629464/137‑1 (MQ=255)
                                 tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccgc                                                                                                     >  2:3223129/1‑140 (MQ=255)
                                          aTAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGaa                                                                                            <  2:1448978/140‑1 (MQ=255)
                                            aGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAaga                                                                                          <  2:2698233/140‑1 (MQ=255)
                                                   gtcgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccgcag                                                                                                   >  2:4233438/1‑122 (MQ=255)
                                                   gtcgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                                                   <  1:3136108/140‑1 (MQ=255)
                                                       tGCCTTCCTGGTTATTCTGTGGGCCATTCTCGGGGAGCTGATCGGGGCTGGTGTCGCATTAGCCCGCCGGTGCGCGAAATAGCAACACTGCTGCGGTGCGCCCAAAGGCGCTCGCCGGTTATTCGAAGAGAATCGAtgtg                                                                               <  2:3175873/140‑1 (MQ=255)
                                                       tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                                                                               >  1:809157/1‑140 (MQ=255)
                                                       tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                                                                               <  2:1004092/140‑1 (MQ=255)
                                                       tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                                                                               >  1:245733/1‑140 (MQ=255)
                                                       tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                                                                               <  1:4862227/140‑1 (MQ=255)
                                                       tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                                                                               >  1:506148/1‑140 (MQ=255)
                                                       tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCACTTATTCGAAGAGAATCGAtgtg                                                                               >  1:4676533/1‑140 (MQ=255)
                                                        gCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                                                                              <  1:3488880/140‑1 (MQ=255)
                                                        gCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                                                                              >  1:4620375/1‑140 (MQ=255)
                                                        gCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                                                                              >  1:3707048/1‑140 (MQ=255)
                                                        gCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                                                                              >  1:1571645/1‑140 (MQ=255)
                                                         ccTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCa                                                                                                                   >  1:1054032/1‑102 (MQ=255)
                                                             ccTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa                                                                         <  1:3667792/140‑1 (MQ=255)
                                                             ccTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa                                                                         >  1:589394/1‑140 (MQ=255)
                                                             ccTGATGATTATGTGCGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa                                                                         <  2:1110055/140‑1 (MQ=255)
                                                                gatgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                                                      <  2:2387143/140‑1 (MQ=255)
                                                                gatgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                                                      <  1:2942481/140‑1 (MQ=255)
                                                                gatgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                                                      <  2:4197560/140‑1 (MQ=255)
                                                                gatgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCAGTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGATTCGATGTGAAAGTactg                                                                      <  2:4078469/140‑1 (MQ=255)
                                                                 atgatTATGTGGGGCATTGTCGGGGTGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                                     <  2:809157/140‑1 (MQ=255)
                                                                 atgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAATAGAATCGATGTGAAAGTactga                                                                     <  2:4047425/140‑1 (MQ=255)
                                                                 atgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                                     <  2:3050005/140‑1 (MQ=255)
                                                                 atgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                                     <  2:941504/140‑1 (MQ=255)
                                                                 atgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                                     <  1:2269272/140‑1 (MQ=255)
                                                                 atgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                                     <  1:1977159/140‑1 (MQ=255)
                                                                 atgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCATTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCTCGTATTCTAAGACAAGCGAGGTGAAAGGactga                                                                     >  2:1086448/1‑140 (MQ=255)
                                                                        tgtggggCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTAtt                                                              >  1:3809385/4‑140 (MQ=255)
                                                                        tgtggggCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTAtt                                                              <  1:1591909/137‑1 (MQ=255)
                                                                             ggCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa                                                         >  2:2013344/1‑140 (MQ=255)
                                                                             ggCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa                                                         <  2:1166115/140‑1 (MQ=255)
                                                                             ggCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa                                                         >  1:3568849/1‑140 (MQ=255)
                                                                               gaTTGTCGGGGGGCTGATCGGGGCTGGGGCCGCATTAACCCGCCGGTGCTCGAAATAGCAAGACTGCTGCGGTGAGGGCAAAGGCGCTCGCCGCTTATTCGAAGAGAGGCTATGTGAAAGTACGGAATGTATTTGGTAcg                                                       >  2:253799/2‑140 (MQ=255)
                                                                               cATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAcg                                                       <  1:4253253/140‑1 (MQ=255)
                                                                               cATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAcg                                                       >  2:3242136/1‑140 (MQ=255)
                                                                                 ttGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTGGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAGGTGAAAGTACTGGCTGTATTTGGTAcgcg                                                     >  2:2184993/1‑140 (MQ=255)
                                                                                   tgtcggggggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACgcgc                                                    >  1:248845/5‑140 (MQ=255)
                                                                                   tgtcggggggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACgcgc                                                    >  2:3858308/5‑140 (MQ=255)
                                                                                       ggggggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCgg                                                <  1:2160481/140‑1 (MQ=255)
                                                                                          gggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAg                                             >  1:2353616/1‑140 (MQ=255)
                                                                                          gggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAg                                             >  2:4897510/1‑140 (MQ=255)
                                                                                           ggCTGATCGGGGTTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTCAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGc                                            <  2:1765881/140‑1 (MQ=255)
                                                                                             cTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTCGTACGCGCCCGGAAGCCa                                          >  1:4005518/1‑140 (MQ=255)
                                                                                               gATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGGGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATc                                        >  1:1491965/1‑140 (MQ=255)
                                                                                               gATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATc                                        >  1:3647668/1‑140 (MQ=255)
                                                                                               gATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATc                                        >  1:3022673/1‑140 (MQ=255)
                                                                                               gATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATc                                        >  2:1252496/1‑140 (MQ=255)
                                                                                               gATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATc                                        >  2:4087749/1‑140 (MQ=255)
                                                                                               gATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATc                                        >  1:1677031/1‑140 (MQ=255)
                                                                                                 tCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCaa                                      <  2:4952875/140‑1 (MQ=255)
                                                                                                  cGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAg                                     <  1:2076172/140‑1 (MQ=255)
                                                                                                     ggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATg                                  >  1:86847/1‑140 (MQ=255)
                                                                                                     ggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATg                                  >  2:255059/1‑140 (MQ=255)
                                                                                                     ggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATg                                  >  1:3430097/1‑140 (MQ=255)
                                                                                                      gCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATgg                                 <  1:1122551/140‑1 (MQ=255)
                                                                                                      gCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATgg                                 <  1:358645/140‑1 (MQ=255)
                                                                                                      gCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATgg                                 <  1:306386/140‑1 (MQ=255)
                                                                                                      ccTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGACGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAAATgg                                 <  2:4695686/139‑1 (MQ=255)
                                                                                                        tGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg                               <  1:3987767/140‑1 (MQ=255)
                                                                                                        tGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg                               <  1:4417026/140‑1 (MQ=255)
                                                                                                        tGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg                               <  1:1959013/140‑1 (MQ=255)
                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGcc                             >  1:4104789/1‑140 (MQ=255)
                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGcc                             >  2:5109127/1‑140 (MQ=255)
                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGcc                             >  1:61764/1‑140 (MQ=255)
                                                                                                           tgtCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCg                            <  1:3317829/140‑1 (MQ=255)
                                                                                                            gtCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGt                           <  1:3936749/140‑1 (MQ=255)
                                                                                                               gCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTgg                        >  2:2490348/1‑140 (MQ=255)
                                                                                                               gCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTgg                        >  2:1698412/1‑140 (MQ=255)
                                                                                                               gCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTgg                        >  2:4553797/1‑140 (MQ=255)
                                                                                                               gCATTAACCCGCCGTTGCTCGAAATAGAAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTgg                        >  2:1899224/1‑140 (MQ=255)
                                                                                                                         gccgTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGtt              >  1:76056/1‑140 (MQ=255)
                                                                                                                          ccgTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTg             >  1:4938777/1‑140 (MQ=255)
                                                                                                                           cgTTGCTGGAAATAGAAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCATATTCGAAGAGAATGGAGGGGAAAGTACTGACTGCATTCGGTCCGCGCACGGAAGCCATGAAGAGGGCGCTGTTCGTGCCTGCGtttg            >  2:633585/1‑138 (MQ=255)
                                                                                                                           cgTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTgg            >  1:437542/1‑140 (MQ=255)
                                                                                                                               gCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCaaa        >  2:4319733/1‑140 (MQ=255)
                                                                                                                                cTCGAAATAGCCACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACGGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCaaaa       >  1:4231589/1‑140 (MQ=255)
                                                                                                                                    aaaTAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATc   >  1:4999997/1‑140 (MQ=255)
                                                                                                                                     aaTAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGGATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATcc  <  1:944936/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCC  >  NZ_CP025268/3971860‑3972132

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