Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,527,329 A→G L28P (CTG→CCG)  nanS ← 9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,527,3290AG95.6% 301.6 / 2.4 91L28P (CTG→CCG) nanS9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (34/53);  total (36/55)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCG  >  NZ_CP009273/4527196‑4527462
                                                                                                                                     |                                                                                                                                     
aaCATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGGTGCGTATGCGCCATTCTCGCTAATTGTTTCATTCTGGGTTGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTGCt                                                                                                                                  >  2:13070/1‑139 (MQ=255)
aaCATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCt                                                                                                                                  >  1:135622/1‑139 (MQ=255)
aaCATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCt                                                                                                                                  <  2:701553/139‑1 (MQ=255)
 aCATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCtt                                                                                                                                 >  2:475455/1‑139 (MQ=255)
 aCATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCtt                                                                                                                                 <  2:1166202/139‑1 (MQ=255)
 aCATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCtt                                                                                                                                 <  1:929138/139‑1 (MQ=255)
   aTCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAGGTGACATGGCGGGCCTTCGGGAAGCAAGTGCGCAAATCGAGCTAATTGTTTCATTCTGGGAAGAGGCACCTGTTCCCGATCCGGTAGTGGAAGGCCTTCg                                                                                                                               >  2:573169/1‑139 (MQ=255)
     cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCtt                                                                                                                                 >  1:1123768/1‑135 (MQ=255)
     cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCtt                                                                                                                                 <  2:1123768/135‑1 (MQ=255)
     cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTTGcc                                                                                                                             >  1:58324/1‑139 (MQ=255)
     cGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTTGcc                                                                                                                             >  2:1199028/1‑139 (MQ=255)
      gtgtGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCa                                                                                                                            >  1:784519/1‑139 (MQ=255)
       tgtgGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTTGCCa                                                                                                                            >  2:933917/1‑138 (MQ=255)
         tgGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATa                                                                                                                          <  2:161556/138‑1 (MQ=255)
         tgGGAAGTGAGTCAGTGTAATAATGTCGTTAAAGTGACATGGCGGCCCTCCGGGATGCGTATGCTCAAATCTCGCTAATTGTTTAACTCTTGGATCAGGCGCATCTTCCCGGTCCGGTAATGGTGGTCCTTCGccctt                                                                                                                          <  2:1205423/138‑4 (MQ=255)
             cAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  1:621727/139‑1 (MQ=255)
             cAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  2:691698/139‑1 (MQ=255)
             cAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  1:719434/139‑1 (MQ=255)
             cAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  2:58324/139‑1 (MQ=255)
             cAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  1:1169378/139‑1 (MQ=255)
             cAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  1:863492/139‑1 (MQ=255)
             cAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  1:1199028/139‑1 (MQ=255)
                tGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTTGGCCTACACCATgg                                                                                                                  >  2:1007697/1‑139 (MQ=255)
                 gAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATgg                                                                                                                  >  1:121808/1‑138 (MQ=255)
                      aGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCttcgccattcgcca                                                                                                                     >  2:1288692/1‑130 (MQ=255)
                      aGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  1:1288692/130‑1 (MQ=255)
                            ataatGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGAc                                                                                                      <  2:115897/139‑1 (MQ=255)
                             taatGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  2:154591/123‑1 (MQ=255)
                             taatGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     >  1:154591/1‑123 (MQ=255)
                             taatGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCCCCCTACGCCCTGGCATTATACTGAcc                                                                                                     >  2:490504/1‑139 (MQ=255)
                                  tCGTTAAAGTGACATGGCGGGCCTCCGGGATTCGTATGCGCAAATCTCTCTAATTGTGTAATTCTGGGATGAGGCTCATCTTCCCTGTCCGGTAATGGCACTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                                                 <  1:573169/138‑1 (MQ=255)
                                  tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  2:946848/118‑1 (MQ=255)
                                  tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     >  1:946848/1‑118 (MQ=255)
                                  tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                                                 <  2:1027594/138‑1 (MQ=255)
                                  tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                                                 <  2:501523/138‑1 (MQ=255)
                                  tCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                                                 <  2:682331/138‑1 (MQ=255)
                                    gTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACg                                                                                              >  2:85266/1‑139 (MQ=255)
                                     ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCCTACGCCATGGCATTAGACTTACCCGCAACgg                                                                                             >  2:1009533/1‑139 (MQ=255)
                                     ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     <  2:418582/115‑1 (MQ=255)
                                     ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                                                 >  1:138259/1‑135 (MQ=255)
                                     ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCa                                                                                                 <  2:138259/135‑1 (MQ=255)
                                     ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACgg                                                                                             >  1:1216878/1‑139 (MQ=255)
                                     ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACgg                                                                                             >  1:659466/1‑139 (MQ=255)
                                     ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGa                                                                                                           <  2:696425/125‑1 (MQ=255)
                                     ttAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGa                                                                                                           >  1:696425/1‑125 (MQ=255)
                                     ttAAAGAGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCa                                                                                                                     >  1:418582/1‑115 (MQ=255)
                                       aaaGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGt                                                                                            <  2:121808/138‑1 (MQ=255)
                                        aaGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTaa                                                                                          <  1:457982/139‑1 (MQ=255)
                                         aGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACGGACCAGCAACGGTAAg                                                                                         <  1:13070/139‑1 (MQ=255)
                                           tGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGaa                                                                                       <  1:187655/139‑1 (MQ=255)
                                           tGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGaa                                                                                       <  2:186727/139‑1 (MQ=255)
                                             aCATGGCGGGCCTCCGGGATGCGTATGCGCAATTCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCAAACGCCATTGGAATAGACTGGCCAGCAACGGTAAGAACa                                                                                     >  1:1117874/1‑139 (MQ=255)
                                             aCATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACa                                                                                     <  1:277198/139‑1 (MQ=255)
                                              cATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAg                                                                                         <  1:95432/134‑1 (MQ=255)
                                              cATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAg                                                                                         >  2:95432/1‑134 (MQ=255)
                                              cATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACAt                                                                                    <  2:1219949/139‑1 (MQ=255)
                                                      gCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGtaataa                                                                            <  2:135622/139‑1 (MQ=255)
                                                           cGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                                       <  1:348037/139‑1 (MQ=255)
                                                           cGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                                       <  1:995612/139‑1 (MQ=255)
                                                           cGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                                       <  2:93433/139‑1 (MQ=255)
                                                           cGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                                       <  2:554822/139‑1 (MQ=255)
                                                           cGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                                       <  1:407802/139‑1 (MQ=255)
                                                           cGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCg                                                                         <  2:855242/137‑1 (MQ=255)
                                                           cGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGAAGAACATAGTAATAATCggg                                                                       <  1:304764/138‑1 (MQ=255)
                                                                tGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGata                                                                  <  1:85266/139‑1 (MQ=255)
                                                                 gCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTTGCCATACGCCATGGCATTAGACTGACCAGCAACGGGAAGAAAATAGTAATAAACGGGCGatat                                                                 >  2:77277/1‑139 (MQ=255)
                                                                    tATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGGAAGAACATAGTAATAATCGGGCGAtattat                                                              >  2:104208/1‑139 (MQ=255)
                                                                        cgcAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTaataat                                                                           <  1:578369/122‑1 (MQ=255)
                                                                        cgcAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTaataat                                                                           >  2:578369/1‑122 (MQ=255)
                                                                        cgcAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCg                                                          <  1:77277/139‑1 (MQ=255)
                                                                            aaTCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGcc                                                                                                                             <  2:223948/68‑1 (MQ=255)
                                                                            aaTCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGcc                                                                                                                             >  1:223948/1‑68 (MQ=255)
                                                                            aaTCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATa                                                                                                                          <  2:210306/71‑1 (MQ=255)
                                                                            aaTCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATa                                                                                                                          >  1:210306/1‑71 (MQ=255)
                                                                                                    ggATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAAt                               <  2:784519/138‑1 (MQ=255)
                                                                                                    ggATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAAt                               <  1:294228/138‑1 (MQ=255)
                                                                                                    ggATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAAt                               >  2:1019551/1‑138 (MQ=255)
                                                                                                    ggATGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAAt                               >  1:688559/1‑138 (MQ=255)
                                                                                                      aTGAGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGAc                            <  1:906822/139‑1 (MQ=255)
                                                                                                         aGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCCTGGCATTAGACTGACCAGCAACGGTAAGAACa                                                                                     >  1:356554/1‑79 (MQ=255)
                                                                                                         aGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACa                                                                                     <  2:356554/79‑1 (MQ=255)
                                                                                                         aGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                                       >  1:400040/1‑93 (MQ=255)
                                                                                                         aGGCGCATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCggg                                                                       <  2:400040/93‑1 (MQ=255)
                                                                                                              cATCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGa                     >  1:501105/1‑138 (MQ=255)
                                                                                                                tCTTCCCTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAg                  >  2:1281582/1‑139 (MQ=255)
                                                                                                                    cccTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGc                 <  1:490504/136‑1 (MQ=255)
                                                                                                                    cccTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAg               <  1:631919/138‑1 (MQ=255)
                                                                                                                     ccTGTCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAgg              <  2:659466/138‑1 (MQ=255)
                                                                                                                         tCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAggg             >  2:829784/1‑135 (MQ=255)
                                                                                                                         tCCGGTAATGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAggg             <  1:829784/135‑1 (MQ=255)
                                                                                                                                aTGGCGGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCg  <  1:1007697/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
AACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGGCGGGCCTCCGGGATGCGTATGCGCAAATCTCGCTAATTGTTTAATTCTGGGATGAGGCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCG  >  NZ_CP009273/4527196‑4527462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: