Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,895,484 G→C 21.6% R375R (CGC→CGG adeP ← adenine permease, high affinity, adenine:H+ symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,895,4840GC21.6% 88.1 / 2.8 47R375R (CGC→CGGadePadenine permease, high affinity, adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base G (17/19);  new base C (9/2);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 8.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAG  >  NC_000913/3895262‑3895707
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
ccagcGATTTTAGTGAGCTTCGAAAAAGACAATCTTCAGGATAAAAAGCAGCTCAACGATGATTACGCACGGTCTTAGATCACGCAAACGTCGGGTACCGATTGTCATCACGCAGGAGGAGATAAAGCCCAGCGCGATACCTGCGGTAATCGAGAGGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGc                                                                                                                                                                                                                                <  2:107117/224‑1 (MQ=255)
   tCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAACCGAGAAGCTGAACGGCATCATCACGGCGGTCATAAACGCCGGAAACGCATCAATAGAATCCTGCCAGATCACCCGTGCCAGa                                                                                                                                                                                                                            >  1:229079/1‑229 (MQ=255)
               aGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACCCCGGAAACGATTTAGTAAGAACCTGCCCGTTCACCCATGcccgcacggaaagcacc                                                                                                                                                                                                                >  1:281541/1‑214 (MQ=255)
               aGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGAGGTAACAAACGCCGGGCCAAGTTCAGAGAACACCTGCCAGTTCACGGCTGCCAGCCTGGGAGacaca                                                                                                                                                                                                                >  1:332859/1‑227 (MQ=255)
                      taaagggcacccTTCAAGAAAAACGGCAGCGAAACGATGATGACGCCGGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGc                                                                                                                                                                                                         <  2:208896/218‑1 (MQ=255)
                      aTAAAGCAAATCTTCAGGATAAACAGCAGCGCAACGATGATTCCGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGc                                                                                                                                                                                                         <  1:556932/229‑1 (MQ=255)
                            aaaactttCAGGAAAAACGGCAGGGCACAGATGATTACGCACGGGCTTAAATCACGCCGACGCGCGGTACCGATTTTCAGCAAGCAGTAGGAGATAAAGCCCAGCGCGATACCTGCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACg                                                                                                                                                                                                    <  2:469656/222‑1 (MQ=255)
                            aCAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACg                                                                                                                                                                                                    <  1:135310/228‑1 (MQ=255)
                                            aCAGCGGCGCAAGGATAATTGCGCCCGGTCTTAGATCACGCGGACGGCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCATCGCGATACCTCCGGTAATCGATAAGTTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGc                                                                                                                                                                                   <  2:592289/229‑1 (MQ=255)
                                            aCAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGATGCTGAACGGCATCATCACGGCGGTAATAAACGCCCGAACAGATTCAGTAAAATCCTGCCAGGTCACCCGTGCCAGACTGGAAGACATCAGCAAGCCCACGTAAATCAACGCACCaa                                                                                                                                                                                    >  1:421355/1‑227 (MQ=255)
                                            aCAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCCGTAAAATCCTGCCAGTTCACGCGTGCCCGACTGGGAGCAATCAAAACGcccaactgaactctggtccttt                                                                                                                                                                                    >  1:208896/1‑208 (MQ=255)
                                                   gcgcaCGATGATTACGCACGGGCTTAGATCCAGCAGACGTCCGGTCCCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCg                                                                                                                                                                           <  1:447210/225‑1 (MQ=255)
                                                   cgcAACGATGATGACGCACGGGCTTAGAGCACGCAGACGTCCGGTTCCGATTTTCATCACTCAGTAGGAGATAAAGCCCAGGGCGATACCTTCGGTAATCGATAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGAGTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCGAACGTAAAGCAGCGCACCAGCTGCAGCg                                                                                                                                                                           <  1:74414/229‑1 (MQ=255)
                                                           tgagTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGAGAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCt                                                                                                                                                                     <  1:328866/227‑1 (MQ=255)
                                                              ttACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTGCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCa                                                                                                                                                                 <  1:134352/228‑1 (MQ=255)
                                                                                  cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTACGTAAGATCCTGCCAGTTCACGCGTGCCCGAATGGAAGGCATCAGCACGCCAACGCAAATAAACGCACCACGTCCAGCGTAGCCTGGGCCCACCCCCCCCACGCGCGACa                                                                                                                                            >  2:137367/1‑229 (MQ=255)
                                                                                  cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAAATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGGCATCAGCACGCCCACATAAATCAGCGCACCACCCGCAAGGTAACCTGGCCCCCTACCCGCCAGCGGCCAc                                                                                                                                             >  1:588321/1‑228 (MQ=255)
                                                                                  cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACCCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAAACCTGAAGTAATCAGCACCCCAACGTAAATCAGCGGACCAGCTGCCGCGGAGCCTGGCACCATCCCCCCCAacggcggc                                                                                                                                             >  1:678999/1‑226 (MQ=255)
                                                                                            cGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACCCCGGAACAGATTCAGTAAGAACCTGCCAGTTCAAGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGCAGCCTGGCACCAACCCCGCCAACGGCGACAGAAAGATAAc                                                                                                                                  >  1:434197/1‑229 (MQ=255)
                                                                                                cccGTTTTTAACCGCGCAGGAGGAAAAAAAGCCACGCGCGATCCCTGCGGAAATCGAGATGCTGAAAGGCATCACACCGGCGGTAATAAACCCCGGAACAGATTCAGTAAGATCCTGCCATGTCCCGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCGGCGCACCAGCTGCAGCGTAGCCTGGCACCAGCCCCGCCAGCGGCGACAGAAAGATAACcagc                                                                                                                              <  2:421355/228‑1 (MQ=255)
                                                                                                  cGATTTTCAACACGCAGTAGGGGATAAAGCCCAGCGCGATCCCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGACCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACcagcag                                                                                                                            <  1:687217/229‑1 (MQ=255)
                                                                                                    aTTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAAAACCGGGAAAGATTCACTAAAATCCTGCCAGGTCACCCGTGCCAGACTGGAAGTCATCAGCACGCCACCGTAAATCACCGCACCAGCTGCAGGGTACCCTGGCACCATCTCCCCCACCGGCGACAGAAAAATAAAcagcagaa                                                                                                                          >  1:482910/1‑227 (MQ=255)
                                                                                                          ctcccgccatggaggaTAAAACCGCGCGCGATCCATTGGGTAAACGAAAAGCAAACGGGTACTCACAGGGGGGAAAAAACCGCCGAAAAAGATTAGTTAAGCTCCGGCCATTTCCCGCGTGCCAGGCTGGAAGTCATCAGCCCCCCACCGTAAACCAGCGCACCAGCTGCAGCGTAGCCGGGCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAAcagca                                                                                                                   <  2:355738/215‑1 (MQ=255)
                                                                                                                           aaaGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAGCGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACCCGTGCCAGACCGGAAAGAATCAGCACCCCCACGTAAACCAGCCCACCCCCTGCCGCGGTACCTGGCACCCTCCCCCCACACGGCGACCAGAAGAAAACCACCAGGAACAGAAAACCCACCAACCAcgcc                                                                                                  >  2:676900/1‑226 (MQ=255)
                                                                                                                             aGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACCCCGGAACAGATTTAGTAAGATCCTTCCCGTTCACCCGGGCCAGACTGGAAGAGCTCACCACGCCAACGTAAATCAGCCCACCCGCCGCAGCGGAACCTGGGACCATCCCCGCCCCCGGCGGACGAAAAAAAACCAACAGGAAcaaaaaacccaaaaaccacgcccg                                                                                                >  2:606471/1‑207 (MQ=255)
                                                                                                                                  aGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGGAATAAACCCCGGGCAAGAATCAGTAGAATTCAGTAAGAGCCACCTTGGCAAAATGGAAAGAATACGCACGCCCAGCGAAATCAGCGCAGCACCGCACACGCAGCAGTGCACCATCCCCGCAAGCCGCGACCGAAAGAGAACCAGCAGAAAACAGCAACCaaacaaaacagccgcaaaac                                                                                           >  2:331802/1‑215 (MQ=255)
                                                                                                                                    cgcgATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAAATCCTGCCAGTTCACGCGTGCTAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCCTCAGCGTAGCCTGGCACCATCCCCGCCACCGGCGACCAAAAGATAACCACCAGGAAAAGCAGACacacaaacccgtctcctatacc                                                                                         >  1:254405/1‑212 (MQ=255)
                                                                                                                                       cttACCTTGGGTAATACAAAAGGTAAAGGGCATCATCCCGGCGGTAAAAAACCGCGAAAAAGATTCAGTAAGATCCTGCCAGTTCAGCGGTGCCAGCCTGGAAGTAACCAGCCCGCCAACGTAAATCAGGGAACCAGCTGAAGCGTAGCGTGCCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTa                                                                                    <  2:283897/227‑1 (MQ=255)
                                                                                                                                                     tCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGGACAGATTCAGTAAGATCCTGCCAGTTCACCCGGGCCAGACTGGAAAACATCAGCACCCCAACGGAAATCAACGCACCAGCTGCAGCGTAGCCTGGGACCCTCCCCCCCAGCGGGGCACGAAAGAAAACCAGCAGGAAAAGGAGACCACAAACCACCGCCGCCAAACCGGAACAGCCCCCAAACGa                                                                      >  2:677167/1‑229 (MQ=255)
                                                                                                                                                     tCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGAATCAGAAAGAGCCTGCCCGTTCACCCGTGCCAGACTGGAAGTCCTCAGAACGCCAACGTAAATCACCGCACCCGCGTCCCCGGACACTGGCCCCCTCCCCGCCAGCGGCGCCCAAAAAAAAACCCGCAGGCACACCAGAACAACACCACACCCCGCCAAACCCGGACGAACGCCAAccaa                                                                      >  2:671825/1‑227 (MQ=255)
                                                                                                                                                     tCGAGAACCTGAACGGCATCACCCCGGCGGTAAAAACCGCCGGAACAAATTCAGTAAGCTCCTGCCAGTTCCCGCGTCCCAGCCTGGACGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGa                                                                      <  1:627263/229‑1 (MQ=255)
                                                                                                                                                       gaaaaGATGAACGGTCTCATCCGGGGGGAAATAACCCCAGAAACAGCTGCAGTAAGACCCTGCCAATTCCGGCGGCCCAGGCGGGAAGCCCTCAGCCCGCAAAGCTAAACCGGCGCACAAGCTGCAGCGTAGCCGGCCCCCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATa                                                                    <  2:281541/226‑1 (MQ=255)
                                                                                                                                                                         catcaCGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCACCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACt                                                     <  1:412512/227‑1 (MQ=255)
                                                                                                                                                                                    aaaaacaCGCAGAAAAAGTTTGAGAAAGATCCAGCGAGTTCCCCCGTGCAAACGTGAAAGTCAACGACACTAAGACGTAAAGCAGGGCACCACGAGGAGCGTAGCCTGGCCCAGTCCCCGCCAGCGGCGACAGAAAGATACCCAGCAGGAACAGCAGACCAACAACCACTGCCGACAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGt                                        <  2:500609/223‑1 (MQ=255)
                                                                                                                                                                                       taaCCGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCCCGCGTGCCAGACTGGAAGCCATCAGCCCGCCAACGTAAACCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAAc                                     <  1:137367/229‑1 (MQ=255)
                                                                                                                                                                                              cGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGACAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGAACGACAGCCAACCGATACGCCGGAAAAGGACTCACTATAAGCCGTAACGgaagaa                              >  1:692834/1‑229 (MQ=255)
                                                                                                                                                                                                gAACAGATTCAGTAAGATCCTGCCACTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACGACCACCGCCGTCAAACCGGTACGACCGCCAACCGATACCCCGACAGAGAACGCAATATAAGCCGGAACGgcagaaaa                            >  2:562993/1‑227 (MQ=255)
                                                                                                                                                                                                          aGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCCGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGAACCAAAATAAGCCGAAACGGAAAAAAGACAGATAAAc                  >  2:489743/1‑229 (MQ=255)
                                                                                                                                                                                                           gTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACCCCGGAAGACGACTCCCACTAAGACGTAACGGAACAAGTACGGATAAACg                 >  1:250871/1‑229 (MQ=255)
                                                                                                                                                                                                                    tcccAGTCCACGCGGGCCAGCCTGGAAGTCACCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCGGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCgg           <  1:564669/224‑1 (MQ=255)
                                                                                                                                                                                                                    tGCCGGTTCACGCGTGCACGCCGGGAAGTCACCAGCACGCCAACTTAAACCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAAcagca                                                                                                                   <  2:569618/124‑1 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAAcagca                                                                                                                   >  1:569619/1‑124 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGCCCGCCAACCGATACGCCGGAACAGGAATCAATATAAGCCGCAACGGAAAAAGGACCCGTAAACGAACAGGACa        >  2:564668/1‑229 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACgg                                   >  1:541182/1‑202 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGCCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGCCTCCATCTACGCCGTACCGGCAGCAGTACCGATAAATGAACCGGTCa        >  1:248453/1‑229 (MQ=255)
                                                                                                                                                                                                                      ccAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCATGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTCCCGTAAGAACGGTACGACCGCCCCCCGATACACCGGCAGCGGACTCCACATACGCCGCAACGGCAGAAGTACCACTAACAGAACCgggaaca      >  2:243748/1‑224 (MQ=255)
                                                                                                                                                                                                                        gtgtcccgcgTGCCAGGCGGGACGACGACAGCACCCCAACGAAAAGCGGCGCAACAGGTGCGGCGAAGTCGGGCACCAACCCCGGCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGCCCAACAACCACGGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGa    <  2:16193/223‑1 (MQ=255)
                                                                                                                                                                                                                          ttCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCACCAACCACTGCCGTCAGACCGGTACGACCGCCAACCGAAACGCCGGAAAGGGACCCAAATTAAGCCGTAACGGAGAAAGAACCAGTAAACAAAACGGGCACGGaaa  >  2:103503/1‑228 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
CTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCATCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCAGCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAACCACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAG  >  NC_000913/3895262‑3895707

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