Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,820,5740CA33.3% 93.3 / 11.6 48G203C (GGC→TGC) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (24/8);  new base A (0/16);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 6.52e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCC  >  NZ_CP009273/1820327‑1820819
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          
gcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCaa                                                                                                                                                                                                                                                          <  1:795723/250‑1 (MQ=255)
          ttCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTgg                                                                                                                                                                                                                                                 <  1:570105/249‑1 (MQ=255)
                       cGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTTTCACGATACCGACACTCTCACCAAATAGCAGCGCAGCAATTTCGCTTGCAGGTAATTCAAACTGGCGAAAAACGGTTTGCCCCAAGGGCCACGCGTTGCCCAGTTCCCGCGTACAGGGCCGCGCGCCCAAAAATCTGGCGCTGCAAAGGGt                                                                                                                                                                                                                                   >  2:70921/1‑249 (MQ=255)
                         gCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAAGGCAACGCTTTGCCCAGTTCCCGCGTACAGGCCAGCGCGCCCAAATGTCTGGCGCTGAAAAGGGGaa                                                                                                                                                                                                                                 >  1:457803/1‑249 (MQ=255)
                              aaaTTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGGTTTCACGATACCGACACCCTCCCCCGATAACAACGCAGCAATTTCGCTGGCAGTTACCTCAAAATGGCGCAGATCGTTTTTCCCAAAAGGCAACGCCTTTCCAAATTCCCACGTCCCGGGCCGCGCGCCCAAATTGCCTGCGCTGAGATGGGTAAACGGac                                                                                                                                                                                                                         >  2:469336/1‑250 (MQ=255)
                                     gCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAgg                                                                                                                                                                                                                  <  1:217161/250‑1 (MQ=255)
                                     gCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAgg                                                                                                                                                                                                                  <  1:36202/250‑1 (MQ=255)
                                                      tGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGGTACCGACACCCTCACCAGATAGCAGCGGCGCAATTGCGCTGGCAGTTTACGCCAACATGCGTAAATCGGTTTTCCCCAAGGGGAAACCTTTGCCAAGTTCCAGCGGACAGGCCAGCGCGCCCAAATGGGTGGCGCCGTAATAGGGAAACGTACCACCCCGGTCCTGATGGGACATCa                                                                                                                                                                                                >  2:159947/1‑250 (MQ=255)
                                                             cGAAGGACGGCGAGTGGCGAGTAAGTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCt                                                                                                                                                                                          <  2:289367/249‑1 (MQ=255)
                                                                     ggcgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCccg                                                                                                                                                                                <  2:537425/250‑1 (MQ=255)
                                                                                                    gATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTTCAAACTGGCGAAGAACGGTTTGCCCAAAGGGCAACGCTTTGTCCAGTTTCCGCGGACAGGGCCGCGCGCCAAAATGTCTCGCGCCGAAATGGGTGAACCTAACAACCGGGTTCTGAATGGAAACCCCCGCCCCCCGGCCCGCCGCCCCCACCCACGTCCGAAAATTCTCgcc                                                                                                                                                 >  2:325793/1‑248 (MQ=255)
                                                                                                      tAACGGAGCGGCGGTGTTCTCACGATACCGCCACTCCCACCAGAAAGGAGCGCAGCAAATGCGCTGGCCGTTACTGCAAAAAGGCGAAGATACTTTTTTCCAAAAGGCAACGCTTTTCCAAGTTCCACCGTACAGGCCAGCGCGCCAAAATATCCGGCGCTGCAAATGGTAAAAGGCACACCACGGTCCCAATGTGACACCAGCGCCTCTCACCCCCCCGCCCCCAACCACGCCACAAATTtctcacccc                                                                                                                                               >  2:385305/1‑245 (MQ=255)
                                                                                                         gggggcggcgGTTTCCCCACGATCCCGACACCCCCACCAGAAAGCAGCGCAGCAATTGCGGTGGCAGTTACTACAAACGGGCGAAGATCTTTTTGCCAAAGGGACAACTCTTTGCAAAGTCCAAGCGAACGGCCCAGCGCGCAAAAATGTCGGGCGCAGAAAGGGGAAACCGCCCCACCGGTTTCCGGAGGAAACCCCGGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAgg                                                                                                                                            <  2:259273/246‑1 (MQ=255)
                                                                                                                gcggTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCCCCCAGCCACGTCCGAAATTTCTCGTCCCAGGGAATgcc                                                                                                                                     >  2:66922/1‑248 (MQ=255)
                                                                                                                        ctcACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAACTTTCTCGTCCCAGGGAATGTCGCGCTGCg                                                                                                                             >  2:174061/1‑250 (MQ=255)
                                                                                                                          cccGACACCGCCCCCCTCCCAAGATAGCAGCGAAGAAATGGCGTGGACGTTTCCGCAAACCGGGGAAAGTCTTTTTTGCCAAAGGGGAACCGCTTTCCAAATTCCCAGCGTACAGGCCGGCGCGCAAAAATGTCGGGCGCAAAAATGGAAAAACGCCCCCCCAGTTCCCTGAGGAACCCCCACCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGt                                                                                                                           <  2:457887/248‑1 (MQ=255)
                                                                                                                                aCCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATaca                                                                                                                     >  1:177864/1‑250 (MQ=255)
                                                                                                                                aCCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGCTGCGGTAATaca                                                                                                                     >  2:53627/1‑250 (MQ=255)
                                                                                                                                 ccGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACGGGCAAAGATCGTTTTGCCAAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCAAAAATGTCTGGCGCTGAAATGGTTAAACGTACCCCCAGTTTCCTGATGAAACCCCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATaca                                                                                                                     <  1:693058/249‑1 (MQ=255)
                                                                                                                                            cccAGAGAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACGGGCGAAGATCGTTTTGCCAAAAGGGCAACGCTTTGCCAAGTCCCAGCGTACAGGCCAGCGCGCAAAAATGTCGGGCGCTAAAATGGGTAAACGTCCCCCCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGc                                                                                                            <  1:291700/246‑1 (MQ=255)
                                                                                                                                                              aaTTGCGCTGGCAGTTACTGCAAACTGGCGCAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAACTGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACCGCGGATGCAAGGAGCCACGCAtt                                                                                       >  2:65188/1‑250 (MQ=255)
                                                                                                                                                              aaTTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATCAACCGAACTGCGGATGCAAGGAGCCACGAAtt                                                                                       >  2:254048/1‑250 (MQ=255)
                                                                                                                                                               aTTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGCTCCTGATGGCACACCAGCGCCTCCAGCCCCCCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGACTGTCGCGCTTCGGTAATACACCGAACTGCGGACGACAGGAGCCACGAATTg                                                                                      >  2:502524/1‑250 (MQ=255)
                                                                                                                                                                ttGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCTCCGTCAGCAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGCACGCAAGGAGCCACGAATTGc                                                                                     >  2:176628/1‑250 (MQ=255)
                                                                                                                                                                        gCAGTTCCTGAAAACTGGCGAAGATCGTTTGGCCAAAAGGGCACCGCTTTGCCAAGTCCCAGCGTACAGGCCAGCGCGCAAAAATGTCTGGCGCTGAAAGGGGTAACCGTCCCCCCAGTTCCCTGAGGGAACACAAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTa                                                                             <  1:664734/250‑1 (MQ=255)
                                                                                                                                                                        gCAGTTACGGCAAGCTGGCAAAGTGCGTTTTGCCCAAGGGCCAACGCTTTGCCAAGTCCCAGCGAACAGCCCCGCGCAAAAAAATGTCCGGCGCGAAAATGGAAAAACGCACCCCCGGGTTCCTGATGAACCACCCGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTa                                                                             <  2:69925/250‑1 (MQ=255)
                                                                                                                                                                            gttgctgaaaCTGGGCAAAGATGTTTTGGCCAAAGGGACACCCTTTTCCAAAGTCCCGGCTACAAGGCCCGCCCCAAAAATGTCCGGGCGTGAAATGGGTAAACCGCCCCCCAGGGTCCTTGAGGAAACCCCCCCCCCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGa                                                                           <  2:626438/242‑1 (MQ=255)
                                                                                                                                                                              aCTGAAAACTGGCGAAGATCGTTTTGCCCAAAGGCCACCGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCAAAAATGTCGGGCGCTAAAATGGGTAAACGTCCCCCCAGGTTCCTGATGAAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCa                                                                        <  1:401333/249‑1 (MQ=255)
                                                                                                                                                                               cTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCGGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAg                                                                      >  2:198578/1‑250 (MQ=255)
                                                                                                                                                                                tGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAgg                                                                     >  1:621832/1‑250 (MQ=255)
                                                                                                                                                                                cgCAAACCGGGGACAGTTCGTTTGACCAAAAGGGCAACGTTTTGCCAGGTCCCAGCGAAGAGGCCGGCGCCAAAAAATGTGTGGCGCTAAAAGGGTTAACCGCCCCACCGGTGCCCTGAGGGAACCCCAGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCaa                                                                       <  2:457803/248‑1 (MQ=255)
                                                                                                                                                                                   aaaCCGGGGAAGGACGGTTTGCCCCAAGGGCAACCCTTTGGCAAGTTGCCAGCTACAGGCCCGCGCGGCAAAATGGCTGGGCCCGAAAAGGGTAAACCGACCCCCCAGGTTCTTGATGAACAACCGCCCCCCCCGCCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTg                                                                   <  2:392176/250‑1 (MQ=255)
                                                                                                                                                                                      cTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCAAAGTTCCAGCGTACAGGCCAGCGCGCAAAAATGTCGGGCGCTGAAAGGGGTAAACGTACCACCAGGTCCCGGATGAAACACCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAg                                                               <  1:791710/250‑1 (MQ=255)
                                                                                                                                                                                                 aTCGTTTGCCCAAAAGGCCAACTCTTTGCAAAGTTCCAGCGTACAGGCCAGCGCGCAAAAATGTCGGGCGCTGAAAGGGGAAAACGTACCACCGGGTTCCTGATGAAACACCAGCGCCCCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGAt                                                      <  1:234528/250‑1 (MQ=255)
                                                                                                                                                                                                         gCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGCTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGCATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGcc                                              >  2:693404/1‑250 (MQ=255)
                                                                                                                                                                                                          cccAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGGTTCTTTGccc                                             >  2:13435/1‑250 (MQ=255)
                                                                                                                                                                                                            aaaaGGGCCAGCGCTTTGCAAAGTGCGCGCGCAGAGGCCGGCGCGACAAAATGTCGGGGCGAGAATGGGATAAACGTCCCACCAGGTTCCTGAGGGAACCCCCGCCCCTCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCt                                            <  2:378670/249‑1 (MQ=255)
                                                                                                                                                                                                                    ccaccGCTTGGCAAACTGCCAGCGTCCAGGCCAGCGCGCAAAATTGCGTGGCGTGAAAAGGGGAAACCCCCCCCCCAGTTTCCGGATGGACCACCGGCGCCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATa                                     <  2:341671/246‑1 (MQ=255)
                                                                                                                                                                                                                         gCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGCTCAAGGTGCCAGCGCACAGATTCTTTGCCTCGGTCATAAAAAt                                >  2:413442/1‑250 (MQ=255)
                                                                                                                                                                                                                         gCTTTGCCAAGTCCAAGCGTACAGGCCGGCGCGCAAAAATGTCGGGCGCTGAAATGGATAAACGTACCACCAGGTTCCTGAGGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAAt                                <  1:670049/250‑1 (MQ=255)
                                                                                                                                                                                                                            ttGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGAAGATCAAGGAGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCtt                             >  2:548938/1‑250 (MQ=255)
                                                                                                                                                                                                                             tGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTc                            >  1:710446/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ccAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATACACATCTTCCAGGCACTg                   >  2:67538/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ccACCGTACGGCCCACCGCGCAAAAATGTGTGGCGCTGAAATGGGAAAACTCCCCACCAGGTCCCTGAAGGAACACCAGCCCCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACt                    <  1:670893/249‑1 (MQ=255)
                                                                                                                                                                                                                                                  gcccAGCGCGCAAAATGGTCGGGCGCTGAAAGGGGAAAACGTCCCACCAGTTCCCTGAGGAAACCCCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTcg         <  1:466303/248‑1 (MQ=255)
                                                                                                                                                                                                                                                   gCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTcgc        >  2:841695/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    ccgggcgcacAAAAGGGCTGGCGCGAAAAGGGGTAAACTGCCCACCAGGTTCCTGATGGAACCCCCGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTcgcg       <  1:431970/241‑1 (MQ=255)
                                                                                                                                                                                                                                                        cgcgcCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGcc   >  2:427312/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         gcgcAAAAAGGTCTGGGCATGAAAGGGGAAACCGCACCCCCGGGTCCGGGATGAAACCCCAGCCCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGccc  <  1:490603/250‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          
GCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCC  >  NZ_CP009273/1820327‑1820819

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