Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_LN832404 693,411 T→C 100% intergenic (+194/+72) ubiF → / ← EK12AG100_RS03430 2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol hydroxylase/tRNA‑Gln

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_LN832404693,4110TC100.0% 185.3 / ‑6.7 54intergenic (+194/+72)ubiF/EK12AG100_RS034302‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol hydroxylase/tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (24/29);  minor base G (0/1);  total (24/30)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGC  >  NZ_LN832404/693127‑693700
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                 
acTGCGTTTTAGGCGTAATTCCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                                                                                                                                                                                                                                                         <  1:49172/295‑1 (MQ=255)
    cGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGCCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGGTTTTTTTATTAGGATGGGGGCAGAGGATTCGAACCTCGGAAGGCCGGATTCAGACTCCGGTGCCt                                                                                                                                                                                                                                                                                     >  2:452904/1‑295 (MQ=255)
      ttttATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTa                                                                                                                                                                                                                                                                                   >  1:231025/1‑295 (MQ=255)
       tttATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAc                                                                                                                                                                                                                                                                                  <  1:395027/295‑1 (MQ=255)
          tggggTAATCTGGGGGTAAGGGGGGGGTAGCGTGCTGGCGTGTTGAAACGTCAGTCGCGTAAATATTCGTTAGGGTGGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGc                                                                                                                                                                                                                                                                               <  2:570824/291‑1 (MQ=255)
              gTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGCATGCCGGAATCAGAACCCGGTGCCTTACCGCTCgg                                                                                                                                                                                                                                                                           >  2:120421/1‑295 (MQ=255)
                  tctcGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGAt                                                                                                                                                                                                                                                                       >  1:23703/1‑295 (MQ=255)
                   cgcGGGTTAATGGCGGCGGATCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATa                                                                                                                                                                                                                                                                      <  1:663839/293‑1 (MQ=255)
                     tgggTGAAGGGCGGCGGGAGGTGGTGGCGAGGTAAAACGGAAGGCGCGGGAGTCTGAGTTAGGGTTGTAGGCTGACAACATTGGCGGATGCGTGCAAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATGCGACATTGAAGGATAAGACGTGTCAACATCGCATTGGACATTGAATGAGCGGAGAAAAGCAAAAAGTTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATAcc                                                                                                                                                                                                                                                                    <  2:527939/294‑1 (MQ=37)
                     cGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCACGGCATGCCGGGAGCAGAACCCGGGGCCTTACCggcgggcgaaccc                                                                                                                                                                                                                                                                    >  2:133369/1‑283 (MQ=255)
                       ggTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGCATGCCGGAATCAGAATCCGTTGCCTTACCGCTTGGCTATAcccc                                                                                                                                                                                                                                                                  >  1:67049/1‑295 (MQ=255)
                        gTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCa                                                                                                                                                                                                                                                                 >  1:671208/1‑295 (MQ=255)
                            aTGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACCTTGAATGAACGCAGAAAAGCCAAAAGATCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCTGAATCCGGTGCCTTACCGCTTGGCGATACCCCACCTg                                                                                                                                                                                                                                                             >  1:657512/1‑295 (MQ=255)
                                      aGCGTGCGGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTa                                                                                                                                                                                                                                                   <  1:120421/295‑1 (MQ=255)
                                                 gtgtTGCAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTAGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAACGACCGCAGAAAAGCAACAAGCTCGCCGAAGCGAGCTTTTTTCATATGGCTGCGGTACTAGGATTCGAACCTCGGAATGCCGGAAGCAGCATCTGGTCCCTTACCGCTCGGCGCTTCCACAACTCGGTACTCTTACCAgtacagca                                                                                                                                                                                                                                        >  2:794413/1‑288 (MQ=255)
                                                   gtTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATACGACATTGAAGGATAAGACGTGTCAACATAGCATTAGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTCATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATACGGTGACTTACCGCTAGGCGATACCACAACAGAGTGGAGTTACAAGGGAAGcaga                                                                                                                                                                                                                                      >  1:475625/1‑292 (MQ=38)
                                                          ggTCAGGGGCGGAAAGTGGCGTTAGGGTTGTATCCGTACAACAGTGCAGGATGCGTGACAACCGTAGGTCGGATAAGACGCGGCAGCGTCACATCCGACAATGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTGAAGGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATa                                                                                                                                                                                                                               <  2:724290/294‑1 (MQ=38)
                                                          cGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACACGGTAAGCGTCTTGAATa                                                                                                                                                                                                                               >  1:463904/1‑295 (MQ=255)
                                                            tCAGGCGATGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCATCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATaaa                                                                                                                                                                                                                             <  1:490234/295‑1 (MQ=255)
                                                              aGGCGCTGAAATCTGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTATGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGCATCCGGTGCCTTACCGCTTGGCCATACCCCAACTGCGTGAGCTTACAATGTAAGCGCCTTGACTaaact                                                                                                                                                                                                                           >  2:24456/1‑293 (MQ=39)
                                                               ggCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACCGGTTGCACTTCCAAGGTCCGCGCCTTGCATAAATTg                                                                                                                                                                                                                          >  2:140923/1‑295 (MQ=255)
                                                                gcgcTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTgg                                                                                                                                                                                                                         >  1:180782/1‑295 (MQ=255)
                                                                gcgcTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAGGGTGCACTTGCAAGGTACGCGTCTAGAAGAAATggg                                                                                                                                                                                                                         >  1:295374/1‑295 (MQ=255)
                                                                          atGCTTTAGGGTTCTAGCCTGACAACATGCGCCGATGCGTGCCAACCGTAGGTCGGATAATACGCTGCAGCGTCGCATCCGACATTGAACGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACg                                                                                                                                                                                                               <  2:476328/295‑1 (MQ=39)
                                                                               ttAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTTGCTGGGGTACgatgat                                                                                                                                                                                                          >  1:625747/1‑295 (MQ=255)
                                                                                     ctGTAGCCATACAACATTGCCGGGTGCGTGCCAACCGTAGGGCGGATAAGACGCGGGAGCGGCGCATCCGACATTGAAGGATAAGAGGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAAc                                                                                                                                                                                                    <  2:595642/294‑1 (MQ=255)
                                                                                                cgagatggcgggAGGCGTGCCGAGCGGAGGTGTGATAAGACGCGGCAGCGGCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGGATTCGACAGTGAATGAACGCAGAAAAGCAAAAAGGTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGcc                                                                                                                                                                                         <  2:241908/285‑1 (MQ=38)
                                                                                                      tGCCGGATGCGTGCCAACCGTAGGTCGGAGAAGACGCGTCAGCGTCGCATCCGACATGGAAGGATAAGACGTGTCAACAGCGCATTCGACATTGAATGAACGCAGAAAAGAAAAAAGCTCGCCGAAGCGAGCTTGTTTAATGTGGCTGGGGTACGAGGATTCTAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATc                                                                                                                                                                                   <  1:475286/295‑1 (MQ=255)
                                                                                                           gATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTACGCGTCTTGAATAAATCTGCCGGGTTACGAGGCTTCGAACGTCGGAATGCCGGCATCAGCAt                                                                                                                                                                              >  2:315747/1‑295 (MQ=39)
                                                                                                           gATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAt                                                                                                                                                                              >  1:578542/1‑295 (MQ=255)
                                                                                                                gTGCCAAACGTAGGTAGGGTAAGACTCAGCAGAGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCTCATTCGACATTGAAGGAACGCAGAAAAGCAAAAAGCTCGCCGGAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                                                                                                                                         <  1:671328/295‑1 (MQ=39)
                                                                                                                   ccgacgtgaGGACGGGTGAGACGCGGGAGCGTCGCATACGACAGTGAAAGATAAGAGATGTCTGCATGGAAATGAACATTGAAAGAACGCAGAAAAGCAAAAAGCTGGCCGAAGCGAGCTTTGTTAATGGGGATGGGGTACGAGGATTCAAAGCACAGAATGCCAGAAAGAGAATCCGGTGCTTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGGATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                                                                                                                                      <  2:747420/287‑1 (MQ=255)
                                                                                                                     aaCCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTCCACGGTAAGCGTCTAGAAGAAATTGGCTGGGGTACGAGGAATGGAAGCTCGGAAGGCCGGCATCAGAATCCCGTGCCtt                                                                                                                                                                    >  2:327617/1‑295 (MQ=39)
                                                                                                                      aCCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTATTGAATAAATTGGAAGGGGTACGAGGCTTCGAACGTCGGAATGCCTGGACTAGAAACCGGTCGCTTa                                                                                                                                                                   >  2:709160/1‑295 (MQ=38)
                                                                                                                                     aGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGTTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATTCCGGAATCAGAATCCGTTGCCTTACCGCTTCGCGATAcc                                                                                                                                                    >  1:4670/1‑295 (MQ=255)
                                                                                                                                      gACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCCCTTACAAGGTAAGCGTCTTGAATATATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGTCTTACCGCTTGGCGATAccc                                                                                                                                                   >  1:515308/1‑295 (MQ=255)
                                                                                                                                           ggCAGCGTCGCATCCGACATTGGAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaaca                                                                                                                                              <  1:596760/295‑1 (MQ=255)
                                                                                                                                                 gTCGCATCCGACATTGAAGGATAGGACGTGTCAACTTTGCAATCGACATTGAATGAACGCAGAAAAGCAAAAAGCACGCTGAACCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTgg                                                                                                                                        <  1:751635/295‑1 (MQ=255)
                                                                                                                                                  tCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCCCGGAATGCCGGCATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGCGCACTTACACGGTAAGCGTCTTGACTAAATTGGCTGGGGTACGCGGCTTCGAAGCTCGTAATGCCGGAATAAGACTCCTGTGCCCTACCGTATTGCGATACACCAACCAATTgga                                                                                                                                       >  2:590933/1‑294 (MQ=38)
                                                                                                                                                   cGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGAGTCTTGAATAAATTGGCTGGTGTACGAGTCTTCGGAGCTCGGCATGCCGGGATCAGCATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAtcgctt                                                                                                                                      >  1:255384/1‑290 (MQ=38)
                                                                                                                                                    ggaaccGACATGGAAGGGTAAGACGTGGCAACTGCGCATTCGGCAATGAAGGAACGCAGAAAAGCAAAAAGCTGGGCGGAGCGAGCTGTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAGCGCTTGGAGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGttt                                                                                                                                     <  2:651317/291‑1 (MQ=38)
                                                                                                                                                    gCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGCACCTCGGCATGCCGGAATCAGAATCCGGTGCCTTCCAGCTTGGCGATACCCCAACCAATTTGttt                                                                                                                                     >  1:755221/1‑295 (MQ=255)
                                                                                                                                                                 aaGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGaa                                                                                                                        <  1:327617/295‑1 (MQ=255)
                                                                                                                                                                                                    aatgaaCGCGGAAAAGAAAAAAGGTCGCCGAAGCGAGCGTTGTTAATGGGGGTGGGGTACGGGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCGACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTacgac                                                                                     <  2:763281/295‑1 (MQ=39)
                                                                                                                                                                                                    aaagaacaacgAAAAAGAAAACGCTCGCCGAAGCGAGCGTTTTTAGTGTGGCTGGGGGACGGGGAGTCGAACCTCGGAATGCCGGAAGCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTacgac                                                                                     <  2:566956/285‑1 (MQ=38)
                                                                                                                                                                                                                  aGCAAAAAGCGCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCtg                                                                       <  1:333989/295‑1 (MQ=255)
                                                                                                                                                                                                                   gCAAAAAGCTCGCCGAAGCGAGCTTGTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCtgt                                                                      <  1:796581/295‑1 (MQ=255)
                                                                                                                                                                                                                     aaaaaGCACGCCGGAGCGAGCTTTGTTAATGTGGCTGTGGTACGAGGATTCGAGCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAAGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGa                                                                    <  1:140923/295‑1 (MQ=255)
                                                                                                                                                                                                                      aaaaGCGCGCCGAAGCTGGCTTTTTTAATGTGGGTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGAc                                                                   <  1:545756/295‑1 (MQ=255)
                                                                                                                                                                                                                          gCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCca                                                               >  1:633051/1‑295 (MQ=255)
                                                                                                                                                                                                                           cTCGCCGAAGCGAGCTTGTTTAATGAGGCTGGAGTACGAGGATTCGAACCTCGGAAAGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCcat                                                              <  1:452904/295‑1 (MQ=255)
                                                                                                                                                                                                                             cgccgcAACGAGGTGTGTGAATGTGAATGGGGTACGGAGATTCGAACCTCGGAATGACAGAAGCAGAATACGGTGCCTTACCGCTTGGCGATACGCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCcatca                                                            <  2:220391/295‑1 (MQ=38)
                                                                                                                                                                                                                                       gagtgtttgaggggggCGGGGGGAGGGGGATTGGCACGTTGGAATGCCGGAATAAGAATCCGGTGCCTTACCGCTTGGCGAGACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGa                                                  <  2:724926/281‑1 (MQ=38)
                                                                                                                                                                                                                                                 aTGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTaa                                        <  1:635589/295‑1 (MQ=255)
                                                                                                                                                                                                                                                        tccgatctgAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACGAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTATGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTg                                 <  1:46901/287‑1 (MQ=255)
                                                                                                                                                                                                                                                              aCGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTAc                           <  1:315747/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ccTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCGACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGt              <  1:697243/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ccTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGt              >  1:531074/1‑295 (MQ=255)
                                                                                                                                                                                                                                                                             ttggAATGACGAAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAAGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGttt            <  1:709160/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 aaTGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTttcttc        <  1:590933/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ggAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGTACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGc  >  1:886361/1‑295 (MQ=255)
                                                                                                                                                                                                                                                                                       agAATCAGAATGCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGAGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGc  <  1:296214/294‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                 
ACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGC  >  NZ_LN832404/693127‑693700

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