Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_LN832404 575,186 A→G 24.1% intergenic (‑134/‑83) EK12AG100_RS02845 ← / → rzoD IS5‑like element IS5 family transposase/prophage lysis lipoprotein RzoD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_LN832404575,1860AG24.1% 52.7 / 14.4 29intergenic (‑134/‑83)EK12AG100_RS02845/rzoDIS5‑like element IS5 family transposase/prophage lysis lipoprotein RzoD
Reads supporting (aligned to +/- strand):  ref base A (14/8);  new base G (4/3);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

CATGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACT  >  NZ_LN832404/575050‑575470
                                                                                                                                        |                                                                                                                                                                                                                                                                                            
ggagCAGAGTAACCTCGATTATATCCGCTCTGATTATCGGCATCATCGTCAGCCTGGTATGTGCGTGCAATCATTACCGGGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCa                                                                                                                                     <  2:627554/287‑1 (MQ=255)
       agTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCGGGTCGTCGTCGGTTGCACATCAAAGCAGTATgtc                                                                                                                              >  1:87129/1‑290 (MQ=255)
          aaCCGCGATTATTTCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTgtcagt                                                                                                                           <  1:394902/290‑1 (MQ=255)
                           ctctGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGACGGAGTTAGCTGATGCGAAAGCTGAAAATGATTCTCTGCGTGATGATGTTGCCGCTGGTCGTAGTCGGTTGCACATCACAGAAGTCTGTCAGTCAGTGCGCGAAcccccc                                                                                                          >  2:748414/1‑290 (MQ=255)
                             ctGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCcaccac                                                                                                        <  1:647314/290‑1 (MQ=255)
                                    tCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTcc                                                                                                 >  1:1042012/1‑290 (MQ=255)
                                    tCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTcc                                                                                                 >  1:829155/1‑290 (MQ=255)
                                                gcatccTTTCAAGGGCGGTCAATCATTACCTGGTTAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAAt                                                                                     <  2:412425/286‑1 (MQ=255)
                                                   gCCTGGAATGTGCGGTAAATCATTACCGAGTGAACGCAATGGCCGAAAAAGTCCAGCGCGACAAAAATGCCAGAGAACGGAAGCGAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCGGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCa                                                                                  <  2:765571/290‑1 (MQ=11)
                                                            gggCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccccc                                                                         >  1:382186/1‑290 (MQ=39)
                                                                 gttaaTCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGCAGGAGTTAGATGATGAGAAAGCTGAAAATGATACTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTTCAAATCAAAGCGGTCTGTCATTAAGTGGGTGAGGCCACCACCCCCTACGCCGTGGCTCATGCAGCATCCCCCAGACt                                                                    >  2:903995/4‑290 (MQ=1)
                                                                 gTCAATCATTACCGGGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGACAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACt                                                                    <  2:37445/290‑1 (MQ=255)
                                                                         ttACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGaca                                                            >  1:130812/1‑290 (MQ=9)
                                                                          tACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGacac                                                           >  1:853636/1‑290 (MQ=9)
                                                                              gTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACCAAGGAGTTAGCTGATGCGAAAGCTGAAACTGATGCTCTGAGTGATGATGTTTCCGCTGGTCGTCCGCGGTTGCACATCAAACCAGTCTGTCAGTCAGTTCTTGCAGCCAACCACGCCTCCGCCGTTGATAATGCCGACTCCCCCCAACTTGCacaccacgac                                                       >  2:786334/1‑284 (MQ=11)
                                                                                      gccaTCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTCCTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGa                                               <  2:358166/287‑1 (MQ=9)
                                                                                      gCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCCGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCCGAACGGGa                                               >  2:990362/1‑290 (MQ=25)
                                                                                         aTCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGAGAAGCCACCACGGCATCCGGCGTGGATAATGCAGCCTCCACCCTACTGCCATACAACGCGAACCGGGatta                                            >  2:641681/1‑290 (MQ=11)
                                                                                         aTCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGAAGTTGCCGCTGGTCGTCTGCGGTTGGACATCCAAGCCATCTGTCAGGCAGTGCGCGGAGCCACCACCGGCTCGGGCGTGGATAATGCAACCTGCCCCCGAATGGCAGACCCCGCTGGACGGGgttg                                            >  2:298233/1‑289 (MQ=11)
                                                                                            gCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTCGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACCTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGCCTGGCCGACCCCCCTGACCGGGATTAttt                                         >  2:794934/1‑290 (MQ=255)
                                                                                                  aaaagcgagggcgACAAAAATGTCAGAGAACGGAAGCGAGCGAACGCGGCAATTACTGACAGGCAGAGGCGTCAGCGTGATGTTGCGGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCt                                   <  2:588280/280‑1 (MQ=11)
                                                                                                  aaaGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAcccc                                   >  1:393745/1‑289 (MQ=18)
                                                                                                   aatgctagagtgACAAAAATGCCATGGAACTGAAGCTAGTGAAAGGGGCAATGACTGACATGAAGATGCTTCGTGGGTATGTGGCGGCGCTCGATGCAAAATACACGATGGAGTTATCGGATGCGAAAGCTGAAAATGATGCGCTGCGTTATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGGGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                  <  2:629287/279‑1 (MQ=11)
                                                                                                   aaGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                  <  1:791300/290‑1 (MQ=18)
                                                                                                    aGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCa                                 <  2:1007524/290‑1 (MQ=18)
                                                                                                    aGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGATGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGACTCCGGCGTGGCTAATGCAGCCTCCCCCCGACTGCCAGACCCCGCGTAACGGTATTATTTCACCCTCa                                 >  2:38496/1‑290 (MQ=255)
                                                                                                           cgcgACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTTCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCCCCCTCagagagag                          >  1:723219/1‑290 (MQ=18)
                                                                                                            aggACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACAGGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAATTGATGCTCGTCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAgg                         <  2:142001/288‑1 (MQ=16)
                                                                                                             tgtcAAAAATGCCAGAGACCTGAAGCTGCGGAAAGCGGAAATTACTGACATGCAGATGCGTCAGCGTGAGGGTGCTTGGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTGGCCGCGGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGc                        <  2:889400/287‑1 (MQ=11)
                                                                                                             cgACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACGGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGc                        <  1:1030003/290‑1 (MQ=21)
                                                                                                                  aaacagCCAGAGAAGTGGAGCTGGCGGAAGTGGCAGTTGCTGACACGCTAAGGGGGCCGCGTGGTGGTGCTGCGCTGCATGCCACATCCCCGAAGGCGTTATCGGATGCGAAAGCTGAAAATGATGTCCGGCGGGATGCTGTTGCCGCTGGTCGGCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATc                   <  2:165749/285‑1 (MQ=11)
                                                                                                                      ggCCAGAGAATGGAAGCGAGGGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGGTGTGTTGCGGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCGGAAAATGATGCTCTGCGTGATGATGTGGCCGCTGGTCGGCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTa               <  2:713023/289‑1 (MQ=11)
                                                                                                                           gagaACTGAAGCTGGCGAACGCGTCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaa          <  1:852102/290‑1 (MQ=25)
                                                                                                                                cTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaaca     <  1:1002966/290‑1 (MQ=37)
                                                                                                                                   aaGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGCCTGGCAGACCCCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCAGATACCTATGCacaaacaaat  >  2:886868/1‑288 (MQ=255)
                                                                                                                                   aaGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGTGAGAGGCTGATCACTATGCAAAAACAACt  >  1:37616/1‑290 (MQ=255)
                                                                                                                                        |                                                                                                                                                                                                                                                                                            
CATGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACT  >  NZ_LN832404/575050‑575470

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