Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,516,705 T→C I159T (ATC→ACC)  mrcA → peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,516,7050TC97.1% 231.6 / ‑0.6 70I159T (ATC→ACC) mrcApeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (39/29);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

CGTCACAAGGGGCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGATCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCTGAACGAAATGGCGGTGATAGCCGGGCT  >  NZ_CP009273/3516570‑3516834
                                                                                                                                       |                                                                                                                                 
cgTCACAAGGGGCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCt                                                                                                                                <  1:856668/139‑1 (MQ=255)
cgTCACAAGGGGCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCt                                                                                                                                <  2:502445/139‑1 (MQ=255)
cgTCACAAGGGGCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCt                                                                                                                                <  2:52439/139‑1 (MQ=255)
    acaAGGGGCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAg                                                                                                                            >  2:342635/1‑139 (MQ=255)
          ggCAAGTATCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTAt                                                                                                                      <  2:244355/139‑1 (MQ=255)
          ggCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTAt                                                                                                                      <  2:240521/139‑1 (MQ=255)
          ggCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTAt                                                                                                                      <  1:687246/139‑1 (MQ=255)
           gCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATc                                                                                                                     >  2:501807/1‑139 (MQ=255)
           gCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATc                                                                                                                     >  1:669409/1‑139 (MQ=255)
                    aTTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAAc                                                                                                                >  2:229511/1‑135 (MQ=255)
                    aTTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAAc                                                                                                                >  2:147605/1‑135 (MQ=255)
                    aTTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAAc                                                                                                                >  1:792191/1‑135 (MQ=255)
                          cagcagCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTAcc                                                                                                      >  2:596081/1‑139 (MQ=255)
                          cagcagCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTAcc                                                                                                      <  1:733421/139‑1 (MQ=255)
                          cagcagCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTAcc                                                                                                      >  2:617306/1‑139 (MQ=255)
                          cagcagCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTAcc                                                                                                      >  2:495677/1‑139 (MQ=255)
                          cagcagCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTAcc                                                                                                      >  2:73546/1‑139 (MQ=255)
                           agcagcTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCt                                                                                                     >  1:65779/1‑139 (MQ=255)
                                   gCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTAc                                                                                              <  1:816323/138‑1 (MQ=255)
                                   gCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTAc                                                                                              <  1:741804/138‑1 (MQ=255)
                                           cttcttCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTAt                                                                                     <  1:17041/139‑1 (MQ=255)
                                                 ccTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGgtgt                                                                                <  1:501807/138‑1 (MQ=255)
                                                  cTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCg                                                                              >  2:167542/1‑139 (MQ=255)
                                                      gTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTAc                                                                                              <  2:186115/119‑1 (MQ=255)
                                                      gTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTAc                                                                                              >  1:186115/1‑119 (MQ=255)
                                                            aaCGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCgg                                                                     >  2:181422/1‑138 (MQ=255)
                                                             acgcacgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGca                                                                   <  1:495677/139‑1 (MQ=255)
                                                              cgcacgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAATATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGcac                                                                  <  1:393527/139‑1 (MQ=255)
                                                              cgcacgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGca                                                                   >  2:656667/1‑138 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGAt                                                                                                           <  1:843943/95‑1 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGAt                                                                                                           >  2:843943/1‑95 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGgtgt                                                                                >  2:272490/1‑122 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGgtgt                                                                                <  1:272490/122‑1 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGtgctg                                                                        <  2:524132/130‑1 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGtgctg                                                                        >  1:524132/1‑130 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAg                                                               >  1:850756/1‑139 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAg                                                               >  2:749029/1‑139 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAg                                                               >  1:89857/1‑139 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAg                                                               >  2:305024/1‑139 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAg                                                               >  2:522732/1‑139 (MQ=255)
                                                                 acgcTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAg                                                               >  1:729062/1‑139 (MQ=255)
                                                                     tGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTa                                                           >  2:158324/1‑139 (MQ=255)
                                                                      gATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTAt                                                          >  1:387030/1‑139 (MQ=255)
                                                                                aTTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGaaaaa                                                <  2:792191/139‑1 (MQ=255)
                                                                                    aGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGaa                                                                                                                 >  2:44642/1‑70 (MQ=255)
                                                                                    aGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGaa                                                                                                                 <  1:44642/70‑1 (MQ=255)
                                                                                    aGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGt                                            >  2:226634/1‑139 (MQ=255)
                                                                                     ggAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTc                                           <  1:654248/139‑1 (MQ=255)
                                                                                       aaGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGa                                         <  1:226634/139‑1 (MQ=255)
                                                                                       aaGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGa                                         <  1:676994/139‑1 (MQ=255)
                                                                                        aGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGAc                                        <  1:181422/139‑1 (MQ=255)
                                                                                         gTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGAcc                                       <  1:688930/139‑1 (MQ=255)
                                                                                         gTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGAc                                        >  2:331576/1‑138 (MQ=255)
                                                                                          tCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGa                                                                                                                  >  2:251593/1‑63 (MQ=255)
                                                                                          tCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGa                                                                                                                  <  1:251593/63‑1 (MQ=255)
                                                                                          tCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGAcc                                       >  2:531386/1‑138 (MQ=255)
                                                                                          tCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGAcc                                       >  2:63066/1‑138 (MQ=255)
                                                                                           cTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCaa                                     >  1:200996/1‑139 (MQ=255)
                                                                                                   cgATTCGCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCt                             >  2:698546/1‑139 (MQ=255)
                                                                                                    gATTCGCATTGAACAGCTGCTGACGAAAGACGAGATCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACg                               >  2:131364/1‑136 (MQ=255)
                                                                                                    gATTCGCATTGAACAGCTGCTGACGAAAGACGAGATCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACg                               <  1:131364/136‑1 (MQ=255)
                                                                                                    gATTCGCATTGAAAAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGAGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCt                             >  2:491143/1‑138 (MQ=255)
                                                                                                         gCATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCTGAACg                        <  2:89857/138‑1 (MQ=255)
                                                                                                          cATTGAACAGCTGCTGACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCTGAACGaa                      <  2:850756/139‑1 (MQ=255)
                                                                                                                    ctgctgACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCTGAACGAAATGGCGGTGa            >  1:864278/1‑139 (MQ=255)
                                                                                                                         gACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTAc                                                                                              <  1:723091/52‑1 (MQ=255)
                                                                                                                         gACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTAc                                                                                              >  2:723091/1‑52 (MQ=255)
                                                                                                                         gACGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCTGAACGAAATGGCGGTGATAGcc       >  1:271594/1‑139 (MQ=255)
                                                                                                                          aCGAAAGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCTGAACGAAATGGCGGTGATAGcc       <  1:305024/138‑1 (MQ=255)
                                                                                                                              aaGACGAGACCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCTGAACGAAATGGCGGGGAGAGACGGGcg  >  2:300323/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                                 
CGTCACAAGGGGCAAGTACCATTACCCAGCAGCTGGCGAGAAACTTCTTCCTCAGTCCAGAACGCACGCTGATGCGTAAGATTAAGGAAGTCTTCCTCGCGATTCGCATTGAACAGCTGCTGACGAAAGACGAGATCCTCGAGCTTTATCTGAACAAGATTTACCTTGGTTACCGCGCCTATGGTGTCGGTGCTGCGGCACAAGTCTATTTCGGAAAAACGGTCGACCAACTGACGCTGAACGAAATGGCGGTGATAGCCGGGCT  >  NZ_CP009273/3516570‑3516834

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