Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,891,283 A→G D32G (GAT→GGT)  sdaA → L‑serine ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,891,2830AG95.7% 303.8 / 3.7 92D32G (GAT→GGT) sdaAL‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (49/39);  total (51/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

GTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGC  >  NZ_CP009273/1891151‑1891415
                                                                                                                                    |                                                                                                                                    
gTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGt                                                                                                                                >  2:1192027/1‑139 (MQ=255)
          gTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGc                                                                                                                      >  1:1265804/1‑139 (MQ=255)
              aGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTAc                                                                                                                  >  2:1016540/1‑139 (MQ=255)
                            gAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACt                                                                                                    >  1:689344/1‑139 (MQ=255)
                                 ttatCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGc                                                                                                                      >  2:1149450/1‑116 (MQ=255)
                                 ttatCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGc                                                                                                                      <  1:1149450/116‑1 (MQ=255)
                                 ttatCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTcgcg                                                                                                <  2:689344/138‑1 (MQ=255)
                                  tatCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGtt                                                                                              >  1:308976/1‑139 (MQ=255)
                                      gTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCg                                                                                          >  1:1057102/1‑139 (MQ=255)
                                       tGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCATCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGt                                                                                         >  1:1226865/1‑139 (MQ=255)
                                         aTTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTg                                                                                        >  2:606071/1‑138 (MQ=255)
                                         aTTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTg                                                                                        >  2:673016/1‑138 (MQ=255)
                                             gTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGt                                                                                   <  2:985238/139‑1 (MQ=255)
                                              tCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGtt                                                                                  >  1:862164/1‑139 (MQ=255)
                                                tATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTg                                                                                        <  1:505276/131‑1 (MQ=255)
                                                tATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTg                                                                                        >  2:505276/1‑131 (MQ=255)
                                                  ttCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCtt                                                                                                                    <  1:854124/101‑1 (MQ=255)
                                                  ttCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCtt                                                                                                                    >  2:854124/1‑101 (MQ=255)
                                                  ttCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTa                                                                                <  1:539259/137‑1 (MQ=255)
                                                  ttCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTa                                                                                <  1:673016/137‑1 (MQ=255)
                                                   tCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATgg                                                                             >  1:107872/1‑139 (MQ=255)
                                                        aTGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTg                                                                                        >  2:1043349/1‑123 (MQ=255)
                                                        aTGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTg                                                                                        <  1:1043349/123‑1 (MQ=255)
                                                          gTTTAAGGTGGGGATTGGTCCCTCCTCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTg                                                                                        >  2:747728/1‑121 (MQ=255)
                                                          gTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTg                                                                                        <  1:747728/121‑1 (MQ=255)
                                                              aaGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCg                                                                   <  1:782922/138‑1 (MQ=255)
                                                              aaGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCg                                                                   <  1:1192027/138‑1 (MQ=255)
                                                                ggTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTg                                                                <  2:98147/139‑1 (MQ=255)
                                                                ggTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTg                                                                <  1:1149223/139‑1 (MQ=255)
                                                                    gggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACgg                                                            >  1:1009357/1‑139 (MQ=255)
                                                                    gggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACgg                                                            >  2:14263/1‑139 (MQ=255)
                                                                     ggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTcgcg                                                                                                <  2:1064249/102‑1 (MQ=255)
                                                                     ggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTcgcg                                                                                                >  1:1064249/1‑102 (MQ=255)
                                                                     ggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACggg                                                           >  2:670988/1‑139 (MQ=255)
                                                                     ggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACggg                                                           >  2:1189028/1‑139 (MQ=255)
                                                                     ggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACggg                                                           >  1:293924/1‑139 (MQ=255)
                                                                       aTTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTa                                                         <  2:293924/139‑1 (MQ=255)
                                                                           gTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGa                                                                                                              >  1:701961/1‑82 (MQ=255)
                                                                           gTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGa                                                                                                              <  2:701961/82‑1 (MQ=255)
                                                                           gTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAgg                                                     >  1:539023/1‑139 (MQ=255)
                                                                               cTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                 <  1:1016540/139‑1 (MQ=255)
                                                                                tCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGccac                                                 >  2:480915/1‑138 (MQ=255)
                                                                                 cATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGa                                                                                                              >  1:684558/1‑76 (MQ=255)
                                                                                 cATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGa                                                                                                              <  2:684558/76‑1 (MQ=255)
                                                                                       cccATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTa                                                                                                                   >  2:476728/1‑65 (MQ=255)
                                                                                       cccATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTa                                                                                                                   <  1:476728/65‑1 (MQ=255)
                                                                                       cccATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCAc                                                                        <  1:511417/108‑1 (MQ=255)
                                                                                       cccATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCAc                                                                        >  2:511417/1‑108 (MQ=255)
                                                                                       cccATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACcaca                                             <  2:481412/135‑1 (MQ=255)
                                                                                       cccATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACcaca                                             <  1:898708/135‑1 (MQ=255)
                                                                                       cccATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACcaca                                             <  1:606071/135‑1 (MQ=255)
                                                                                          aTACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGatat                                      <  2:107872/139‑1 (MQ=255)
                                                                                          aTACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGatat                                      <  1:644492/139‑1 (MQ=255)
                                                                                              cGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGcc                                  <  2:274976/139‑1 (MQ=255)
                                                                                               gTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCtt                                                                                                                    >  1:287842/1‑56 (MQ=255)
                                                                                               gTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCtt                                                                                                                    <  2:287842/56‑1 (MQ=255)
                                                                                                tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGcc                                                                                           >  1:449565/1‑80 (MQ=255)
                                                                                                tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGcc                                                                                           <  2:449565/80‑1 (MQ=255)
                                                                                                tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTa                                                                                >  1:435565/1‑91 (MQ=255)
                                                                                                tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTa                                                                                <  2:435565/91‑1 (MQ=255)
                                                                                                tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                >  1:1100444/1‑139 (MQ=255)
                                                                                                tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                >  2:991508/1‑139 (MQ=255)
                                                                                                  gggCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCatta                              >  2:1041333/1‑139 (MQ=255)
                                                                                                    gCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGttt                                                                                 <  2:802751/86‑1 (MQ=255)
                                                                                                    gCCTATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGttt                                                                                 >  1:802751/1‑86 (MQ=255)
                                                                                                      cTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATtattat                          <  1:1042488/139‑1 (MQ=255)
                                                                                                       tATGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATg                         >  2:299515/1‑139 (MQ=255)
                                                                                                        aTGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTg                                                                >  1:691318/1‑99 (MQ=255)
                                                                                                        aTGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTg                                                                <  2:691318/99‑1 (MQ=255)
                                                                                                        aTGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATgg                        >  2:868107/1‑139 (MQ=255)
                                                                                                         tGAAGGCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg                       >  2:389763/1‑139 (MQ=255)
                                                                                                              gCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAg                                                                                                                        <  2:403440/37‑1 (MQ=255)
                                                                                                              gCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAg                                                                                                                        >  1:403440/1‑37 (MQ=255)
                                                                                                              gCAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                  >  1:531441/1‑139 (MQ=255)
                                                                                                               cAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg                       <  1:685054/133‑1 (MQ=255)
                                                                                                               cAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATggg                       >  2:685054/1‑133 (MQ=255)
                                                                                                               cAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                 <  1:480915/139‑1 (MQ=255)
                                                                                                               cAGGTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                 <  2:295733/139‑1 (MQ=255)
                                                                                                                 ggTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGAc                                                              <  2:611541/92‑1 (MQ=255)
                                                                                                                 ggTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGAc                                                              >  1:611541/1‑92 (MQ=255)
                                                                                                                 ggTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg               >  2:366826/1‑139 (MQ=255)
                                                                                                                  gTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCa                                 <  2:1184180/120‑1 (MQ=255)
                                                                                                                  gTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCa                                 >  1:1184180/1‑120 (MQ=255)
                                                                                                                  gTAAACAGTTCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg              >  1:8188/1‑139 (MQ=255)
                                                                                                                          ttCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCg                                                                   >  1:498009/1‑78 (MQ=255)
                                                                                                                          ttCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCg                                                                   <  2:498009/78‑1 (MQ=255)
                                                                                                                          ttCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATtattat                          <  1:419833/119‑1 (MQ=255)
                                                                                                                          ttCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATtattat                          >  2:419833/1‑119 (MQ=255)
                                                                                                                          ttCGTCGATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc      >  2:829664/1‑139 (MQ=255)
                                                                                                                            cgtcgATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt    >  1:536621/1‑139 (MQ=255)
                                                                                                                             gtcgATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTg   <  1:1230491/139‑1 (MQ=255)
                                                                                                                               cgATGGTCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc  >  2:1031235/1‑138 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGC  >  NZ_CP009273/1891151‑1891415

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