Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,501,8650AC33.9% 113.8 / 13.1 59L123V (TTG→GTG) fecDFe(3+) dicitrate ABC transporter permease subunit FecD
Reads supporting (aligned to +/- strand):  ref base A (36/3);  new base C (0/20);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 6.34e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGT  >  NZ_CP009273/4501620‑4502110
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
gagaTCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGGCCTTG‑CACGCATCTTCAGCAAAACCAAaccc                                                                                                                                                                                                                                                   >  2:649863/1‑250 (MQ=255)
 agaTCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCg                                                                                                                                                                                                                                                  >  1:130349/1‑250 (MQ=255)
                cTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCa                                                                                                                                                                                                                                    >  1:160759/1‑249 (MQ=255)
                  acggCACAGGGGGAGAAATAAAAACCCCAGGGGGGTGGCAACCTGTCACAAA‑CCCCAGTCCCGGCCCCCATAGCGGGCGGGTCAGCCACAGCAGGGGGTGGTCCACACCCGGGGGGCGCGAGAGCACCAGATAACCCGTCAGG‑CTGCCCCAGCAGGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCACCGGCGGGTGGGTCTTTG‑CCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCaaa                                                                                                                                                                                                                                  <  2:48014/248‑1 (MQ=255)
                           gggCAAAAAAAAAATCACAAGCGGAATGGAAATCT‑CCAAAAGGCCCAATCCCCGC‑CCCCAAAGGCGGCCGTTCACCCAAAGCGGGGCGTTGTCCACACCCGGGGGGCGCGAGAGCATCAGATAACCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCCCCCCGCGGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCagcagc                                                                                                                                                                                                                        <  2:669167/249‑1 (MQ=255)
                                 aaaTAAAATCATCAGCGGGATGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGTCTTTG‑CCAGCATTCTCAGCAACATCAAACCCGCCAAGCCGCCCGCAAAGGCCAGCCGCCgcat                                                                                                                                                                                                                   >  2:377482/1‑248 (MQ=255)
                                  aaTAAAATCATCAGCGGGATGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCCAGCCCCCCGCAAAGGGCCGCAGCGgcagc                                                                                                                                                                                                                  >  2:509699/1‑249 (MQ=255)
                                    aaaacTCACCGGCGGGATGCCAATCT‑ACCAAAACCCCCAGCCCCGC‑CCCCATAGGCGCCCGTGCAGCCACAGCGGGGCGTTGTTCACACCCGGCGGGCGCGAGAGCATCAGATAACCCGCCAGG‑CGGGCCCAGCAGGCAAAAGGCGCCACGCCGGTGAGCGCCAGCTTCACCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAcc                                                                                                                                                                                                               <  2:197643/249‑1 (MQ=255)
                                                   gATGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc                                                                                                                                                                                                >  1:408154/1‑250 (MQ=255)
                                                     tGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGGCTTTG‑CCAGCATCTTCAGCAATATCCACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCCTCCCGGGCAGTGACgg                                                                                                                                                                                              >  2:356678/1‑250 (MQ=255)
                                                       gCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGCCTTTG‑CCAGCCTCTTCAGCAATATCAACCCCGCCCTTCCGCCCGCCACAGCCCGCAGCGGGCGCACCAACACCGGCCATGACGGcc                                                                                                                                                                                            >  2:343358/1‑249 (MQ=255)
                                                          aTCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGGGAGCGCCAGCTTCATCGGCTGGGGGGTCTTTG‑CCAGCCTCGTCCACAATATCAACCCCCCCATGCCGCCCGCAAAAGCCAGCCGCGGCACCACCATCACGAGCAGAggcggccg                                                                                                                                                                                           >  2:23811/1‑246 (MQ=255)
                                                             t‑tCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGAGGCAGCAACATCACGGGCAGTGACGGCATaag                                                                                                                                                                                        >  2:281448/1‑248 (MQ=255)
                                                                  caAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaagaag                                                                                                                                                                                     >  1:662011/1‑247 (MQ=255)
                                                                           aGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCCGGTACGGCCCAAGCAGTAAAGccccccc                                                                                                                                                                         >  2:633842/1‑248 (MQ=255)
                                                                                  gccccaaaaaGCTGCGGGCCGGCCAAGGAGGGGCGTTGTCCACACCCTGGGGGCGGGAGAGCAT‑AAAAAACCCCGCAGGCCTGCCCCAGCATGCAGAAGGCCCCACGCCGGTGAGCGCCAGCTTCACCGGCGGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGc                                                                                                                                                                   <  2:189724/242‑1 (MQ=255)
                                                                                      ccATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CGGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAgg                                                                                                                                                               <  1:428351/250‑1 (MQ=255)
                                                                                        aTAAGCTGCGGTCCACCCCCAGCGGGGCTTTTTCCCCACCCGGCGGGCGGGAGAGCATCAGATATCCCTCCGGG‑CTGGCCCAGCATGCAGAAGGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTTTTTG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCt                                                                                                                                                             <  1:618375/250‑1 (MQ=255)
                                                                                        aTAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCt                                                                                                                                                             >  1:577381/1‑250 (MQ=255)
                                                                                        aTAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGACAGGCt                                                                                                                                                             >  1:490817/1‑250 (MQ=255)
                                                                                        aTAAGCTGCCGGCAAGCCACAGCGGGGCGTTTTCCCCATCCTGCGGGCGCGAGCGCATCAGATAATCCGTCAGG‑CTGGCCCACCATGCAAAAAGCCCCACGCCGGGGAGCGCCAGCTTCATCGTCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCt                                                                                                                                                             <  2:353972/250‑1 (MQ=255)
                                                                                            gcggCCGGTCAGCCGCAGCGGGGCTTTGTCCACACCCTGCGGGCGCGAGAGAATCAGATACTCCGTCAGG‑CGGGCCAAGCATGCAGAAGGCGCCCCGCCGGTGAGCGCCAGCTTCATCGGCGGGTGGGTTTTTG‑CCAGCTTCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTggcg                                                                                                                                                         <  2:41125/247‑1 (MQ=255)
                                                                                                 cGGTCACCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATg                                                                                                                                                    <  1:20615/250‑1 (MQ=255)
                                                                                                               agGGCGTTGTTCCCACCCTGGGGGCGCGAGAGCATCAAATAACCCGTCAGG‑CGGGCCCAGCATGCAGAAAGCGCCCCGCCGGTGAGCGCCAGCTTAATCGGCTGGTGGGTCTTGG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGaa                                                                                                                                      <  1:35478/250‑1 (MQ=255)
                                                                                                                     tggtccacaCCCGGGGGGCGCGAGCGCTCAAGAAAACCCGTCGGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCGGGGGAGCGCCAGCTTCATCGGCGGGTGGTTCTTGG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCg                                                                                                                                <  2:237099/245‑1 (MQ=255)
                                                                                                                      tgttCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCTAGAATATCCgg                                                                                                                               >  1:630981/1‑250 (MQ=255)
                                                                                                                                  gCGGGCGCGGGAGCATCAGATACCCCGCAGGG‑CCGGCCCGGCATGAAGAAGGCGCCACGCCGGTGAGCGCCAGCTTCACCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAgg                                                                                                                   <  2:275754/250‑1 (MQ=255)
                                                                                                                                   cGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAggg                                                                                                                  >  1:246574/1‑250 (MQ=255)
                                                                                                                                      gcggGAGAGCATCAGATACTCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTg                                                                                                               <  1:393661/250‑1 (MQ=255)
                                                                                                                                             agCATCAGAAACCCCGCCGGG‑CTGGCCCGGAAGGAAAAAAGCGCCACGCCGGGGAGCGCCAGCTTCATCGGCGGGTGGGTCTTTG‑CCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACaa                                                                                                        <  2:183757/250‑1 (MQ=255)
                                                                                                                                                            gTCGGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCATCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCa                                                                                         <  1:392714/250‑1 (MQ=255)
                                                                                                                                                              cAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCACAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACCGCGGCCCGGCCGGCGGCATGGGTTACACCGCGAACATCCCGGGAAGCCAGAGCGCTGCGGCCAACACCCTGtccgcacacg                                                                                       >  2:149054/1‑241 (MQ=255)
                                                                                                                                                               aGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCAcgc                                                                                      >  1:554587/1‑250 (MQ=255)
                                                                                                                                                                     ggCCCAGCAGGCAGAACGCCCCCCGCGGGTAGGCCCCAGCTCCCTCGTCGGGTGGGTCTTTG‑CCAACACTTCCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCAcgcccgcc                                                                                 <  2:408154/250‑1 (MQ=255)
                                                                                                                                                                       cccAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCAAGCCGCCCGCACAGGCCAGCCGCGGCAGCACCATCACGGGCAGGGACGGCATGAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCCGCCCGGTTAACACCGCGAATATCCGGAAGAGCCCGAGGGGTGCGCACCAGCCCCGGAAACAGCActcccccccc                                                                               >  2:607442/1‑248 (MQ=255)
                                                                                                                                                                          aGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGGGATGCCAGAGGGCTGCGCACAATCCCCCGTATCAGCCCCCCCGCCCCGGc                                                                            >  2:334917/1‑250 (MQ=255)
                                                                                                                                                                             gtgaaGAACGCCCCGCGCGGGGGGGCGCAAGTTTCATGGGGGGGGGGGTTTTGG‑CCGGCATCTCCAGAAATACACCCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCCCGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAg                                                                         <  2:21966/246‑1 (MQ=255)
                                                                                                                                                                                cAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGc                                                                      >  1:397105/1‑250 (MQ=255)
                                                                                                                                                                                   aaaGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTg                                                                    >  1:442405/1‑249 (MQ=255)
                                                                                                                                                                                       cgcCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGCGATGCCAGAGGGCTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCAcc                                                                >  1:429612/1‑249 (MQ=255)
                                                                                                                                                                                          cccgcgGGGGAGCGCCACCTTCACGGGCTGGGGGGTTTTTG‑CCAGCCTCTCCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACa                                                            <  1:232547/248‑1 (MQ=255)
                                                                                                                                                                                           aCGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACa                                                            >  1:502502/1‑249 (MQ=255)
                                                                                                                                                                                                     gcgcCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCcag                                                 >  2:255219/1‑250 (MQ=255)
                                                                                                                                                                                                      ccccAGCTCCCTGGGCTGG‑GGGTTTTTGCCCAGCATCTCCACAAATACAAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCcagc                                                <  1:640509/248‑1 (MQ=255)
                                                                                                                                                                                                          aGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCagcagc                                             >  1:404176/1‑249 (MQ=255)
                                                                                                                                                                                                             ttCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAgc                                         >  1:180450/1‑250 (MQ=255)
                                                                                                                                                                                                             ttCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAAAGCCCCCACAGAAGCCAGGCTGGCGGCACGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTCTATCAGCACGCCCCCCACGGCGACGGCTGCCCCGACAAACAGTGCCACCACCCAgc                                         >  2:645280/1‑250 (MQ=255)
                                                                                                                                                                                                                     cTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGt                                 >  1:202383/1‑250 (MQ=255)
                                                                                                                                                                                                                      ggtggggTCTTGC‑CCACAATCTACAAAAATCCCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCCCCATCACGGGCAGTGACGGCAT‑AGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCg                               <  1:471174/246‑1 (MQ=255)
                                                                                                                                                                                                                            gTCTTTG‑CCAGCATCTCCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACt                          <  1:413222/250‑1 (MQ=255)
                                                                                                                                                                                                                                 tG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCAt                      >  1:587258/1‑249 (MQ=255)
                                                                                                                                                                                                                                  g‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATc                     >  2:292286/1‑249 (MQ=255)
                                                                                                                                                                                                                                    ccAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGAACTCCATCa                    >  1:588226/1‑249 (MQ=255)
                                                                                                                                                                                                                                     cAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCaa                   >  2:435841/1‑249 (MQ=255)
                                                                                                                                                                                                                                                  aaTATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCt     >  1:270144/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  aaTAACACCCCCGCCATGCCGCCCGAAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTACCACCGAGAATACCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCt     <  2:343794/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   atatCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTc    <  1:255219/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    tatTAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCg   >  1:98605/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     atCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCgt  >  1:502113/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
GAGATCGCGGCAAAAACTCAGGCTCAGCGGCAGAAATAAAATCATCAGCGGGATGGCAATCT‑TCACAAAGCTCCAGTCACGG‑CCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGG‑CTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTG‑CCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGT  >  NZ_CP009273/4501620‑4502110

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