Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,005,890 C→A intergenic (+138/+83) arcC → / ← yqeB carbamate kinase/selenium‑dependent molybdenum cofactor biosynthesis protein YqeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,005,8900CA97.1% 259.2 / 0.2 69intergenic (+138/+83)arcC/yqeBcarbamate kinase/selenium‑dependent molybdenum cofactor biosynthesis protein YqeB
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (32/35);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

GGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAA  >  NZ_CP009273/3005763‑3006023
                                                                                                                               |                                                                                                                                     
gggtgggtAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCg                                                                                                                             >  2:1205102/1‑138 (MQ=255)
       tAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCGTTGTTGATGCCGGCTGCGGCTCGAGCTCCTTCTCCGTGCTACAACACATCGCAAATTCAATAGATTGCAATGCACTTATACGAGTGATAAGCGTAGAGCAt                                                                                                                     >  2:647258/1‑139 (MQ=37)
          tACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAg                                                                                                                  <  2:643749/139‑1 (MQ=255)
               tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAAt                                                                                                             >  2:417067/1‑139 (MQ=255)
               tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAAt                                                                                                             >  1:11804/1‑139 (MQ=255)
               tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAAt                                                                                                             >  1:341878/1‑139 (MQ=255)
               tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAAt                                                                                                             >  1:371144/1‑139 (MQ=255)
               tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAAt                                                                                                             >  2:763470/1‑139 (MQ=255)
               tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAAt                                                                                                             >  2:891936/1‑139 (MQ=255)
                ccGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAAtt                                                                                                            >  2:619436/1‑139 (MQ=255)
                 cGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAAtt                                                                                                            <  2:11804/138‑1 (MQ=255)
                   gCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTgg                                                                                                         <  1:527534/139‑1 (MQ=255)
                       aaGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGtt                                                                                                     <  1:1205102/139‑1 (MQ=255)
                         gATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTg                                                                                                   <  1:619436/139‑1 (MQ=255)
                              tCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                                                               <  1:1000852/138‑1 (MQ=255)
                              tCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                                                               <  2:21108/138‑1 (MQ=255)
                              tCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                                                               >  1:561088/1‑138 (MQ=255)
                              tCAAGTAAGTGCGCTTTGATGATGACGGATGCGGCGCGAGCGCCTCATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                                                               <  2:616182/138‑1 (MQ=255)
                                        gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAgg                                                                                     <  1:1163631/138‑1 (MQ=255)
                                        gcgcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAgg                                                                                     >  2:1163631/1‑138 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:371144/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:322210/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:729085/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:561088/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:66930/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  1:41828/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  1:417067/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  1:909319/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:839566/139‑1 (MQ=255)
                                           cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:593709/139‑1 (MQ=255)
                                                                cgcgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCa                                                            >  2:631948/1‑139 (MQ=255)
                                                                     gcgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAtttt                                                        >  1:651116/1‑138 (MQ=255)
                                                                                 tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:539044/101‑1 (MQ=255)
                                                                                 tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 >  1:539044/1‑101 (MQ=255)
                                                                                 tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                                                               <  2:694991/119‑1 (MQ=255)
                                                                                 tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                                                               >  1:694991/1‑119 (MQ=255)
                                                                                 tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACt                                           >  1:1241076/1‑139 (MQ=255)
                                                                                   cTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTc                                         >  1:253907/1‑139 (MQ=255)
                                                                                           cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                              >  1:343674/1‑126 (MQ=255)
                                                                                           cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                              <  2:343674/126‑1 (MQ=255)
                                                                                           cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActt                                 >  2:541173/1‑139 (MQ=255)
                                                                                           cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActt                                 >  2:964424/1‑139 (MQ=255)
                                                                                           cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActt                                 >  1:964598/1‑139 (MQ=255)
                                                                                              cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                                                               <  2:812367/74‑1 (MQ=255)
                                                                                              cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGtcat                                                                                               >  1:812367/1‑74 (MQ=255)
                                                                                                cAAATTCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCtttt                            >  2:479491/1‑139 (MQ=255)
                                                                                                    ttCAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg                        <  1:891936/139‑1 (MQ=255)
                                                                                                      cAATAGATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCt                      >  1:887186/1‑139 (MQ=255)
                                                                                                           gATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                  <  2:506369/138‑1 (MQ=255)
                                                                                                           gATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                  <  2:607643/138‑1 (MQ=255)
                                                                                                           gATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                  <  1:631948/138‑1 (MQ=255)
                                                                                                           gATTGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                  <  2:881747/138‑1 (MQ=255)
                                                                                                              tGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGcc              >  2:772734/1‑139 (MQ=255)
                                                                                                              tGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGcc              >  2:769034/1‑139 (MQ=255)
                                                                                                              tGCAATGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGcc              >  2:618762/1‑139 (MQ=255)
                                                                                                                   tGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 >  1:1055294/1‑67 (MQ=255)
                                                                                                                   tGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                 <  2:1055294/67‑1 (MQ=255)
                                                                                                                   tGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGAt         <  2:260246/139‑1 (MQ=255)
                                                                                                                   tGCACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAGCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGAt         >  1:574180/1‑139 (MQ=255)
                                                                                                                    gCAATTGTAGAACTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTCATCAGTCTCAGGCCGGGTGGTTAATACCATCCGGCCATTTTCATTTAATCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGAt         <  2:411334/138‑1 (MQ=37)
                                                                                                                     cACTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGc       <  2:362140/139‑1 (MQ=255)
                                                                                                                       cTTGTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAt     <  2:253907/139‑1 (MQ=255)
                                                                                                                          gTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGtct                                                                                         >  2:611418/1‑52 (MQ=255)
                                                                                                                          gTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGtct                                                                                         <  1:611418/52‑1 (MQ=255)
                                                                                                                          gTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTc                   <  1:817471/122‑1 (MQ=255)
                                                                                                                          gTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTc                   >  2:817471/1‑122 (MQ=255)
                                                                                                                          gTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCaa  >  1:860080/1‑139 (MQ=255)
                                                                                                                          gTAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCaa  >  2:1145991/1‑139 (MQ=255)
                                                                                                                           tAGGACTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCaa  <  1:999069/138‑1 (MQ=255)
                                                                                                                               |                                                                                                                                     
GGGTGGGTAATACCATCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAA  >  NZ_CP009273/3005763‑3006023

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: