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→CGG) adePadenine 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

CTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAG  >  NC_000913/3895262‑3895707
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
ccagcGATTTTAGTGAGCTTCGAAAAAGACAATCTTCAGGATAAAAAGCAGCTCAACGATGATTACGCACGGTCTTAGATCACGCAAACGTCGGGTACCGATTGTCATCACGCAGGAGGAGATAAAGCCCAGCGCGATACCTGCGGTAATCGAGAGGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGc                                                                                                                                                                                                                                <  2:107117/224‑1 (MQ=255)
   tCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAACCGAGAAGCTGAACGGCA‑TCATCACGGCGGTCATAAACGCCGGAAACGCATCAATAGAATCCTGCCAGATCACCCGTGCCAGa                                                                                                                                                                                                                            >  1:229079/1‑229 (MQ=255)
               aGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACCCCGGAAACGATTTAGTAAGAACCTGCCCGTTCACCCATGcccgcacggaaagcacc                                                                                                                                                                                                                >  1:281541/1‑214 (MQ=255)
               aGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGAGGTAACAAACGCCGGGCCAAGTTCAGAGAACACCTGCCAGTTCACGGCTGCCAGCCTGGGAGacaca                                                                                                                                                                                                                >  1:332859/1‑227 (MQ=255)
                      taaagggcacccTTCAAGAAAAACGGCAGCGAAACGATGATGACGCCGGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGc                                                                                                                                                                                                         <  2:208896/218‑1 (MQ=255)
                      aTAAAGCAAATCTTCAGGATAAACAGCAGCGCAACGATGATTCCGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGc                                                                                                                                                                                                         <  1:556932/229‑1 (MQ=255)
                            aaaactttCAGGAAAAACGGCAGGGCACAGATGATTACGCACGGGCTTAAATCACGCCGACGCGCGGTACCGATTTTCAGCAAGCAGTAGGAGATAAAGCCCAGCGCGATACCTGCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACg                                                                                                                                                                                                    <  2:469656/222‑1 (MQ=255)
                            aCAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACg                                                                                                                                                                                                    <  1:135310/228‑1 (MQ=255)
                                            aCAGCGGCGCAAGGATAATTGCGCCCGGTCTTAGATCACGCGGACGGCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCATCGCGATACCTCCGGTAATCGATAAGTTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGc                                                                                                                                                                                   <  2:592289/229‑1 (MQ=255)
                                            aCAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGATGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCCGAACAGATTCAGTAAAATCCTGCCAGGTCACCCGTGCCAGACTGGAAGACATCAGCAAGCCCACGTAAATCAACGCACCaa                                                                                                                                                                                    >  1:421355/1‑227 (MQ=255)
                                            aCAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCCGTAAAATCCTGCCAGTTCACGCGTGCCCGACTGGGAGCAATCAAAACGcccaactgaactctggtccttt                                                                                                                                                                                    >  1:208896/1‑208 (MQ=255)
                                                   gcgcaCGATGATTACGCACGGGCTTAGATCCAGCAGACGTCCGGTCCCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCg                                                                                                                                                                           <  1:447210/225‑1 (MQ=255)
                                                   cgcAACGATGATGACGCACGGGCTTAGAGCACGCAGACGTCCGGTTCCGATTTTCATCACTCAGTAGGAGATAAAGCCCAGGGCGATACCTTCGGTAATCGATAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGAGTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCGAACGTAAAGCAGCGCACCAGCTGCA‑GCg                                                                                                                                                                           <  1:74414/229‑1 (MQ=255)
                                                           tgagTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGAGAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCt                                                                                                                                                                     <  1:328866/227‑1 (MQ=255)
                                                              ttACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTGCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCa                                                                                                                                                                 <  1:134352/228‑1 (MQ=255)
                                                                                  cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTACGTAAGATCCTGCCAGTTCACGCGTGCCCGAATGGAAGGCATCAGCACGCCAACGCAAATAAACGCACCACGTCCA‑GCGTAGCCTGGGCCCACCCCCCCCACGCGCGACa                                                                                                                                            >  2:137367/1‑229 (MQ=255)
                                                                                  cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAAATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGGCATCAGCACGCCCACATAAATCAGCGCACCACCCGCA‑AGGTAACCTGGCCCCCTACCCGCCAGCGGCCAc                                                                                                                                             >  1:588321/1‑228 (MQ=255)
                                                                                  cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACCCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAAACCTGAAGTAATCAGCACCCCAACGTAAATCAGCGGACCAGCTGCC‑GCGGAGCCTGGCACCATCCCCCCCAacggcggc                                                                                                                                             >  1:678999/1‑226 (MQ=255)
                                                                                            cGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACCCCGGAACAGATTCAGTAAGAACCTGCCAGTTCAAGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGCAGCCTGGCACCAACCCCGCCAACGGCGACAGAAAGATAAc                                                                                                                                  >  1:434197/1‑229 (MQ=255)
                                                                                                cccGTTTTTAACCGCGCAGGAGGAAAAAAAGCCACGCGCGATCCCTGCGGAAATCGAGATGCTGAAAGGCA‑TCACACCGGCGGTAATAAACCCCGGAACAGATTCAGTAAGATCCTGCCATGTCCCGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCGGCGCACCAGCTGCA‑GCGTAGCCTGGCACCAGCCCCGCCAGCGGCGACAGAAAGATAACcagc                                                                                                                              <  2:421355/228‑1 (MQ=255)
                                                                                                  cGATTTTCAACACGCAGTAGGGGATAAAGCCCAGCGCGATCCCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGACCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACcagcag                                                                                                                            <  1:687217/229‑1 (MQ=255)
                                                                                                    aTTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAAAACCGGGAAAGATTCACTAAAATCCTGCCAGGTCACCCGTGCCAGACTGGAAGTCATCAGCACGCCACCGTAAATCACCGCACCAGCTGCA‑GGGTACCCTGGCACCATCTCCCCCACCGGCGACAGAAAAATAAAcagcagaa                                                                                                                          >  1:482910/1‑227 (MQ=255)
                                                                                                          ctcccgccatggaggaTAAAACCGCGCGCGATCCATTGGGTAAACGAAAAGCAAACGGGTACTCA‑CAGGGGGGAAAAAACCGCCGAAAAAGATTAGTTAAGCTCCGGCCATTTCCCGCGTGCCAGGCTGGAAGTCATCAGCCCCCCACCGTAAACCAGCGCACCAGCTGCA‑GCGTAGCCGGGCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAAcagca                                                                                                                   <  2:355738/215‑1 (MQ=255)
                                                                                                                           aaaGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAGCGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACCCGTGCCAGACCGGAAAGAATCAGCACCCCCACGTAAACCAGCCCACCCCCTGCC‑GCGGTACCTGGCACCCTCCCCCCACACGGCGACCAGAAGAAAACCACCAG‑GAACAGAAAACCCACCAACCAcgcc                                                                                                  >  2:676900/1‑226 (MQ=255)
                                                                                                                             aGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACCCCGGAACAGATTTAGTAAGATCCTTCCCGTTCACCCGGGCCAGACTGGAAGAGCTCACCACGCCAACGTAAATCAGCCCACCCGCCGCA‑GCGGAACCTGGGACCATCCCCGCCCCCGGCGGACGAAAAAAAACCAACAG‑GAAcaaaaaacccaaaaaccacgcccg                                                                                                >  2:606471/1‑207 (MQ=255)
                                                                                                                                  aGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGGAATAAACCCCGGGCAAGAATCAGTAGAATTCAGTAAGAGCCACCTTGGCAAAATGGAAAGAATACGCACGCCCAGCGAAATCAGCGCAGCA‑CCGCACACGCAGCAGTGCACCATCCCCGCAAGCCGCGACCGAAAGAGAACCAGCAGAAAACAGCAAC‑Caaacaaaacagccgcaaaac                                                                                           >  2:331802/1‑215 (MQ=255)
                                                                                                                                    cgcgATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAAATCCTGCCAGTTCACGCGTGCTAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCCTCA‑GCGTAGCCTGGCACCATCCCCGCCACCGGCGACCAAAAGATAACCACCAG‑GAAAAGCAGA‑Cacacaaacccgtctcctatacc                                                                                         >  1:254405/1‑212 (MQ=255)
                                                                                                                                       cttACCTTGGGTAATACAAAAGGTAAAGGGCA‑TCATCCCGGCGGTAAAAAACCGCGAAAAAGATTCAGTAAGATCCTGCCAGTTCAGCGGTGCCAGCCTGGAAGTAACCAGCCCGCCAACGTAAATCAGGGAACCAGCTGAA‑GCGTAGCGTGCCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTa                                                                                    <  2:283897/227‑1 (MQ=255)
                                                                                                                                                     tCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGGACAGATTCAGTAAGATCCTGCCAGTTCACCCGGGCCAGACTGGAAAACATCAGCACCCCAACGGAAATCAACGCACCAGCTGCA‑GCGTAGCCTGGGACCCTCCCCCCCAGCGGGGCACGAAAGAAAACCAGCAG‑GAAAAGGAGA‑CCACAAACCAC‑CGCCGCCAAACCGGAACAGCCCCCAAACGa                                                                      >  2:677167/1‑229 (MQ=255)
                                                                                                                                                     tCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGAATCAGAAAGAGCCTGCCCGTTCACCCGTGCCAGACTGGAAGTCCTCAGAACGCCAACGTAAATCACCGCACCCGCGTCC‑CCGGACACTGGCCCCCTCCCCGCCAGCGGCGCCCAAAAAAAAACCCGCAG‑GCACACCAGA‑ACAAC‑ACCACACCCCGCCAAACCCGGACGAACGCCAAccaa                                                                      >  2:671825/1‑227 (MQ=255)
                                                                                                                                                     tCGAGAACCTGAACGGCA‑TCACCCCGGCGGTAAAAACCGCCGGAACAAATTCAGTAAGCTCCTGCCAGTTCCCGCGTCCCAGCCTGGACGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGa                                                                      <  1:627263/229‑1 (MQ=255)
                                                                                                                                                       gaaaaGATGAACGGTC‑TCATCCGGGGGGAAATAACCCCAGAAACAGCTGCAGTAAGACCCTGCCAATTCCGGCGGCCCAGGCGGGAAGCCCTCAGCCCGCAAAGCTAAACCGGCGCACAAGCTGCA‑GCGTAGCCGGCCCCCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGATa                                                                    <  2:281541/226‑1 (MQ=255)
                                                                                                                                                                         catcaCGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCACCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACt                                                     <  1:412512/227‑1 (MQ=255)
                                                                                                                                                                                    aaaaacaCGCAGAAAAAGTTTGAGAAAGATCCAGCGAGTTCCCCCGTGCAAACGTGAAAGTCAACGACACTAAGACGTAAAGCAGGGCACCACGAGGA‑GCGTAGCCTGGCCCAGTCCCCGCCAGCGGCGACAGAAAGATACCCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGACAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGt                                        <  2:500609/223‑1 (MQ=255)
                                                                                                                                                                                       taaCCGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCCCGCGTGCCAGACTGGAAGCCATCAGCCCGCCAACGTAAACCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAAc                                     <  1:137367/229‑1 (MQ=255)
                                                                                                                                                                                              cGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGACAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGAACGACAGCCAACCGATACGCCGGAAAAGGACTCACTATAAGCCGTAACGgaagaa                              >  1:692834/1‑229 (MQ=255)
                                                                                                                                                                                                gAACAGATTCAGTAAGATCCTGCCACTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACGACCAC‑CGCCGTCAAACCGGTACGACCGCCAACCGATACCCCGACAGAGAACGCAATATAAGCCGGAACGgcagaaaa                            >  2:562993/1‑227 (MQ=255)
                                                                                                                                                                                                          aGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCCGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGAACCAAAATAAGCCGAAACGGAAAAAAGACAGATAAAc                  >  2:489743/1‑229 (MQ=255)
                                                                                                                                                                                                           gTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGATACCCCGGAAGACGACTCCCACTAAGACGTAACGGAACAAGTACGGATAAACg                 >  1:250871/1‑229 (MQ=255)
                                                                                                                                                                                                                    tcccAGTCCACGCGGGCCAGCCTGGAAGTCACCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCGGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCgg           <  1:564669/224‑1 (MQ=255)
                                                                                                                                                                                                                    tGCCGGTTCACGCGTGCACGCCGGGAAGTCACCAGCACGCCAACTTAAACCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAAcagca                                                                                                                   <  2:569618/124‑1 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAAcagca                                                                                                                   >  1:569619/1‑124 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGCCCGCCAACCGATACGCCGGAACAGGAATCAATATAAGCCGCAACGGAAAAAGGACCCGTAAACGAACAGGACa        >  2:564668/1‑229 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACgg                                   >  1:541182/1‑202 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGCCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGCCTCCATCTACGCCGTACCGGCAGCAGTACCGATAAATGAACCGGTCa        >  1:248453/1‑229 (MQ=255)
                                                                                                                                                                                                                      ccAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCATGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TCCCGTAAGAACGGTACGACCGCCCCCCGATACACCGGCAGCGGACTCCACATACGCCGCAACGGCAGAAGTACCACTAACAGAACCgggaaca      >  2:243748/1‑224 (MQ=255)
                                                                                                                                                                                                                        gtgtcccgcgTGCCAGGCGGGACGACGACAGCACCCCAACGAAAAGCGGCGCAACAGGTGCG‑GCGAAGTCGGGCACCAACCCCGGCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGC‑CCAACAACCAC‑GGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGa    <  2:16193/223‑1 (MQ=255)
                                                                                                                                                                                                                          ttCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCACCAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGAAACGCCGGAAAGGGACCCAAATTAAGCCGTAACGGAGAAAGAACCAGTAAACAAAACGGGCACGGaaa  >  2:103503/1‑228 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
CTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAACAGCAGCGCAACGATGATTACGCACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCA‑TCATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCAGTTCACGCGTGCCAGACTGGAAGTCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGCA‑GCGTAGCCTGGCACCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAG‑GAACAGCAGA‑CCAACAACCAC‑TGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAG  >  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: ‑