Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,531,275 G→A L98L (CTG→TTG)  narW ← nitrate reductase molybdenum cofactor assembly chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,531,2750GA99.0% 370.8 / ‑4.4 97L98L (CTG→TTG) narWnitrate reductase molybdenum cofactor assembly chaperone
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (64/32);  total (65/32)
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.79e-01

ATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAAC  >  NZ_CP009273/1531191‑1531360
                                                                                    |                                                                                     
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  1:1187274/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  2:1878498/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  1:1277640/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  1:1075585/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  2:1677091/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  2:120229/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  2:2174973/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  1:2443501/1‑90 (MQ=255)
    ttcttTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGcc                                                                              <  2:2868517/90‑1 (MQ=255)
     tcttTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCg                                                                             <  1:2567102/90‑1 (MQ=255)
       tttCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           <  1:19936/90‑1 (MQ=255)
        ttCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAcc                                                                          <  1:2082942/90‑1 (MQ=255)
          cGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCtt                                                                        >  2:1403430/1‑90 (MQ=255)
          cGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           <  1:160986/87‑1 (MQ=255)
          cGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           >  2:160986/1‑87 (MQ=255)
          cGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCtt                                                                        >  2:1315940/1‑90 (MQ=255)
               gATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCAt                                                                   >  2:1995241/1‑90 (MQ=255)
                          aGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagc                                                        <  1:1230637/90‑1 (MQ=255)
                          aGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAACCCAACTGCAGGCCGACCTTTTCATACTCCGCcagc                                                        <  2:677200/90‑1 (MQ=255)
                           gTACGCTTAAATACTCCGGATACAGCGGGATATAATCGGGCATTGCCCGACAATCCAACAGCAGGCCGACCTTTTCATACTCCGccacca                                                       >  1:1886719/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCt                                                                         <  1:2121459/72‑1 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCt                                                                         >  2:2121459/1‑72 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTc                                                               >  2:590903/1‑82 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTc                                                               <  1:590903/82‑1 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       >  1:2220595/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       >  1:2334610/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       >  1:2439262/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       <  1:1995241/90‑1 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       >  2:2670833/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       >  2:704747/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       <  1:960465/90‑1 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       >  1:907712/1‑90 (MQ=255)
                             aCGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAgg                                                     >  2:2198517/1‑90 (MQ=255)
                                 ttAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTcca                                                 >  1:1447561/1‑90 (MQ=255)
                                  tAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTccac                                                >  2:2180035/1‑90 (MQ=255)
                                    aaTACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTccacca                                              >  2:114888/1‑90 (MQ=255)
                                     aTACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCAt                                             <  2:2736287/90‑1 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                          >  1:345758/1‑90 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                          >  2:806097/1‑90 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                          >  2:2955524/1‑90 (MQ=255)
                                         tCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGcc                                         <  2:1447561/90‑1 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           <  1:13124/55‑1 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           >  2:13124/1‑55 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                        >  1:2057328/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                        >  1:2211724/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                        >  1:2447369/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                        >  2:153527/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                        >  1:2503806/1‑90 (MQ=255)
                                            aGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTgg                                      >  2:1858857/1‑90 (MQ=255)
                                             gATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGc                                     <  1:2198908/90‑1 (MQ=255)
                                             gATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGc                                     >  2:2286553/1‑90 (MQ=255)
                                                 cAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCcgcg                                 <  1:2174973/90‑1 (MQ=255)
                                                 cAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCcgcg                                 >  1:772548/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           >  2:2063544/1‑47 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           <  1:2063544/47‑1 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                                >  2:2974906/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                                >  2:2115640/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                                >  1:314208/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                                >  2:1542899/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                                >  2:96049/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                                >  2:2148838/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                                >  2:2341178/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                                >  1:723246/1‑90 (MQ=255)
                                                      gTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATcgc                            <  2:1157711/90‑1 (MQ=255)
                                                      gTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATcgc                            >  2:74629/1‑90 (MQ=255)
                                                       tAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATcgcg                           <  1:1504962/90‑1 (MQ=255)
                                                           tAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGAc                       <  2:2211724/90‑1 (MQ=255)
                                                           tAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGAc                       >  1:35481/1‑90 (MQ=255)
                                                            aaTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACt                      <  1:74629/90‑1 (MQ=255)
                                                             aTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                     >  1:541958/1‑90 (MQ=255)
                                                              tCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCg                    <  2:381800/90‑1 (MQ=255)
                                                               cGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCgg                   >  2:893318/1‑90 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                  >  2:2800218/1‑90 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                  >  1:425662/1‑90 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                  >  1:954252/1‑90 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                  >  1:2811120/1‑90 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGCCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                  >  1:1465516/1‑90 (MQ=255)
                                                                 ggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCa                 <  2:999788/90‑1 (MQ=255)
                                                                 ggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCa                 <  2:425662/90‑1 (MQ=255)
                                                                 ggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCa                 <  2:1220635/90‑1 (MQ=255)
                                                                  gCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGAc                       >  1:2561161/1‑83 (MQ=255)
                                                                  gCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGAc                       <  2:2561161/83‑1 (MQ=255)
                                                                  gCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCAt                >  2:111585/1‑90 (MQ=255)
                                                                  gCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCAt                >  2:197738/1‑90 (MQ=255)
                                                                     gTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                     >  2:1358008/1‑82 (MQ=255)
                                                                     gTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                     <  1:1358008/82‑1 (MQ=255)
                                                                      ttCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAAc            <  2:2334610/90‑1 (MQ=255)
                                                                          cGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGt        >  1:2795776/1‑90 (MQ=255)
                                                                           gACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGtt       <  2:848012/90‑1 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATcgc                            >  1:1864712/1‑65 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATcgc                            <  2:1864712/65‑1 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                     <  2:826716/72‑1 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                     >  1:826716/1‑72 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa   >  1:2526123/1‑90 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa   >  2:1735453/1‑90 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa   >  1:1604035/1‑90 (MQ=255)
                                                                                tCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAAc  <  1:2281495/90‑1 (MQ=255)
                                                                                    |                                                                                     
ATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAAC  >  NZ_CP009273/1531191‑1531360

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: