Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 3,791,922 T→G 5.7% E126D (GAA→GAC EO946_RS19535 ← DNA‑3‑methyladenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0349433,791,9220TG 5.7% 935.2 / 7.7 401E126D (GAA→GACEO946_RS19535DNA‑3‑methyladenine glycosylase
Reads supporting (aligned to +/- strand):  ref base T (160/200);  new base G (3/37);  total (164/237)
Fisher's exact test for biased strand distribution p-value = 1.82e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 433 total aligned reads displayed.

TGCCGATGCGCGGGCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACGT  >  NZ_CP034943/3791783‑3792061
                                                                                                                                           |                                                                                                                                           
 gccgatgcGCGGGCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCAtt                                                                                                                                            >  1:361304/1‑140 (MQ=255)
          cgGGCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGAt                                                                                                                                   >  2:353002/1‑140 (MQ=255)
          cgGGCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGAt                                                                                                                                   >  2:3270609/1‑140 (MQ=255)
          cgGGCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGAt                                                                                                                                   >  1:1677880/1‑140 (MQ=255)
          cgGGCCGGTTGAAAACTCTTCCGGTGTATACCCGCTTTCGATGTAAACCGGCTGTTACGTTAACCAGCGTCAGTAGTCATACAACGAGACGACAAGCTCTATCGTCAGATTTCCGGGCCCGTTTGTCCATTCCATCGGAt                                                                                                                                   >  2:1841997/1‑140 (MQ=255)
           gggCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATg                                                                                                                                  <  2:3384531/140‑1 (MQ=255)
            ggCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGc                                                                                                                                 >  2:3122428/1‑140 (MQ=255)
             gCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTGTCCGGGCCCGTTTGTCCCTTCCCTCGGATGcc                                                                                                                                >  2:191378/1‑140 (MQ=255)
                ggTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGccggc                                                                                                                             <  1:3443862/140‑1 (MQ=255)
                 gTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGccggcc                                                                                                                            >  2:3434172/1‑140 (MQ=255)
                 gTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGccggcc                                                                                                                            >  2:2419919/1‑140 (MQ=255)
                            ttCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCCAGCGCTTTCGTCAGCTTTCCGGGCCCGGTTGTCCATTCCCTCGGATGCCGGCCCGGCCTtctct                                                                                                                 >  2:3616814/1‑140 (MQ=255)
                            ttCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTtctct                                                                                                                 <  1:2301631/140‑1 (MQ=255)
                            ttCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTtctct                                                                                                                 >  1:2607902/1‑140 (MQ=255)
                            ttCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTtctct                                                                                                                 >  2:1297120/1‑140 (MQ=255)
                            ttCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTtctct                                                                                                                 <  2:751620/140‑1 (MQ=255)
                            ttCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTtctct                                                                                                                 <  1:1030355/140‑1 (MQ=255)
                              ccGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTctcc                                                                                                               <  1:2348823/140‑1 (MQ=255)
                               cGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTctcct                                                                                                              <  2:3005812/140‑1 (MQ=255)
                               cGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTctcct                                                                                                              <  2:887282/140‑1 (MQ=255)
                                    tataCCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCAtt                                                                                                         >  1:816573/1‑140 (MQ=255)
                                    tataCCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCAtt                                                                                                         >  1:2486558/1‑140 (MQ=255)
                                          cGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCGCCTCCATTAACAGc                                                                                                   <  2:695776/140‑1 (MQ=255)
                                            cTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCCTTCATCGTGACGACAATCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGCTGCCGGCCCGGCCTGCTCTCCTCCATTAACAGCTg                                                                                                 >  2:1077839/1‑140 (MQ=255)
                                            cTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTg                                                                                                 >  1:2230422/1‑140 (MQ=255)
                                            cTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTg                                                                                                 >  1:1565879/1‑140 (MQ=255)
                                            cTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGCCCTTCTCTCCTCCATTAACAGCTg                                                                                                 >  1:3138927/1‑140 (MQ=255)
                                            cTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCAAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCATCGGATGCCGGCCCGGCCTTCTCTCCTCCAATAACACCTg                                                                                                 >  1:325836/1‑140 (MQ=255)
                                              ttCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGc                                                                                               >  1:1048707/1‑140 (MQ=255)
                                              ttCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGc                                                                                               <  1:1601996/140‑1 (MQ=255)
                                              ttCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGc                                                                                               <  1:306961/140‑1 (MQ=255)
                                              ttCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGc                                                                                               <  1:3296840/140‑1 (MQ=255)
                                              ttCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGc                                                                                               >  2:640524/1‑140 (MQ=255)
                                                     tAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCAtt                                                                                                         <  2:2159040/123‑1 (MQ=255)
                                                             cTGTTCCGTGATCCAGCGTCCGTAGTCATTCAGCGTGACGCCAAGCGCTTTCGTCCGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCaa                                                                                >  1:3662684/1‑140 (MQ=255)
                                                                 tCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGc                                                                            >  1:99734/1‑140 (MQ=255)
                                                                  ccGTGCTCCCGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTCTCCGGGCCCGGGGGTCCAGTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACACCTGGCCGTCATGAGGCTCAATCGcc                                                                           <  2:283372/140‑1 (MQ=255)
                                                                  ccGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCACGCGCTTTCGTCAGCTGTCCGGGCCCGTTTGTCCATTCCCGCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           <  2:3682197/140‑1 (MQ=255)
                                                                  ccGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           <  1:2541325/140‑1 (MQ=255)
                                                                  ccGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           <  1:175123/140‑1 (MQ=255)
                                                                  ccGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           <  2:1780500/140‑1 (MQ=255)
                                                                  ccGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           <  2:2510202/140‑1 (MQ=255)
                                                                  ccGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           <  1:1300482/140‑1 (MQ=255)
                                                                  ccGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTGTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           <  2:1695997/140‑1 (MQ=255)
                                                                   cGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           <  1:2151990/139‑1 (MQ=255)
                                                                     tGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCa                                                                                          >  1:1821720/1‑122 (MQ=255)
                                                                     tGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCa                                                                                          >  1:352504/1‑122 (MQ=255)
                                                                     tGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTg                                                                        >  1:1265796/1‑140 (MQ=255)
                                                                            gCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATa                                                                 >  2:30967/1‑140 (MQ=255)
                                                                            gCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCGCAATCGCCCTGATTAATa                                                                 >  2:2691436/1‑140 (MQ=255)
                                                                                  gTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCtt                                                           >  2:2656540/1‑140 (MQ=255)
                                                                                       cattcatCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCggg                                                      <  1:1530910/140‑1 (MQ=255)
                                                                                        attcatCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGc                                                             <  1:1133335/132‑1 (MQ=255)
                                                                                         ttcatCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCaa                                                                                >  2:381865/1‑112 (MQ=255)
                                                                                         ttcatCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGac                                                    >  1:2832723/1‑140 (MQ=255)
                                                                                              cGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCACTCGCCCGGGTTAATACCGCTTGCGGGACACCTg                                               >  2:248989/1‑140 (MQ=255)
                                                                                              cGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTg                                               >  2:827445/1‑140 (MQ=255)
                                                                                                tGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGcc                                                                           >  2:1612953/1‑110 (MQ=255)
                                                                                                tGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCt                                             >  1:446788/1‑140 (MQ=255)
                                                                                                  aCGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGc                                                                                               >  2:2212633/1‑88 (MQ=255)
                                                                                                        aGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGc                                     >  1:633920/1‑140 (MQ=255)
                                                                                                          cgcTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTCACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCaa                                   <  1:2739404/140‑1 (MQ=255)
                                                                                                          cgcTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCaa                                   <  1:954784/140‑1 (MQ=255)
                                                                                                          cgcTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCaa                                   <  2:1133501/140‑1 (MQ=255)
                                                                                                          cgcTTTCGTCAGCTTTCCGGGCCCGTTTGTCCAGTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCaa                                   <  2:2332295/140‑1 (MQ=255)
                                                                                                          cgcTGTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCaa                                   <  1:3744192/140‑1 (MQ=255)
                                                                                                               tCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACa                              >  1:3487284/1‑140 (MQ=255)
                                                                                                                cGTCAGCGTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACAt                             <  1:2928025/140‑1 (MQ=255)
                                                                                                                 gTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACAtt                            <  2:3198700/140‑1 (MQ=255)
                                                                                                                 gTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCGGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACAtt                            <  2:52772/140‑1 (MQ=255)
                                                                                                                    aGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAg                         <  2:1396173/140‑1 (MQ=255)
                                                                                                                    aGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGATTGCGGGACACCTGCTTCTGCGGCAACGACATTCAg                         <  2:397299/140‑1 (MQ=255)
                                                                                                                      ctgtCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGGCTTATCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCa                       <  1:1040192/137‑1 (MQ=255)
                                                                                                                      cTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCa                       <  1:3544544/140‑1 (MQ=255)
                                                                                                                      cTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCa                       <  1:3806332/140‑1 (MQ=255)
                                                                                                                      cTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCa                       <  1:384495/140‑1 (MQ=255)
                                                                                                                      cTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCa                       >  1:1974750/1‑140 (MQ=255)
                                                                                                                      cTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCa                       <  1:687093/140‑1 (MQ=255)
                                                                                                                      cTTTCCGGGCCCGTTTGTCCATTCCCGCGGATGCCGGCCCGGCCGTCTCTCCGCCAGTCACAGCTGGCCTTCATGACGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCa                       <  2:458672/140‑1 (MQ=255)
                                                                                                                      cTTTCCGGGCCCGTTTGTCCAGTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATGAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCa                       <  1:1396575/140‑1 (MQ=255)
                                                                                                                         tCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCgg                                      >  1:3201726/1‑122 (MQ=255)
                                                                                                                                 cGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTGCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCACTGTGTGGgtg            <  2:1474885/140‑1 (MQ=255)
                                                                                                                                   tttGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGtgtg          >  2:2878937/1‑140 (MQ=255)
                                                                                                                                      ggcccttCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCAt       <  2:804458/135‑1 (MQ=255)
                                                                                                                                      gTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCAt       <  1:1002345/140‑1 (MQ=255)
                                                                                                                                      gTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCAt       <  1:755738/140‑1 (MQ=255)
                                                                                                                                        ccATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCa     <  2:3510177/140‑1 (MQ=255)
                                                                                                                                        ccATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCa     <  2:3349095/140‑1 (MQ=255)
                                                                                                                                        ccATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCa     <  1:1439845/140‑1 (MQ=255)
                                                                                                                                        ccATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCa     <  1:2208420/140‑1 (MQ=255)
                                                                                                                                        ccATGCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCa     <  1:3425127/140‑1 (MQ=255)
                                                                                                                                        ccAGTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCa     <  2:3744329/140‑1 (MQ=255)
                                                                                                                                        ccAGTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCa     <  1:2434867/140‑1 (MQ=255)
                                                                                                                                          aTTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACg   <  1:2663618/140‑1 (MQ=255)
                                                                                                                                           ttCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTCACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACgt  <  1:1933732/140‑1 (MQ=255)
                                                                                                                                           ttCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACgt  <  1:1184021/140‑1 (MQ=255)
                                                                                                                                           ttCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACgt  <  2:1565879/140‑1 (MQ=255)
                                                                                                                                           ttCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACgt  <  1:868270/140‑1 (MQ=255)
                                                                                                                                           ttCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACgt  <  1:50109/140‑1 (MQ=255)
                                                                                                                                           ttCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCGGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACgt  <  1:2393664/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TGCCGATGCGCGGGCCGGTTGAAATGTCTTCCGGTGTATACCCGCTTTCGATGTAAAGCGGCTGTTCCGTGATCCAGCGTCCGTAGTCATTCATCGTGACGCCAAGCGCTTTCGTCAGCTTTCCGGGCCCGTTTGTCCATTCCCTCGGATGCCGGCCCGGCCTTCTCTCCTCCATTAACAGCTGGCCTTCATGAGGCTCAATCGCCCTGATTAATACCGCTTGCGGGACACCTGCTTCTGCGGCAACGACATTCAGCAATGTGTGGGTGTGCATCACGT  >  NZ_CP034943/3791783‑3792061

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