Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 1,576,151 T→G I238I (ATA→ATC gadB ← glutamate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252681,576,1510TG98.8% 183.8 / NA 84I238I (ATA→ATCgadBglutamate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (43/40);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATG  >  NZ_CP025268/1576084‑1576286
                                                                   |                                                                                                                                         
cccAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCa                                                                   <  1:2982135/140‑1 (MQ=11)
cccAGACGACATCCGGTGCGACGAACGGTGCGAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACTGGTATCGGCCTGAAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCa                                                                   <  1:2982136/140‑1 (MQ=11)
 ccAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑g                                                                >  2:3445197/1‑140 (MQ=14)
  cAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTAACA‑‑gt                                                               >  2:2055194/1‑140 (MQ=17)
   aGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑gtg                                                              >  2:987086/1‑140 (MQ=18)
        atatCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGt                                                         <  2:1269056/140‑1 (MQ=18)
        atatCCGGGGCGACGAAAGGTGCCAGGACGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGt                                                         <  2:1268044/140‑1 (MQ=12)
         tatCCGGGGTGACGAACGGTGCCAGGAAGCCGCCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACTA‑‑GTGTAGGTc                                                        >  2:1136933/1‑140 (MQ=255)
         tatCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGCTGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTc                                                        >  2:1136828/1‑140 (MQ=9)
         tatCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTc                                                        >  1:99009/1‑140 (MQ=18)
         tatCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTc                                                        >  1:5327977/1‑140 (MQ=18)
         tatCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTc                                                        >  2:4910726/1‑140 (MQ=14)
         tatCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTc                                                        >  1:2173158/1‑140 (MQ=18)
         tatCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTc                                                        >  1:4490832/1‑140 (MQ=18)
         tagccGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGTGAACTCATAGTTACCA‑‑GTGTAGGTc                                                        >  1:418283/4‑140 (MQ=9)
          atCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCa                                                       >  2:1774758/1‑140 (MQ=17)
           tCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATGCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCAc                                                      >  2:5345933/1‑140 (MQ=12)
             cGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTCGGTCACGc                                                    <  2:3361510/140‑1 (MQ=18)
               gggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCg                                                  >  2:1827933/1‑140 (MQ=18)
               gggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCg                                                  >  1:2160838/1‑140 (MQ=18)
               gggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCg                                                  >  1:4508268/1‑140 (MQ=14)
               gggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCg                                                  <  2:4481662/140‑1 (MQ=18)
               gggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCg                                                  >  1:4085501/1‑140 (MQ=18)
               gggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCg                                                  >  1:3084733/1‑140 (MQ=18)
                ggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGa                                                 >  1:435415/1‑140 (MQ=17)
                ggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGa                                                 >  1:5402078/1‑140 (MQ=18)
                ggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGa                                                 <  2:1868629/140‑1 (MQ=18)
                  cgacgaACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGaaa                                               >  1:385813/1‑140 (MQ=18)
                   gacgaACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAg                                              >  2:2565865/1‑140 (MQ=18)
                     cgaACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTc                                            >  1:4868917/1‑140 (MQ=18)
                     cgaACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTc                                            >  2:888892/1‑140 (MQ=18)
                          ggTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGcac                                       <  2:3855928/140‑1 (MQ=17)
                            tGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGcacca                                     >  1:2366682/1‑140 (MQ=17)
                              ccAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACg                                   <  1:4287453/140‑1 (MQ=18)
                               cAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGc                                  <  2:3868383/140‑1 (MQ=17)
                                  gAAGCTACCGCGGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGa                               <  2:1123748/140‑1 (MQ=12)
                                    aGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATg                             <  1:1462442/140‑1 (MQ=18)
                                    aGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATg                             <  2:4941289/140‑1 (MQ=18)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGCTGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGCACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:4653860/1‑140 (MQ=255)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCAGTGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATg                             >  1:2589018/1‑140 (MQ=11)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  2:4280310/1‑140 (MQ=18)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  2:658590/1‑140 (MQ=17)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:1954259/1‑140 (MQ=18)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  2:241360/1‑140 (MQ=14)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:1980018/1‑140 (MQ=17)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:5390275/1‑140 (MQ=18)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:572491/1‑140 (MQ=17)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:4563306/1‑140 (MQ=18)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:4378944/1‑140 (MQ=18)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:4102947/1‑140 (MQ=18)
                                      ccACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgt                           >  1:2325228/1‑140 (MQ=18)
                                      ccACCGCTGGCAGCGTCGATGTGAATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGGCGAAAGTCGGCACCACACCAATggg                           >  2:2225246/1‑139 (MQ=11)
                                        aCCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGgtgt                         >  1:2730495/1‑140 (MQ=18)
                                        aCCGCTGGCAGCGTAGATGTGGATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGGTGTGGGAACTCATAGTCACCA‑‑GTGGAGGTCACGCCGAAAGTCGGCACCACGCCGATGgtgt                         >  2:3337395/1‑140 (MQ=11)
                                         ccGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGtt                        <  2:312112/140‑1 (MQ=18)
                                         ccGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGtt                        <  2:3522452/140‑1 (MQ=18)
                                              ggCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGt                   >  2:1218806/1‑140 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:5402078/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:4531443/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:4531442/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:5577470/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:418283/140‑1 (MQ=17)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:2181352/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:3877044/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:3384007/140‑1 (MQ=17)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:1637112/140‑1 (MQ=17)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:2386446/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:2473317/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:2473316/140‑1 (MQ=18)
                                                cAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:2384774/140‑1 (MQ=18)
                                                cAGAGTCGATGTGAAAGTCGATCTAGCTACCGGTATCGGCCTGGACTTTATCCAGCGCATAGTACAGCGGTTGTGGGAACTCATAGTTACAA‑‑GTGTAGGTAACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTca                 <  2:4278482/140‑1 (MQ=11)
                                                 aGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTcac                <  2:5136472/140‑1 (MQ=18)
                                                  gCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCCTAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTcaca               <  1:2398373/140‑1 (MQ=18)
                                                  gCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTcaca               <  2:4798208/140‑1 (MQ=14)
                                                  gCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTcaca               <  2:2430845/140‑1 (MQ=18)
                                                  gCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTcaca               <  2:2446944/140‑1 (MQ=18)
                                                  gCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTcaca               <  2:5106045/140‑1 (MQ=14)
                                                  gCGTCGATGTGCACGTCGATGTCGATACCGGTATCGGCCTGGCATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAATTCATAGTTACCA‑‑GTGTAGGTCCCACCGAAAGTCGGCACCACGCCGCTGGTGTTTTCGTcaca               <  2:4798222/140‑1 (MQ=11)
                                                     ccGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGc            <  2:2118683/139‑1 (MQ=18)
                                                      cGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCt           <  2:5571295/140‑1 (MQ=18)
                                                         tgtgCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCa        >  2:942877/1‑140 (MQ=18)
                                                          gtgCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCaa       <  1:564496/140‑1 (MQ=17)
                                                             cATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAGGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCa    <  2:5566821/140‑1 (MQ=12)
                                                              atgtcgatgtcgATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCAt   >  1:2838952/1‑140 (MQ=17)
                                                               tgtcgatgtcgaTACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGTATGTCATAGTTATCA‑‑GTGTAGGTCACGGTGATAGTCGGCACCACGCCGATGGTGTTTTTGTGATACTCTTCaacgatt  >  2:5334922/1‑135 (MQ=11)
                                                               tgtcgatgtcgaTACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATg  <  2:2131647/140‑1 (MQ=18)
                                                               tctCTATGTCGGTATCGGTATCGGCCCTTAATTTATCCAGCGCATTGTTCAGCCGTTGTGGGACATCATACTTACCA‑‑GTGTATGTCACGCAGACAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATg  <  2:4486277/138‑1 (MQ=11)
                                                                   |                                                                                                                                         
CCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCA‑‑GTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATG  >  NZ_CP025268/1576084‑1576286

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: