Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,681,622 T→A 13.9% L216Q (CTG→CAG)  folM → dihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,6220TA13.9% 109.0 / 9.4 43L216Q (CTG→CAG) folMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (20/17);  new base A (2/4);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGT  >  NC_000913/1681328‑1681922
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                             
aTGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGa                                                                                                                                                                                                                                                                                                         <  1:472796/301‑1 (MQ=255)
     cTGGTAAGATTACTGCTTGGTCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcact                                                                                                                                                                                                                                                                                                    <  1:167679/301‑1 (MQ=255)
         tttcagacgtgtgctcttccgatctcgccgcCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttctt                                                                                                                                                                                                                                                                                                <  2:26240/276‑4 (MQ=255)
         aaaTATTACTGCGTGGGCACTGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttctt                                                                                                                                                                                                                                                                                                <  2:530539/301‑4 (MQ=255)
          aaGATTACTGCGTGGGCACGGTCACGCCGCCTGCGATATCAGTCACTTTACCGATTATGTGGGGGAGTGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttcttt                                                                                                                                                                                                                                                                                               <  1:407001/301‑5 (MQ=255)
                          gaCGGACACGCCGCCAGGGATATCATTCACTTTACGGAGTATGTGGGGGAGCGCGGTAGCGAAAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCa                                                                                                                                                                                                                                                                               <  1:396339/300‑1 (MQ=255)
                                  cgccgcCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTAGAtcacttcttaaatcggaaaagcgtcggggaaggg                                                                                                                                                                                                                                                                       >  1:26189/1‑273 (MQ=255)
                                        aggCGATATCATTCACTTTACCGATGATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTAGGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTc                                                                                                                                                                                                                                                                 <  1:105410/299‑1 (MQ=255)
                                            tttattatTCACTTTACCGATTGTGTGGTGGGGCGCGGTAGCGACAAACATTTTGCGTATGCTGCAAGCAAAGCGGCACTGGTTAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCAc                                                                                                                                                                                                                                                             <  1:339968/299‑1 (MQ=255)
                                              tatCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGATTTGTCACAGGACGCAGTTTCCCaaaa                                                                                                                                                                                                                                                           >  2:340431/1‑298 (MQ=255)
                                                    tctcTTTACTGATTGTGTGGTGGGGGGCGGTAGCGACAAACATTTTGCTTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGc                                                                                                                                                                                                                                                     <  1:112846/298‑1 (MQ=255)
                                                           aCCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCGAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgtcgtc                                                                                                                                                                                                                                              >  2:299075/1‑301 (MQ=255)
                                                                    gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTCATGAACTTGATGATGCCGCATCTCTACCACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGCCCTGGTCGATTACTTACTTACCAGTTGCTTTGTCCCCGGACGCAGTTTCCCCCTTTATGGCGATCGTCCTCTCCGtt                                                                                                                                                                                                                                     >  1:556552/1‑301 (MQ=255)
                                                                    gggggtgggcggGGTAGCGACAAACATATTGTGTATGCTGAAAGCAAAGCGGCACTGGATTATATGACCCGCTCGTTTGCCCGCATGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGtt                                                                                                                                                                                                                                     <  1:165435/293‑1 (MQ=255)
                                                                     tggtggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTa                                                                                                                                                                                                                                    >  2:165587/1‑301 (MQ=255)
                                                                                   aGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGTCTGTCGTCATCTGCGTTAATACAGCTTTAtac                                                                                                                                                                                                                      >  1:106498/1‑299 (MQ=255)
                                                                                       aCAAACATATTGCGTATGTTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCa                                                                                                                                                                                                                 <  2:541294/301‑1 (MQ=255)
                                                                                          aaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATATGCGTTAATGCAACTTTAACCAACaaat                                                                                                                                                                                                              >  1:514219/1‑301 (MQ=255)
                                                                                                                        cACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAGGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAaca                                                                                                                                                                                <  2:148028/301‑1 (MQ=255)
                                                                                                                               taataTGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACa                                                                                                                                                                         <  2:123707/301‑1 (MQ=255)
                                                                                                                                          gCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGTATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCAAGCATATAATCAACACCAATGCGCTTATCCCCCAACACGAAACAAAAACaaaaaa                                                                                                                                                              >  1:49440/1‑301 (MQ=255)
                                                                                                                                                  gCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATACGCTTATCCCCAAACACGCAACAGAACCAACAAAAAcggcgg                                                                                                                                                      >  1:477272/1‑301 (MQ=255)
                                                                                                                                                  gCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAACGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGACCAACAAAAGCGGCggg                                                                                                                                                     >  2:636534/1‑301 (MQ=255)
                                                                                                                                                           cTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGANAGAAGTGATCGANNNNNNNNATTACNNNNNNNNCAGTTGCTTNNNCACCGGACGCANNNNNNNACTTGATGGAGNNNNNNNTCTACNNNNNNACACATTATNCNNACATATAATACAAACCCAACCGCATATCCCACAACACACCACAAAACCAACaaaaacccacgcacaacgac                                                                                                                                             >  1:131616/1‑285 (MQ=255)
                                                                                                                                                               cACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGcatcat                                                                                                                                         >  2:236562/1‑301 (MQ=255)
                                                                                                                                                                        tgaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGCAAAGCGCATCATTCCGGTGaa                                                                                                                                >  2:517912/1‑301 (MQ=255)
                                                                                                                                                                                 tctttccctacactacgctcttccgttctAATGAACATGATGATGCCGAAATTCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAg                                                                                                                       <  1:209159/273‑1 (MQ=255)
                                                                                                                                                                                  ttCTATTGAGCCATGGCGTATCCTGGTTAATGAACATGATGATGCAGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGGGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTGCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGc                                                                                                                      <  1:236296/301‑1 (MQ=255)
                                                                                                                                                                                     tATTGCGCCTTCGCTGTTCCTTTTTAATGAACATGTTGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGGGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGac                                                                                                                   <  1:298671/301‑1 (MQ=255)
                                                                                                                                                                                              aTCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTaa                                                                                                          >  2:468237/1‑301 (MQ=255)
                                                                                                                                                                                                cGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATAACACAGCGAAACCAGGTaaaa                                                                                                        >  2:171567/1‑300 (MQ=255)
                                                                                                                                                                                                              aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCAAAgagatcggaaaagcacacgtctgaactcca                                                                                          >  2:209402/1‑275 (MQ=255)
                                                                                                                                                                                                              aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGCTGAAAACCCACAGCGACACCAAGAAAACAAAATAGGAaataa                                                                                          >  1:52404/1‑301 (MQ=255)
                                                                                                                                                                                                                 gAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACGTATGAAAACCGCGCCGGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCg                                                                                       <  1:517240/301‑1 (MQ=255)
                                                                                                                                                                                                                  aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTAGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTCCTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGATAATGCAGCCTTATCCAGCATATAATCACCACCCATGCGCTTATCCCTCAACACGCAACAGACCCAACAAAAGCGGCGGGCACCACCACCATTCCGTTGCAATACCACAGCAACACCACGTAAACAAAAACGGGacaaaccgc                                                                                      >  1:488363/1‑293 (MQ=255)
                                                                                                                                                                                                                                gAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGCAAACAAAAAAGGGAATAATCAACAAAAAACTaaaaa                                                                        >  1:123327/1‑301 (MQ=255)
                                                                                                                                                                                                                                 aatgtggACAACAGGTGCTGAATAAATCTCTGATGAAATCCGCGTCTGGCGGGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAat                                                                       <  1:467673/295‑1 (MQ=255)
                                                                                                                                                                                                                                           acaGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTgcg                                                             <  2:514891/301‑1 (MQ=255)
                                                                                                                                                                                                                                                       tAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAACGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATAGTTGCGCGTAAGACTTCa                                                 >  1:426402/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                           cgcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCCCTGGcaac                             >  1:414612/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     aaaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgc                   <  2:49508/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        gtaagTAATCGACCAGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgcgca                >  2:392724/3‑300 (MQ=255)
                                                                                                                                                                                                                                                                                              gATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGCTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCAATAAgc          >  2:592961/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                  gtccTGGTCTATTACTTACTTACCATTTGCTTTGTCACCGGACGCAGTTTCCCACTGGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGCTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCaa      <  1:592200/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      tGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgcgcttatcgcgcattata  <  2:477873/301‑16 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                             
ATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGT  >  NC_000913/1681328‑1681922

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