Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,681,614 G→A 14.3% V213V (GTG→GTA 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,6140GA14.3% 99.5 / 8.6 42V213V (GTG→GTAfolMdihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base A (2/4);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

ATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGC  >  NC_000913/1681328‑1681914
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                             
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)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                             
ATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAGCATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGC  >  NC_000913/1681328‑1681914

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: