Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,744,804 T→C E156G (GAG→GGG)  ydhX ← 4Fe‑4S dicluster domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,744,8040TC98.4% 220.3 / ‑3.9 63E156G (GAG→GGG) ydhX4Fe‑4S dicluster domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (30/32);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

GGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCA  >  NZ_CP009273/1744669‑1744938
                                                                                                                                       |                                                                                                                                      
ggTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAg                                                                                                                                     >  1:602343/1‑139 (MQ=255)
    gCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAgcgc                                                                                                                                  <  1:630702/138‑1 (MQ=255)
     cTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAgcg                                                                                                                                   >  2:1055519/1‑136 (MQ=255)
     cTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAgcg                                                                                                                                   >  1:704834/1‑136 (MQ=255)
     cTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAgcg                                                                                                                                   >  2:370040/1‑136 (MQ=255)
       gATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGa                                                                                                                               >  1:283460/1‑138 (MQ=255)
          atTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAAt                                                                                                                            >  1:1097733/1‑138 (MQ=255)
             gATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATcac                                                                                                                         <  1:47314/138‑1 (MQ=255)
                  tatttattTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACt                                                                                                                                        <  1:731831/118‑1 (MQ=255)
                  tatttattTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACt                                                                                                                                        >  2:731831/1‑118 (MQ=255)
                  tatttattTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAAt                                                                                                                            <  2:482872/130‑1 (MQ=255)
                  tatttattTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAAt                                                                                                                            >  1:482872/1‑130 (MQ=255)
                  tatttattTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTAt                                                                                                                   >  2:554705/1‑139 (MQ=255)
                  tatttattTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTAt                                                                                                                   >  1:119159/1‑139 (MQ=255)
                  tatttattTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTAt                                                                                                                   >  2:1140051/1‑139 (MQ=255)
                  tatttattTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTAt                                                                                                                   >  1:981943/1‑139 (MQ=255)
                        ttttGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGcc                                                                                                              <  1:370040/138‑1 (MQ=255)
                            gctgcAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACtt                                                                                                          >  2:192497/1‑138 (MQ=255)
                            gctgcAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACtt                                                                                                          >  2:633568/1‑138 (MQ=255)
                            gctgcAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACtt                                                                                                          >  2:686265/1‑138 (MQ=255)
                             ctgcAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACtt                                                                                                          >  2:438717/1‑137 (MQ=255)
                                 aaCCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGa                                                                                                    <  2:1097733/139‑1 (MQ=255)
                                  aCCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGAc                                                                                                   <  1:438717/139‑1 (MQ=255)
                                  aCCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGAc                                                                                                   <  2:1271614/139‑1 (MQ=255)
                                  aCCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGAc                                                                                                   <  2:119159/139‑1 (MQ=255)
                                       gCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGt                                                                                              >  1:26954/1‑139 (MQ=255)
                                         cTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTa                                                                                            <  1:557514/139‑1 (MQ=255)
                                                 cTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGc                                                                                     <  2:704834/138‑1 (MQ=255)
                                                 cTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTACCTGCCACTTTAGTGACAGGGTTAAGGTAGc                                                                                     <  2:1119398/138‑1 (MQ=255)
                                                    ggCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAAc                                                                                 <  1:686265/139‑1 (MQ=255)
                                                     gCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAAcc                                                                                <  1:26820/139‑1 (MQ=255)
                                                         tCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTg                                                                                                                <  2:970334/103‑1 (MQ=255)
                                                         tCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTg                                                                                                                >  1:970334/1‑103 (MQ=255)
                                                         tCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTgg                                                                             >  2:1223488/1‑138 (MQ=255)
                                                           ttCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTAc                                                                          >  1:1099080/1‑139 (MQ=255)
                                                                gACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACa                                                                     <  2:602343/139‑1 (MQ=255)
                                                                gACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACa                                                                     <  1:1140051/139‑1 (MQ=255)
                                                                  ccAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAg                                                                   <  1:1223488/139‑1 (MQ=255)
                                                                   cAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAgc                                                                  >  1:419263/1‑139 (MQ=255)
                                                                   cAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAgc                                                                  >  2:1237037/1‑139 (MQ=255)
                                                                      aGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGcc                                                               <  2:1099080/139‑1 (MQ=255)
                                                                      aGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGcc                                                               <  2:345084/139‑1 (MQ=255)
                                                                        aTCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCaa                                                             <  1:1186728/139‑1 (MQ=255)
                                                                            aTGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGc                                                          >  2:24551/1‑138 (MQ=255)
                                                                            aTGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGc                                                          >  1:585175/1‑138 (MQ=255)
                                                                                  gtTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACt                                                   >  1:473508/1‑139 (MQ=255)
                                                                                   ttCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACt                                                   <  1:554705/138‑1 (MQ=255)
                                                                                     ccGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGc                                                 <  1:24551/138‑1 (MQ=255)
                                                                                             gcgcTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGcc                                                               >  2:198460/1‑116 (MQ=255)
                                                                                             gcgcTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGcc                                                               <  1:198460/116‑1 (MQ=255)
                                                                                                  gACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGc                                   >  2:227655/1‑139 (MQ=255)
                                                                                                             gCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTAc                        >  1:691490/1‑139 (MQ=255)
                                                                                                             gCGGAAAGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTAc                        >  2:609935/1‑139 (MQ=255)
                                                                                                              cGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTAcc                       >  2:358322/1‑139 (MQ=255)
                                                                                                                  aaGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCCCTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCa                   <  1:227655/139‑1 (MQ=255)
                                                                                                                  aaGCCTTTTGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCa                   <  2:388741/139‑1 (MQ=255)
                                                                                                                       ttttGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATc              >  1:974144/1‑139 (MQ=255)
                                                                                                                         ttGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATc              <  1:315249/137‑1 (MQ=255)
                                                                                                                          tGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATccc            <  2:582549/138‑1 (MQ=255)
                                                                                                                          tGCCAGTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATccc            <  2:701463/138‑1 (MQ=255)
                                                                                                                               gTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAAt                                                       >  2:519206/1‑90 (MQ=255)
                                                                                                                               gTCGGGACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAAt                                                       <  1:519206/90‑1 (MQ=255)
                                                                                                                                  gggACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCtcat    <  2:26954/138‑1 (MQ=255)
                                                                                                                                  gggACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCtcat    <  2:974144/138‑1 (MQ=255)
                                                                                                                                   ggACCCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCtcatca  <  2:655993/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GGTAGCTGATATTGATAGTATTTATTTTGCTGCAACCACGCCTGAATTTCTGGGCTATCTTCGCGACCAAAGATCAATGCGTGTTCCGGGCAGGCGCTGACGCAAATGGGCGGAAAGCCTTTTGCCAGTCGGGACTCAGCGCAGAAATCACATTTATCTGCCACTTTAGTGACAGGGTTAAGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCA  >  NZ_CP009273/1744669‑1744938

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