Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,119,962 T→C E274G (GAA→GGA)  mdtH ← multidrug efflux MFS transporter MdtH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,119,9620TC98.7% 289.6 / ‑3.9 77E274G (GAA→GGA) mdtHmultidrug efflux MFS transporter MdtH
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (39/37);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CACGCGCAGGCTCGGCAATGATCGACCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTTCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAATTGGCAGCAT  >  NZ_CP009273/1119828‑1120095
                                                                                                                                      |                                                                                                                                     
cACGCGCAGGCTCGGCAATGATCGACCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCa                                                                                                                                    <  2:436337/138‑1 (MQ=255)
                   gATCGACCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTc                                                                                                                >  1:1173233/1‑139 (MQ=255)
                   gATCGACCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTc                                                                                                                >  2:1239626/1‑139 (MQ=255)
                       gACCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAg                                                                                                             <  2:381899/138‑1 (MQ=255)
                        aCCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCg                                                                                                           <  1:610997/139‑1 (MQ=255)
                        aCCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCg                                                                                                           <  1:908121/139‑1 (MQ=255)
                        aCCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCg                                                                                                           <  1:1217122/139‑1 (MQ=255)
                          ccGATATAAAACAGACAAGTCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCggg                                                                                                         <  2:1184990/139‑1 (MQ=255)
                            gATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCAACCGGCATCATGCTTAATGCCATTATCAACAGCCAAGCCATAAACCGGTGTCCAAGACGAAATTGCTTTTACCTCAAGCGGGc                                                                                                        >  2:729164/1‑138 (MQ=255)
                              tataAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGa                                                                                                      >  1:661670/1‑138 (MQ=255)
                                  aaaCAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAggg                                                                                                 <  2:664014/139‑1 (MQ=255)
                                  aaaCAGACAAAACAGGGGTAAAAGTTTTTGCAGGCCGCCTACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAggg                                                                                                 <  1:729164/139‑1 (MQ=255)
                                      aGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTa                                                                                               >  1:122377/1‑137 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCt                                                                                                                        >  1:13347/1‑110 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCt                                                                                                                        <  2:1049383/110‑1 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCt                                                                                                                        <  2:13347/110‑1 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCt                                                                                                                        >  1:1049383/1‑110 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGa                                                                                                      <  2:1200237/128‑1 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGa                                                                                                      >  1:39023/1‑128 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGa                                                                                                      >  1:1200237/1‑128 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGa                                                                                                      <  2:39023/128‑1 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTagag                                                                                            >  2:1015771/1‑138 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTagag                                                                                            >  2:100594/1‑138 (MQ=255)
                                        aCAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTagag                                                                                            >  1:914766/1‑138 (MQ=255)
                                         cAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCa                                                                                          <  1:974411/139‑1 (MQ=255)
                                                    gAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACga                                                                                >  2:118922/1‑138 (MQ=255)
                                                     aaaaGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgac                                                                               <  1:1128359/138‑1 (MQ=255)
                                                       aaGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacaga                                                                            <  2:661670/139‑1 (MQ=255)
                                                         gttgttGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacagaca                                                                          <  2:122377/139‑1 (MQ=255)
                                                         gttgttGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacagaca                                                                          <  2:1179612/139‑1 (MQ=255)
                                                         gttgttGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacagaca                                                                          <  1:1015771/139‑1 (MQ=255)
                                                         gttgttGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTACCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacagaca                                                                          <  2:671714/139‑1 (MQ=255)
                                                           tgttgCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACg                                                                        >  2:121136/1‑139 (MQ=255)
                                                             ttgCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGc                                                                       <  2:391744/138‑1 (MQ=255)
                                                                     cGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGc                                                              <  1:770035/139‑1 (MQ=255)
                                                                      gCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCa                                                             <  2:690308/139‑1 (MQ=255)
                                                                       cTGACCATGCCCACCGGCATCATGCTTAATGATATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCa                                                             >  1:1129441/1‑138 (MQ=255)
                                                                             aTGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATc                                                       >  1:44154/1‑138 (MQ=255)
                                                                               gCCCACCGGCATCATGCTTAATTACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacagaca                                                                          >  1:1050079/1‑117 (MQ=255)
                                                                               gCCCACCGGCATCATGCTTAATTACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacagaca                                                                          <  2:1050079/117‑1 (MQ=255)
                                                                               gCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCAt                                                    >  2:228241/1‑139 (MQ=255)
                                                                                    ccGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGc                                                                                                            >  1:1285678/1‑78 (MQ=255)
                                                                                    ccGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGc                                                                                                            <  2:1285678/78‑1 (MQ=255)
                                                                                    ccGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAAc                                               >  2:61837/1‑139 (MQ=255)
                                                                                     cGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACg                                              >  2:392302/1‑139 (MQ=255)
                                                                                        catcatGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCa                                           <  2:1173233/139‑1 (MQ=255)
                                                                                             tGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAggg                                      >  2:1249556/1‑139 (MQ=255)
                                                                                                  aaTGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGa                                                                                                      <  1:777207/70‑1 (MQ=255)
                                                                                                  aaTGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGa                                                                                                      >  2:777207/1‑70 (MQ=255)
                                                                                                  aaTGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGgcgcgc                                 >  1:48084/1‑139 (MQ=255)
                                                                                                  aaTGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGgcgcgc                                 >  2:49857/1‑139 (MQ=255)
                                                                                                   aTGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGgcgcgc                                 >  1:593302/1‑138 (MQ=255)
                                                                                                    tGACATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGgcgcgc                                 >  2:713250/1‑137 (MQ=255)
                                                                                                       cATTATCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCgccagcc                            <  1:153990/139‑1 (MQ=255)
                                                                                                            tCAACAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCAcgtc                       >  2:435678/1‑139 (MQ=255)
                                                                                                               aCAGCCCAGCCATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGtt                    <  1:435678/139‑1 (MQ=255)
                                                                                                                cagcccagccATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacagaca                                                                          <  1:985894/84‑1 (MQ=255)
                                                                                                                cagcccagccATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacagaca                                                                          >  2:985894/1‑84 (MQ=255)
                                                                                                                  gcccagccATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGAc                 <  2:914766/139‑1 (MQ=255)
                                                                                                                   cccagccATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGgcg                                    <  1:254192/119‑1 (MQ=255)
                                                                                                                   cccagccATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGgcg                                    >  2:254192/1‑119 (MQ=255)
                                                                                                                     cagccATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCa               >  1:617497/1‑138 (MQ=255)
                                                                                                                      agccATCAACCGGTGTTCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCAt              >  2:111114/1‑138 (MQ=255)
                                                                                                                        ccATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAAt           >  1:977028/1‑139 (MQ=255)
                                                                                                                        ccATCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAAt           >  1:1242376/1‑139 (MQ=255)
                                                                                                                           tCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAATTgg        <  2:617497/139‑1 (MQ=255)
                                                                                                                           tCAACCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAATTgg        <  1:118922/139‑1 (MQ=255)
                                                                                                                             aaCCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacaga                                                                            >  2:1080268/1‑69 (MQ=255)
                                                                                                                             aaCCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACgacaga                                                                            <  1:1080268/69‑1 (MQ=255)
                                                                                                                             aaCCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGa                  <  2:1143710/127‑1 (MQ=255)
                                                                                                                             aaCCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGa                  >  1:1143710/1‑127 (MQ=255)
                                                                                                                             aaCCGGTGTCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAATTGgc       <  2:593302/138‑1 (MQ=255)
                                                                                                                                  gtgtCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCa                                                     >  1:451713/1‑87 (MQ=255)
                                                                                                                                  gtgtCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCa                                                     <  2:451713/87‑1 (MQ=255)
                                                                                                                                  gtgtCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAATTGGcagcat  >  1:297024/1‑138 (MQ=255)
                                                                                                                                  gtgtCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAATTGGcagcat  >  1:1211097/1‑138 (MQ=255)
                                                                                                                                   tgtCCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTACCGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAATTGGcagcat  >  2:1003275/1‑137 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CACGCGCAGGCTCGGCAATGATCGACCCGATATAAAACAGACAAATCAGGGTGAAAAGTTGTTGCAGGCCGCTGACCATGCCCACCGGCATCATGCTTAATGACATTATCAACAGCCCAGCCATCAACCGGTGTTCCAGACGAAAATGCTTTTCACTCCAGCGGGCGATAGGGTAGAGCAACGTTAACGACAGACACGCTTCAATGGCATACATCCATTTAACGGCAGAGGGCGCGCCAGCCACGTCGTTGACCATAATTGGCAGCAT  >  NZ_CP009273/1119828‑1120095

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