Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,820,6360TC32.1% 97.0 / 12.0 56E182G (GAG→GGG) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (29/9);  new base C (1/17);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 8.56e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTC  >  NZ_CP009273/1820388‑1820875
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
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)
                                                                                                                                                                                                             gcTGAAAGGGGAAACCGCCCCCCCAGGTCCTGGAGGAAACCCCAGCGCCCCCCGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCtt                                         <  1:376565/250‑1 (MQ=255)
                                                                                                                                                                                                                   aTGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGc                                    >  1:644541/1‑249 (MQ=255)
                                                                                                                                                                                                                    ttggTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCa                                   >  1:135989/3‑249 (MQ=255)
                                                                                                                                                                                                                    gggggAAACGTCCCCCCCGGGTTCCTGTGGGACCCCCGCGCCCCCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCa                                   <  2:139807/249‑1 (MQ=255)
                                                                                                                                                                                                                         aaaaCCTCCCCCCAGTTTCCTGATGAACCCCCCGCGCCCCCCCCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAg                              <  2:644541/249‑1 (MQ=255)
                                                                                                                                                                                                                                    caccaGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAAc                      >  1:656360/1‑249 (MQ=255)
                                                                                                                                                                                                                                        catgTTCCGGATGAACCCCCGGCCCCTCCAGCTCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGAcc                  <  1:423545/247‑1 (MQ=255)
                                                                                                                                                                                                                                        aGGTTCCTGATGGAACACCAGCGCCTTCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCg                 >  1:73769/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCg                 >  2:231593/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTACAGCAAATAGCTGCCAACGACCGc                >  1:161016/1‑249 (MQ=255)
                                                                                                                                                                                                                                             tCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCa              >  1:318910/1‑249 (MQ=255)
                                                                                                                                                                                                                                              ccTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCaaa            >  2:402012/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    gggacccccAGCCCCCCCCGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTc       <  1:121640/243‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGTTGCCAACGACCGCCAAACATTCg      >  2:107735/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       aaCACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCAAACCATTCGAt    >  2:39467/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         cacCAGCCCCCCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTc  <  1:290067/250‑1 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
CGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTC  >  NZ_CP009273/1820388‑1820875

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