Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,004,029 C→A 8.0% M50I (ATG→ATT)  rarD ← putative chloramphenical resistance permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,0290CA 8.0% 957.3 / 6.4 296M50I (ATG→ATT) rarDputative chloramphenical resistance permease
Reads supporting (aligned to +/- strand):  ref base C (177/84);  new base A (2/33);  total (179/117)
Fisher's exact test for biased strand distribution p-value = 6.05e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

AGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCC  >  NC_000913/4003880‑4004176
                                                                                                                                                      |                                                                                                                                                   
   atatGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATa                                                                                                                                                  <  2:789313/150‑1 (MQ=255)
     atGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATaaa                                                                                                                                                <  2:4393105/150‑1 (MQ=255)
       ggTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGACTCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGa                                                                                                                                              >  2:5790090/1‑150 (MQ=255)
       ggTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGa                                                                                                                                              <  2:4041699/150‑1 (MQ=255)
        gTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGaa                                                                                                                                             >  1:5945370/1‑150 (MQ=255)
         tGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGaaa                                                                                                                                            >  2:2542488/1‑150 (MQ=255)
          gATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGaaaa                                                                                                                                           >  1:6046654/1‑150 (MQ=255)
            ttGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAAc                                                                                                                                         <  2:6171177/150‑1 (MQ=255)
             tGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACg                                                                                                                                        <  2:832785/150‑1 (MQ=255)
              gTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGa                                                                                                                                       >  1:5484793/1‑150 (MQ=255)
                tCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGAcc                                                                                                                                     >  2:373800/1‑150 (MQ=255)
                    cGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGAt                                                                                                                                 >  2:3902274/1‑150 (MQ=255)
                     gCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGGGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATc                                                                                                                                >  2:5108367/1‑150 (MQ=255)
                     gCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATc                                                                                                                                >  1:4016549/1‑150 (MQ=255)
                     gCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATc                                                                                                                                >  2:3417294/1‑150 (MQ=255)
                         ataTAAACAGTAGCCAGTTGCAACCAATCAGCACGGCAGAGACTGCCAGAATAAAAATTTTCTGTGGCGTCTTAATCAGCTTTTTTAAATAGGACCACTGGCGGAAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCAcgc                                                                                                                            <  2:1662321/148‑1 (MQ=255)
                              aaCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGc                                                                                                                       >  1:2964022/1‑150 (MQ=255)
                              aaCACTAGCAATTTCCCACCAAACACCACGGCACACACAGCCACAAAAAAAATTTTACGTGCCGACACAAACAGCTTTTTAAAAAAGCACCACTGGCGCAAAATTAACAT‑AACAACAACAATAAACAAAAACGACCAGATCACGCGATGc                                                                                                                       <  2:2784592/150‑1 (MQ=255)
                               aCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCg                                                                                                                      >  2:3285241/1‑150 (MQ=255)
                                cAGTAGCCTGTTGCCACAACTCAGCACGGCAGAGACTGCCAGAATAAAAATTTTCTGTGGCGTCTTAATCAGCTTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGt                                                                                                                     <  2:1519582/150‑1 (MQ=255)
                                cAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CCGCACCACCATAAAGACAGACGACCAGATCACGCGATGCgg                                                                                                                     >  1:2763234/1‑149 (MQ=255)
                                cAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGt                                                                                                                     >  2:2833846/1‑150 (MQ=255)
                                 aGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTc                                                                                                                    >  2:1005774/1‑150 (MQ=255)
                                 aGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTc                                                                                                                    >  2:171451/1‑150 (MQ=255)
                                     gCCTGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAATTTTTCTGGGGCTTCTGAATCAGCTTTTTTAAATAGCACCCCTGGCGGAAAATGCTCAT‑CAGCACCACCATAAACAAAAACGACCAGATCACGCGATGCGTCAGGa                                                                                                                <  2:2113804/150‑1 (MQ=255)
                                     gCCAGTGGCCCCCACTCCGCACGGCAGAGACTGCCCGCATAAAATTTTTCTGGGGCGTCGGACTCCGCGTTTTTAATTCGGACCACTGGCGGCACTTGCTCAT‑CCGCACCCCCATAACGAACAACGACCAGATCACGCGATGCGTCAGGa                                                                                                                <  2:4332336/150‑1 (MQ=255)
                                      ccAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGAt                                                                                                               >  2:6101735/1‑150 (MQ=255)
                                      ccAGTTGCCACCAATCAGCACGGCAGAGACTGCAAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGAACCACAATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGAt                                                                                                               >  1:4919777/1‑150 (MQ=255)
                                              caccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATcgg                                                                                                       >  1:4326710/1‑150 (MQ=255)
                                                 caATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATcggcgg                                                                                                    >  1:718119/1‑150 (MQ=255)
                                                    tCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCa                                                                                                 >  2:3542606/1‑150 (MQ=255)
                                                     cAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAATAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAc                                                                                                <  2:4508253/150‑1 (MQ=255)
                                                          cGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACgtagt                                                                                           >  2:5939512/1‑150 (MQ=255)
                                                           gacAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACgtagta                                                                                          >  2:932564/3‑150 (MQ=255)
                                                             cAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTTGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTaaa                                                                                        >  1:4479144/1‑150 (MQ=255)
                                                             cAGAGACGGCCAGCAAAACATTTTTCGGTGGGGTCTGACTCTGGTTTTTTAATTAGGACCACTGGCGGCAAATGCCCAG‑AACCACCACCTTAAATAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTaaa                                                                                        <  2:6554153/150‑1 (MQ=255)
                                                               tagaCTCCCAGCATAAAATTTTTCTGGGGCGTCTGAATCTGCTTTTTTAAATCGGACCACGGGCGGAAAATCCACAA‑CACCACCACCATAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATc                                                                                      <  2:3097837/149‑1 (MQ=255)
                                                                  aCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAc                                                                                   >  1:1334583/1‑150 (MQ=255)
                                                                      caacaataataTTTTTGTGGGGGGTCTGATTCAGTTTTTTTAATTCGGACCCCGGGCTGAAATGGCTAAT‑CTCCACCAACATAATTAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGa                                                                               <  1:4226014/145‑1 (MQ=255)
                                                                       cAGCATAAAATTTTTGTGTGGCGTCTGAATCAGCTTTTTTAAATAGGACCCCTGGCGGAAAAGGCTCAT‑CCGCACCACCATAAATAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGaa                                                                              <  2:226276/150‑1 (MQ=255)
                                                                       cAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGaa                                                                              >  2:3056182/1‑150 (MQ=255)
                                                                       cAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGaa                                                                              >  1:1564451/1‑150 (MQ=255)
                                                                        tgaaTAAAAATTTTGTGTGGCGTGTGACTCAGCTTTTTTAAAGAGGACCCGTGGCGGAAAATGC‑CATCCCGCACCACCATAAAGACCCACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAg                                                                             <  2:5841363/147‑1 (MQ=255)
                                                                        agaaTAAAATTTTTGTGTGGCGTCTTAATCAGGTTTTTTAAATAGGACCACTGGCGGAAATTGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAg                                                                             <  2:6046654/147‑1 (MQ=255)
                                                                            tAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                         >  1:438579/1‑150 (MQ=255)
                                                                             aaaatTTTTCTGTGGCTTCTTAATCTGCTTTTTAAAATTGCACCCCTGGCGCAAAATCCACAA‑CCCCCCCACCAAAAACAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGc                                                                        <  1:3220320/150‑1 (MQ=255)
                                                                              aaatTTTTTTGTGGCGTATTATTCTGTTTTTTTAAATCGCCCCAGTGGCGCAAATTGAAAAT‑CACCACCACAAAAAAAAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCt                                                                       <  1:5273179/150‑1 (MQ=255)
                                                                               aatttttCTGTGGCGTCTTAATCAGCTTTTTTAAATAGGACCACGGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTg                                                                      <  2:6378282/147‑1 (MQ=255)
                                                                                 aTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATCGGACCACGGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgc                                                                    <  2:2794837/150‑1 (MQ=255)
                                                                                  ttttCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                   >  1:1131004/1‑150 (MQ=255)
                                                                                   tttCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcgc                                                                  >  1:3466009/1‑150 (MQ=255)
                                                                                    ttgtgtgGTTTCTTAATCTGCTTTTTTAATTTGGACCCCTGGCGGAAATTGCCAAT‑CACCACCACCATAAACAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt                                                                 <  1:2838052/147‑1 (MQ=255)
                                                                                         tgGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcc                                                            >  2:3367625/1‑150 (MQ=255)
                                                                                          ggCTTCTTAATCTGCTTTTTTAATTAGCACCCCGGGCGGAAATTGCTAAT‑CAGCCCCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAccc                                                           <  2:190346/150‑1 (MQ=255)
                                                                                          ggCTTCTGAATATGCTTTTTTAATTTGGACCAGGGGCGGAAATTGCTCAT‑AAGCACCACAATAAATAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAccc                                                           <  1:6581448/150‑1 (MQ=255)
                                                                                          ggCGTTTTATTCTGCTTTTTTAATTCGGACCAGTGGCGGCAAATGATCAA‑CACCCCAACAATAAATAACAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAccc                                                           <  1:1708370/150‑1 (MQ=255)
                                                                                              tattaaTCAGCTTTTTTAAATAGGACCACTGGCGGAAACTGCTCAT‑CAGCACCACAATAAATAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaa                                                       <  2:5482252/146‑1 (MQ=255)
                                                                                              tCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaa                                                       >  2:1312430/1‑150 (MQ=255)
                                                                                               cTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaat                                                      >  1:5558219/1‑150 (MQ=255)
                                                                                                  aaTCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaa                                                   <  2:4988435/150‑1 (MQ=255)
                                                                                                       gCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                              >  1:1039324/1‑150 (MQ=255)
                                                                                                       gCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                              >  1:3047672/1‑150 (MQ=255)
                                                                                                         gTTTTTAAATCGGACCAGTGGCGGAAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                            <  2:5619404/150‑1 (MQ=255)
                                                                                                         gTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                            >  2:3223692/1‑150 (MQ=255)
                                                                                                         gTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                            >  1:2277771/1‑150 (MQ=255)
                                                                                                          tttttAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                           >  1:2166284/1‑150 (MQ=255)
                                                                                                          tttttAAATAGGACCACTGGCGGCAAATGCTCAT‑AAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                           >  1:5238904/1‑150 (MQ=255)
                                                                                                             ttAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                        >  2:5499126/1‑150 (MQ=255)
                                                                                                             ttAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                        >  1:508939/1‑150 (MQ=255)
                                                                                                                 aTAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                    >  2:4919717/1‑150 (MQ=255)
                                                                                                                  tAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCACTATACCCCAAATAAAATAAGCGGCAAGAGCg                                   >  1:931635/1‑150 (MQ=255)
                                                                                                                     gACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                >  1:5202947/1‑150 (MQ=255)
                                                                                                                     gACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                >  2:6641898/1‑150 (MQ=255)
                                                                                                                     gACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                >  1:4665576/1‑150 (MQ=255)
                                                                                                                     gACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                >  2:698940/1‑150 (MQ=255)
                                                                                                                     gACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                >  1:1874034/1‑150 (MQ=255)
                                                                                                                          cTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                           >  1:4545505/1‑150 (MQ=255)
                                                                                                                              cggcAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                       >  2:6702054/1‑150 (MQ=255)
                                                                                                                              cggaAATTGCTCAT‑CAGAACCACCATAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                       <  2:2696794/150‑1 (MQ=255)
                                                                                                                               ggcAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                      >  1:2463921/1‑150 (MQ=255)
                                                                                                                                  aaaTGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGcc                   >  2:3972960/1‑150 (MQ=255)
                                                                                                                                   aaTGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                  >  2:3578/1‑150 (MQ=255)
                                                                                                                                     tGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgcg                >  2:5328185/1‑150 (MQ=255)
                                                                                                                                      gCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt               >  2:1038122/1‑150 (MQ=255)
                                                                                                                                      gCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt               >  1:6415501/1‑150 (MQ=255)
                                                                                                                                      gCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt               <  2:4943278/150‑1 (MQ=255)
                                                                                                                                        tcat‑caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGTCGCgttt             >  2:6296996/1‑150 (MQ=255)
                                                                                                                                         caaacaccaccacCAAAAACAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt             <  1:6299525/143‑1 (MQ=255)
                                                                                                                                          at‑caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACg                                                                                               >  1:809812/1‑66 (MQ=255)
                                                                                                                                          at‑caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgt           >  1:4739048/1‑150 (MQ=255)
                                                                                                                                            tccccaccacAATAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt          <  2:2015137/146‑1 (MQ=255)
                                                                                                                                              acccccaccaCAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt        <  1:5797846/146‑1 (MQ=255)
                                                                                                                                               gCACCACCATAAAGAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTg       <  2:4794797/150‑1 (MQ=255)
                                                                                                                                                caccacCTAAACCAAACACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc      <  2:4016633/150‑1 (MQ=255)
                                                                                                                                                caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGGTTGTTTTGc      >  1:3430466/1‑150 (MQ=255)
                                                                                                                                                 accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca     >  2:5308656/1‑150 (MQ=255)
                                                                                                                                                   caccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc   >  1:1942722/1‑150 (MQ=255)
                                                                                                                                                    accaTAAATAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc  <  1:5659037/150‑1 (MQ=255)
                                                                                                                                                    accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc  >  2:2445608/1‑150 (MQ=255)
                                                                                                                                                    accaAAAACAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc  <  1:3631647/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCAT‑CAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCC  >  NC_000913/4003880‑4004176

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