Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_LN832404 4,061,440 A→G 6.1% L267P (CTG→CCG)  EK12AG100_RS20130 ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_LN8324044,061,4400AG 6.1% 219.5 / 4.0 66L267P (CTG→CCG) EK12AG100_RS20130MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (26/36);  new base G (2/2);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

CGGTTTTCCATTCACCGTAATCTACGGTATCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGC  >  NZ_LN832404/4061156‑4061726
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                              
cgtgtttCCAATCACCGTAATCTACGGTATCTTGTACCAGTGCCCAGTTCAGTCTGGTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATa                                                                                                                                                                                                                                                                                         <  1:522009/288‑1 (MQ=255)
  gTTTTCCATTCACCGTAATCTACGGTATCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACt                                                                                                                                                                                                                                                                                       <  1:336400/292‑1 (MQ=255)
     ttCCATTCAACGTAATCTACGGTATCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCCGAAATAATTCAGTATATCACCCACCGCCCCCAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGAAGAGGATGCATCCCATATTGAAGAATCCCATCAACGACAATAAATTAATATCATTCAGCACATACAggg                                                                                                                                                                                                                                                                                    >  1:129852/1‑292 (MQ=255)
         aTTCACCGTAATCTACGGTATCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTa                                                                                                                                                                                                                                                                                >  1:40267/1‑292 (MQ=255)
              ccGTAATCTACGGTATCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGa                                                                                                                                                                                                                                                                           <  1:57320/292‑1 (MQ=255)
                  aaTCTACGGTATCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCTGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGCATATCACTCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGTACCAGTAATACGCCGATGAGGATGCATCCAATACCGAAAAACTCAAATCACCACAATAAATTAATATCATTAAGCCAAAACTGGGTGTAATTGActca                                                                                                                                                                                                                                                                       >  2:789979/1‑289 (MQ=255)
                       aCGGTATCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATACATTAATATCATTAAGCACATACTGGGTGTAATAGCACTGAATCg                                                                                                                                                                                                                                                                  >  2:363342/1‑292 (MQ=255)
                            aTCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCAACAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGt                                                                                                                                                                                                                                                             >  1:213925/1‑292 (MQ=255)
                                      aGTGCCCAGTTCAGGCTGGTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTaa                                                                                                                                                                                                                                                   <  2:790160/292‑1 (MQ=255)
                                       gTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTaaa                                                                                                                                                                                                                                                  <  1:64020/292‑1 (MQ=255)
                                           ccAGTACAGGCGGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCGCCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATgcc                                                                                                                                                                                                                                              <  1:789979/292‑1 (MQ=255)
                                             aGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATgccgc                                                                                                                                                                                                                                            >  1:360616/1‑292 (MQ=255)
                                              gTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATgccgcc                                                                                                                                                                                                                                           >  2:262915/1‑292 (MQ=255)
                                              gTTCAGGCTATTTAAAAAGGCCGGGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATgccgcc                                                                                                                                                                                                                                           <  2:652250/292‑1 (MQ=255)
                                                              aaGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCTCATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGtt                                                                                                                                                                                                                             <  1:765136/290‑1 (MQ=255)
                                                              aaGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGtt                                                                                                                                                                                                                             >  2:765136/1‑290 (MQ=255)
                                                                   cGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATg                                                                                                                                                                                                                      <  1:358685/292‑1 (MQ=255)
                                                                             aaGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCATAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCa                                                                                                                                                                                                            <  2:507734/292‑1 (MQ=255)
                                                                              acgCGACACAAGCGACCATGACCAAAATGACAGAGTTACATCCCCAGAAATTATTCCGTATATCACCCACCCACCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAAAAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAAGGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAg                                                                                                                                                                                                           <  2:129852/290‑1 (MQ=255)
                                                                               ggCGAAACAAGAGAACATGACCAAAGTGAAAGAGTCACTTCCCCAGAAAAAATCCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGc                                                                                                                                                                                                          <  2:40267/292‑1 (MQ=255)
                                                                                gCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCa                                                                                                                                                                                                         <  1:502704/292‑1 (MQ=255)
                                                                                            aaCAGGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTcc                                                                                                                                                                                             <  1:617937/292‑1 (MQ=255)
                                                                                              cATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGa                                                                                                                                                                                           <  1:484871/292‑1 (MQ=255)
                                                                                               aTGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCGTACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCa                                                                                                                                                                                                                                                               <  2:182673/223‑1 (MQ=255)
                                                                                               aTGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCGTACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCa                                                                                                                                                                                                                                                               >  1:182673/1‑223 (MQ=255)
                                                                                                  aCGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAata                                                                                                                                                                                       >  1:106104/1‑292 (MQ=255)
                                                                                                   cGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAatat                                                                                                                                                                                      <  1:363342/292‑1 (MQ=255)
                                                                                                            aaGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCaaaa                                                                                                                                                                             >  2:108613/1‑292 (MQ=255)
                                                                                                                             aaaTAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTg                                                                                                                                                                                                                                                          >  2:649254/1‑198 (MQ=255)
                                                                                                                             aaaTAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTg                                                                                                                                                                                                                                                          <  1:649254/198‑1 (MQ=255)
                                                                                                                               ataatTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCt                                                                                                                                                          >  1:388399/1‑292 (MQ=255)
                                                                                                                                   ttaaGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGAACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTAt                                                                                                                                                      <  2:72202/289‑1 (MQ=255)
                                                                                                                                              acccaccGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATgccgc                                                                                                                                                                                                                                            <  2:28132/195‑1 (MQ=255)
                                                                                                                                              acccaccGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATgccgc                                                                                                                                                                                                                                            >  1:28132/1‑195 (MQ=255)
                                                                                                                                              acccaccGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGt                                                                                                                                                                                                                              >  2:673581/1‑209 (MQ=255)
                                                                                                                                              acccaccGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGt                                                                                                                                                                                                                              <  1:673581/209‑1 (MQ=255)
                                                                                                                                                ccaccGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATg                                                                                                                                                                                                                      <  2:366610/215‑1 (MQ=255)
                                                                                                                                                ccaccGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATg                                                                                                                                                                                                                      >  1:366610/1‑215 (MQ=255)
                                                                                                                                                     gCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCGGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGCAATATCGCGCAAAATGATTTGAGAATACTGGCTTTCTGTCCGGTTTGCGTtgttt                                                                                                                                    >  2:662172/1‑290 (MQ=255)
                                                                                                                                                              aCCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAAt                                                                                                                           <  2:388399/292‑1 (MQ=255)
                                                                                                                                                                     cagccagATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCttt                                                                                                                    <  2:360616/292‑1 (MQ=255)
                                                                                                                                                                       gccagATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGtt                                                                                                                                              <  2:84800/264‑1 (MQ=255)
                                                                                                                                                                       gccagATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGtt                                                                                                                                              >  1:84800/1‑264 (MQ=255)
                                                                                                                                                                           gATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTAGCGTGGTGTCCCCATACTGCCCTTTGCCGct                                                                                                              >  2:241778/1‑291 (MQ=255)
                                                                                                                                                                                 ccTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTa                                                                                                        >  2:142063/1‑292 (MQ=255)
                                                                                                                                                                                  cTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATcgc                                                                                                                                                                                   <  2:156692/216‑1 (MQ=255)
                                                                                                                                                                                  cTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATcgc                                                                                                                                                                                   >  1:156692/1‑216 (MQ=255)
                                                                                                                                                                                        tttACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCtt                                                                                                                                                         >  2:126381/1‑236 (MQ=255)
                                                                                                                                                                                        tttACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCtt                                                                                                                                                         <  1:126381/236‑1 (MQ=255)
                                                                                                                                                                                             cAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGcatcat                                                                                            <  1:142063/292‑1 (MQ=255)
                                                                                                                                                                                                  cATTTTACAGTTAATGGCACCAGGAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCAATCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATaaaaa                                                                                       <  1:419863/292‑1 (MQ=255)
                                                                                                                                                                                                      gatcggTTAATGTAACCAGTAATACGCCGATAAGGATGAATCCCATACTGAAGACACCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                                   <  2:16889/287‑1 (MQ=255)
                                                                                                                                                                                                        aCAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgc                                                                                 <  1:284116/292‑1 (MQ=255)
                                                                                                                                                                                                           gTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaaagatc                                                                              >  1:16889/1‑287 (MQ=255)
                                                                                                                                                                                                           gTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgc                                                                                 >  2:7500/1‑289 (MQ=255)
                                                                                                                                                                                                           gTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgc                                                                                 <  1:7500/289‑1 (MQ=255)
                                                                                                                                                                                                               aTGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCGGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTcc                                                                                                                                  <  1:676788/236‑1 (MQ=255)
                                                                                                                                                                                                               aTGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCGGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTcc                                                                                                                                  >  2:676788/1‑236 (MQ=255)
                                                                                                                                                                                                                    aCCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCGGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGgaga                                                                     <  1:662172/292‑1 (MQ=255)
                                                                                                                                                                                                                                  gATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTcc                                                                                                                                                  >  2:265436/1‑201 (MQ=255)
                                                                                                                                                                                                                                  gATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTcc                                                                                                                                                  <  1:265436/201‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAg                       <  1:108613/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGc                <  1:262915/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aTTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGa          >  1:587347/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                                  aaTATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                                   >  2:198294/1‑216 (MQ=255)
                                                                                                                                                                                                                                                                                  aaTATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                                   <  1:198294/216‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     ttCATTCAGCACATACTGGGGTTAATATGCTTGAATGGCCGGTTTTATATTAAATGCCGCCAGGGTAAACAAGTTAGCAAAGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTTGCTTTATGTCCGGTTGGCTTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAAc    <  2:665479/291‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      tCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGTGTGAGAAAAAAGCGACTGCAACGGAATAAAca   >  1:665479/1‑291 (MQ=255)
                                                                                                                                                                                                                                                                                       cATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCTACATAATGCTCTTTGACGCCTCTGTAGCCCAGCATCATAAAAATAAAGCCGCCAATGTAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGc  >  1:373880/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                              
CGGTTTTCCATTCACCGTAATCTACGGTATCTGGTACCAGTGCCCAGTTCAGGCTGTTTACAAAGGCCGTGCCAAAAAAGGCGACACAAGAGAACATGACGAAAGTGAAAGAGTTACTTCCCCAGAAATAATTCAGTATATCACCCACCGCCCACAGCACCATGCCAGCCAGATAAACCTGTTTTTTACCAAAACATTTTACAGTTAATGGCACCAGTAATACGCCGATGAGGATGCATCCCATACTGAAGAATCCCATCCACGACAATAAATTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGC  >  NZ_LN832404/4061156‑4061726

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