Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,264,031 A→T intergenic (‑232/‑22) uvrA ← / → ssb1 excinuclease ABC subunit UvrA/single‑stranded DNA‑binding protein SSB1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,264,0310AT97.0% 265.0 / ‑0.3 67intergenic (‑232/‑22)uvrA/ssb1excinuclease ABC subunit UvrA/single‑stranded DNA‑binding protein SSB1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (42/23);  minor base . (1/1);  total (43/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

TATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGATTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCC  >  NZ_CP009273/4263901‑4264166
                                                                                                                                   |                                                                                                                                       
tATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGtttttttt                                                                                                                                  >  1:566329/1‑138 (MQ=255)
 aTTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGtttttttt                                                                                                                                  >  1:165065/1‑137 (MQ=255)
  ttGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAg                                                                                                                               >  1:329785/1‑139 (MQ=255)
  ttGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAg                                                                                                                               <  2:676176/139‑1 (MQ=255)
  ttGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAg                                                                                                                               >  2:320714/1‑139 (MQ=255)
   tGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAg                                                                                                                               >  2:535446/1‑138 (MQ=255)
        tGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGaca                                                                                                                         >  2:99625/1‑139 (MQ=255)
          cATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGacac                                                                                                                        >  1:252403/1‑138 (MQ=255)
          cATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGacac                                                                                                                        >  1:610428/1‑138 (MQ=255)
          cATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGacac                                                                                                                        >  2:105979/1‑138 (MQ=255)
          cATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGacac                                                                                                                        >  1:594820/1‑138 (MQ=255)
          cATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGacac                                                                                                                        >  1:391434/1‑138 (MQ=255)
          cATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGacac                                                                                                                        >  1:177620/1‑138 (MQ=255)
          cATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGacac                                                                                                                        >  2:635222/1‑138 (MQ=255)
          cATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGacac                                                                                                                        >  1:574748/1‑138 (MQ=255)
                        ttgtGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCcagca                                                                                                          <  2:252403/138‑1 (MQ=255)
                             tAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCg                                                                                                    >  2:325383/1‑139 (MQ=255)
                                 aaTGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTaaa                                                                                                <  2:594820/139‑1 (MQ=255)
                                  aTGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTaaa                                                                                                >  1:775771/1‑138 (MQ=255)
                                  aTGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTaaa                                                                                                >  1:34174/1‑138 (MQ=255)
                                     tCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAg                                                                                            >  1:531021/1‑139 (MQ=255)
                                     tCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAg                                                                                            >  2:205054/1‑139 (MQ=255)
                                     tCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAg                                                                                            >  1:234499/1‑139 (MQ=255)
                                       tGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGt                                                                                          >  1:624472/1‑139 (MQ=255)
                                       tGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGt                                                                                          >  1:774124/1‑139 (MQ=255)
                                        ggCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGt                                                                                          >  2:101476/1‑138 (MQ=255)
                                         gCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTa                                                                                        <  1:101476/139‑1 (MQ=255)
                                             ggTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATtct                                                                                    <  2:574748/139‑1 (MQ=255)
                                              gtttgtttCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATtctc                                                                                   >  1:673045/1‑139 (MQ=255)
                                                        cGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATc                                                                         <  1:544729/139‑1 (MQ=255)
                                                        cGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATc                                                                         <  1:320714/139‑1 (MQ=255)
                                                               gAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCa                                                                  <  2:391434/139‑1 (MQ=255)
                                                                       aacaTAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCgg                                                          <  2:165065/139‑1 (MQ=255)
                                                                       aacaTAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCgg                                                          <  2:234499/139‑1 (MQ=255)
                                                                            aGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTa                                                     >  1:703104/1‑139 (MQ=255)
                                                                             gTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTa                                                     <  1:483109/138‑1 (MQ=255)
                                                                             gTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTa                                                     <  2:774124/138‑1 (MQ=255)
                                                                              tAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACg                                                   >  1:335753/1‑139 (MQ=255)
                                                                                aaaGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCt                                                 <  1:635222/139‑1 (MQ=255)
                                                                                 aaGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTCAGGAGACACGAAc                                                                                                                    >  1:695173/1‑70 (MQ=255)
                                                                                 aaGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTCAGGAGACACGAAc                                                                                                                    <  2:695173/70‑1 (MQ=255)
                                                                                   gcgcTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACa                                              <  2:703104/139‑1 (MQ=255)
                                                                                      cTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATg                                            >  1:600334/1‑139 (MQ=255)
                                                                                      cTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGc                                           >  1:820492/1‑139 (MQ=255)
                                                                                           ggTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCaaa                                       <  2:820492/138‑1 (MQ=255)
                                                                                              aaTGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAAtggt                                   >  1:328896/1‑139 (MQ=255)
                                                                                              aaTGGTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAAtggt                                   >  2:608486/1‑139 (MQ=255)
                                                                                               aTGGTACAATCGCGCGATTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAAtggtg                                  >  2:409043/1‑139 (MQ=255)
                                                                                                 ggTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGtt                                                                                <  1:501387/91‑1 (MQ=255)
                                                                                                 ggTACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGtt                                                                                >  2:501387/1‑91 (MQ=255)
                                                                                                   tACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGt                                                                    <  2:611128/101‑1 (MQ=255)
                                                                                                   tACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGt                                                                    >  1:611128/1‑101 (MQ=255)
                                                                                                   tACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCgg                                                          <  1:114132/111‑1 (MQ=255)
                                                                                                   tACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCgg                                                          >  2:114132/1‑111 (MQ=255)
                                                                                                   tACAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGgcg                               >  1:523674/1‑138 (MQ=255)
                                                                                                     cAATCGCGCGTTTACGCTTATTCAGAACGTTTTTTTTCAGGAGTCACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGAGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCa                             <  2:610428/138‑1 (MQ=255)
                                                                                                     cAATCGCGCGTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCa                             <  2:328896/138‑1 (MQ=255)
                                                                                                           cgcgTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATTGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTAAATGCCAAATGGTGGCGCAGTTGcc                       >  2:191864/1‑138 (MQ=255)
                                                                                                           cgcgTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGcc                       >  2:838034/1‑138 (MQ=255)
                                                                                                            gcgTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCaa                     >  1:622006/1‑139 (MQ=255)
                                                                                                             cgTTTACACTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGAGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCaa                     >  1:288412/1‑138 (MQ=255)
                                                                                                                   cacTTATTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACg              <  1:838034/139‑1 (MQ=255)
                                                                                                                        aTTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGc         >  1:744834/1‑139 (MQ=255)
                                                                                                                        aTTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGc         >  2:38181/1‑139 (MQ=255)
                                                                                                                        aTTCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGc         >  2:787990/1‑139 (MQ=255)
                                                                                                                         ttCAGAACGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCt        <  2:566329/139‑1 (MQ=255)
                                                                                                                              aaCGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATTCGGCTCAGGACCCGGAAGTACGCTACTTGCCAAATGGCGGCGCAGTCGCCAACATTACGCTGGCTACTTc   >  1:358833/1‑139 (MQ=255)
                                                                                                                              aaCGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTc   >  2:825085/1‑139 (MQ=255)
                                                                                                                               aCGTTTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTcc  <  2:107378/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAAAAGCGCTATTGGTAATGGTACAATCGCGCGTTTACACTTATTCAGAACGATTTTTTTCAGGAGACACGAACATGGCCAGCAGAGGCGTAAACAAGGTTATTCTCGTTGGTAATCTGGGTCAGGACCCGGAAGTACGCTACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCC  >  NZ_CP009273/4263901‑4264166

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