Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,753,672 G→A A313V (GCC→GTC)  selA ← L‑seryl‑tRNA(Sec) selenium transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,753,6720GA97.2% 234.4 / ‑2.1 73A313V (GCC→GTC) selAL‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (34/35);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

GGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGC  >  NZ_CP009273/3753536‑3753807
                                                                                                                                        |                                                                                                                                       
ggTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAggcggacg                                                                                                                                       <  2:21758/139‑1 (MQ=255)
 gTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGa                                                                                                                                      >  1:88414/1‑139 (MQ=255)
        ggCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCAtctc                                                                                                                               >  2:1280010/1‑139 (MQ=255)
        ggCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCAtctc                                                                                                                               >  1:757346/1‑139 (MQ=255)
          cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCtt                                                                                                                             >  2:807466/1‑139 (MQ=255)
                 tCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTAcc                                                                                                                      <  2:110161/139‑1 (MQ=255)
                 tCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTAcc                                                                                                                      <  1:819653/139‑1 (MQ=255)
                    gTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACcaac                                                                                                                   <  1:1032838/139‑1 (MQ=255)
                    gTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACcaac                                                                                                                   <  1:539996/139‑1 (MQ=255)
                           ttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAAtt                                                                                                            >  1:273653/1‑139 (MQ=255)
                                         cAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGcc                                                                                               <  1:142355/138‑1 (MQ=255)
                                           gggTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGc                                                                                            <  1:1071713/139‑1 (MQ=255)
                                            ggTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGcc                                                                                           >  2:580102/1‑139 (MQ=255)
                                              tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGcc                                                                                               >  2:1148155/1‑133 (MQ=255)
                                              tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGcc                                                                                               <  1:1148155/133‑1 (MQ=255)
                                              tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCt                                                                                          >  1:1221543/1‑138 (MQ=255)
                                              tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCt                                                                                          >  1:505871/1‑138 (MQ=255)
                                                taaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTaa                                                                                        >  2:133734/1‑138 (MQ=255)
                                                taaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTaa                                                                                        >  1:998264/1‑138 (MQ=255)
                                                taaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTaa                                                                                        >  2:756044/1‑138 (MQ=255)
                                                 aaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTaac                                                                                       >  2:1131795/1‑138 (MQ=255)
                                                        gACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTc                                                                               >  1:601827/1‑139 (MQ=255)
                                                             aaCGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCgg                                                                          <  1:580102/139‑1 (MQ=255)
                                                             aaCGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCgg                                                                          <  1:201231/139‑1 (MQ=255)
                                                                 tGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGaa                                                                      <  1:1112334/139‑1 (MQ=255)
                                                                   gCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGaaa                                                                     <  2:505871/138‑1 (MQ=255)
                                                                           cgccgcGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGa                                                            >  2:93957/1‑139 (MQ=255)
                                                                              cgcgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACt                                                          >  2:293064/1‑138 (MQ=255)
                                                                                cgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTaac                                                                                       <  1:284358/107‑1 (MQ=255)
                                                                                cgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTaac                                                                                       >  2:284358/1‑107 (MQ=255)
                                                                                cgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGAc                                                       >  1:809611/1‑139 (MQ=255)
                                                                                cgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGAc                                                       >  2:681857/1‑139 (MQ=255)
                                                                                     ttCATTTTATCCGCGCGTAAGGCACGCTTCAGCGGGGGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCGGGAAAAACTCACCAGACTGAcgcccg                                                  <  2:757346/138‑1 (MQ=255)
                                                                                      tccttttACCCGCGCTTAATGCCCGCTTCAGTGGGTGGCTTTGCAGGCGTACGACCATCTTTTTTTGACAAACAAAATCTCCTTCCTGCGCCCCCCCTAACAACTTCTCGCCGGAAAAACCCTCTAAACTGAcgcccgc                                                 <  2:779895/136‑1 (MQ=255)
                                                                                              tCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATc                                          <  1:807466/138‑1 (MQ=255)
                                                                                                cgcgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAAc                                       >  1:896101/1‑139 (MQ=255)
                                                                                                 gcgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACt                                      >  2:793671/1‑139 (MQ=255)
                                                                                                 gcgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGccc                                                                                              <  1:619220/83‑1 (MQ=255)
                                                                                                 gcgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGccc                                                                                              >  2:619220/1‑83 (MQ=255)
                                                                                                  cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGcc                                                                                               >  2:975157/1‑81 (MQ=255)
                                                                                                  cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGcc                                                                                               <  1:975157/81‑1 (MQ=255)
                                                                                                  cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACt                                                                   >  1:354813/1‑109 (MQ=255)
                                                                                                  cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACt                                                                   <  2:354813/109‑1 (MQ=255)
                                                                                                  cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTc                                     >  1:34106/1‑139 (MQ=255)
                                                                                                  cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTc                                     >  2:576465/1‑139 (MQ=255)
                                                                                                     gTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTCCCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTg                                  <  1:459384/139‑1 (MQ=255)
                                                                                                      tAATCCACGCTACAGCGGGTGGCTTTGCAGGCGAACGATCTCCTCTTTTTTACCAACAAAAATCCCTGCCTGCGCCCCGCCAAACAACTTGCCGCCGGAAAACCTCACCAGCCGGACGCCCGCCGCAATCAACTCCTGc                                 <  2:34106/139‑1 (MQ=255)
                                                                                                               cTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGc                        >  1:969364/1‑139 (MQ=255)
                                                                                                                   aGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGaa                                                                      >  2:373275/1‑89 (MQ=255)
                                                                                                                   aGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGaa                                                                      <  1:373275/89‑1 (MQ=255)
                                                                                                                   aGCGGGTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCtt                    >  1:1055245/1‑139 (MQ=255)
                                                                                                                        gTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCgg                >  2:1289551/1‑138 (MQ=255)
                                                                                                                        gTGGCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCgg                >  2:712478/1‑138 (MQ=255)
                                                                                                                          ggCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa              >  1:428291/1‑138 (MQ=255)
                                                                                                                          ggCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa              >  1:1198821/1‑138 (MQ=255)
                                                                                                                          ggCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa              >  2:624683/1‑138 (MQ=255)
                                                                                                                          ggCTTTGCAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa              >  1:58428/1‑138 (MQ=255)
                                                                                                                           gCTTTGCAGGCGGCCGATCATCTCTTTTTTCCCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGCCGCCGGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCaaa            <  1:560096/139‑1 (MQ=255)
                                                                                                                                gCAGGCGGCCGTTCATCTCTTTTTAACCAACAATAATTCCTCCCTGCGGCCCGCCTAAAATCTTGCCCCGGGAGAAACTCCCCAGACTGAGGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  2:809611/139‑1 (MQ=255)
                                                                                                                                gCAGGCGGACGATCATCTCTTTTTTACCAACAATAATCCCTCCCGGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  2:430193/139‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTCCCAAAAATAATTCCTGCCTGCGGCCCGCTTAACACTTTGCCGCCGGAGAAACTCACAAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  1:160360/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  1:1174990/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  1:889709/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  1:293064/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  1:624683/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  1:1224994/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  2:58428/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  2:896101/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  1:638672/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGCCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  1:133734/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCGGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  2:580693/138‑1 (MQ=255)
                                                                                                                                 cAGGCGGACGATCATCTCTTTTTTACCAACAATAATCCCTGCCTGCGCCCCGCCTAACAACTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  2:88414/138‑1 (MQ=255)
                                                                                                                                 aaGGCGGACGATCACCTCTTTTTTCCCAACAATAATTCCGGCCGGCGCCCGGCCAAACAACTTGTCTCCGGAGAAACGCTCAAGACGGGCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa       <  2:568425/137‑1 (MQ=255)
                                                                                                                                  aGGCGGACGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACt     >  2:606481/1‑139 (MQ=255)
                                                                                                                                     cggacgATCATCTCTTTTTTCCCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGc  <  1:1280010/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGC  >  NZ_CP009273/3753536‑3753807

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