Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,095,838 C→A 5.1% E49* (GAA→TAA)  rhaA ← L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,095,8380CA 5.1% 754.0 / 19.0 212E49* (GAA→TAA) rhaAL‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base C (94/103);  new base A (3/7);  total (100/112)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 216 total aligned reads displayed.

TAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGC  >  NC_000913/4095690‑4095986
                                                                                                                                                    |                                                                                                                                                     
tagatGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTc                                                                                                                                                       <  2:211577/149‑1 (MQ=255)
    tGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCTGTCAGAATACCGTCTCAGTTTTCaaaa                                                                                                                                                   <  2:1566510/149‑1 (MQ=255)
      gCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAAcc                                                                                                                                                >  2:1303369/1‑149 (MQ=255)
      gCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCAGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAccc                                                                                                                                                <  1:181526/149‑1 (MQ=255)
        aTGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAAc                                                                                                                                                 <  3:1766869/146‑1 (MQ=255)
        aTGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCgg                                                                                                                                              >  4:1767322/1‑149 (MQ=255)
         tGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGa                                                                                                                                             >  3:627459/1‑149 (MQ=255)
          gTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGaa                                                                                                                                            <  3:1219338/149‑1 (MQ=255)
           tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGaag                                                                                                                                           <  4:627459/149‑2 (MQ=255)
           tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGTTTTGAAAACCCGGaag                                                                                                                                           <  1:1518629/149‑2 (MQ=255)
           tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGaaa                                                                                                                                           <  4:1154475/149‑1 (MQ=255)
           tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGaaa                                                                                                                                           >  2:290578/1‑149 (MQ=255)
           tAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGaaa                                                                                                                                           >  2:840294/1‑149 (MQ=255)
              aTTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcat                                                                                                                                        >  4:1047255/1‑149 (MQ=255)
              aTTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcat                                                                                                                                        >  1:540381/1‑149 (MQ=255)
                tAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcatca                                                                                                                                      <  1:1269851/149‑1 (MQ=255)
                         tCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGc                                                                                                                             <  1:267431/149‑1 (MQ=255)
                          cGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTg                                                                                                                              <  4:1489129/147‑1 (MQ=255)
                          cGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTg                                                                                                                              >  3:1489129/1‑147 (MQ=255)
                             ccccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagc                                                                                                                         >  2:927541/1‑149 (MQ=255)
                              cccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagca                                                                                                                        <  1:840294/149‑1 (MQ=255)
                              cccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagca                                                                                                                        >  1:882753/1‑149 (MQ=255)
                               ccGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagcag                                                                                                                       >  4:519315/1‑149 (MQ=255)
                                cGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGt                                                                                                                      <  1:290578/149‑1 (MQ=255)
                                   aaTCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCa                                                                                                                   <  3:1667554/149‑1 (MQ=255)
                                     tCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCAtt                                                                                                                 >  1:361020/1‑149 (MQ=255)
                                     tCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCAtt                                                                                                                 <  2:357323/149‑1 (MQ=255)
                                       aGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTc                                                                                                                                                       >  2:133753/1‑110 (MQ=255)
                                       aGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagcag                                                                                                                       <  3:1527658/141‑1 (MQ=255)
                                       aGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCcagcag                                                                                                                       >  4:1527658/1‑141 (MQ=255)
                                             aTAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACggg                                                                                                         <  1:429641/149‑1 (MQ=255)
                                             aTAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACggg                                                                                                         <  1:1192539/149‑1 (MQ=255)
                                                gCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGcc                                                                                                                            >  2:1114299/1‑127 (MQ=255)
                                                gCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGcc                                                                                                                            <  1:1114299/127‑1 (MQ=255)
                                                  cTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAAc                                                                                                    <  4:275570/149‑1 (MQ=255)
                                                     ttCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACatcatc                                                                                                                                    >  3:1468171/1‑114 (MQ=255)
                                                     ttCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACatcatc                                                                                                                                    <  4:1468171/114‑1 (MQ=255)
                                                       ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACgg                                                                                                          <  1:672724/138‑1 (MQ=255)
                                                       ccAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACgg                                                                                                          >  2:672724/1‑138 (MQ=255)
                                                         aGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGaaa                                                                                                             >  3:1050194/1‑133 (MQ=255)
                                                           aTCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCaa                                                                                              >  2:1495004/1‑146 (MQ=255)
                                                           aTCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCaa                                                                                              <  1:1495004/146‑1 (MQ=255)
                                                               gCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcatca                                                                                                                                      >  3:863048/1‑102 (MQ=255)
                                                                cACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAAcc                                                                                                                                                <  4:110287/91‑1 (MQ=255)
                                                                   gTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCt                                                                                                                               <  2:266879/105‑1 (MQ=255)
                                                                        tCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTCCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGcctc                                                                              >  2:482013/1‑149 (MQ=255)
                                                                         cACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAAcatca                                                                                                                                      >  1:240381/1‑92 (MQ=255)
                                                                            tGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCg                                                                          >  3:794127/1‑149 (MQ=255)
                                                                                aTTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAgcgc                                                                                 >  4:185374/1‑138 (MQ=255)
                                                                                aTTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAgcgc                                                                                 <  3:185374/138‑1 (MQ=255)
                                                                                       gcTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCAtt                                                                                                                 >  4:737261/1‑99 (MQ=255)
                                                                                       gcTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCAtt                                                                                                                 <  3:737261/99‑1 (MQ=255)
                                                                                       gcTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGa                                                                                                  <  1:1646394/114‑1 (MQ=255)
                                                                                          ttGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCAc                                                            <  4:794127/149‑1 (MQ=255)
                                                                                          ttGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTCCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCAc                                                            <  1:482013/149‑1 (MQ=255)
                                                                                               cGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCaa                                                                                              >  1:3454/1‑110 (MQ=255)
                                                                                               cGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCaa                                                                                              <  2:3454/110‑1 (MQ=255)
                                                                                               cGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCAccgccg                                                       >  1:1876834/1‑149 (MQ=255)
                                                                                                     aTTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAAc                                                 >  1:108333/1‑149 (MQ=255)
                                                                                                      ttGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAAc                                                                                                    >  1:723243/1‑97 (MQ=255)
                                                                                                      ttGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAAc                                                                                                    <  2:723243/97‑1 (MQ=255)
                                                                                                      ttGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAg                                                                                             >  3:11483/1‑104 (MQ=255)
                                                                                                      ttGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACa                                                                       >  4:883911/1‑126 (MQ=255)
                                                                                                        gCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCt                                              <  3:563073/149‑1 (MQ=255)
                                                                                                          ctatgGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCa                                                             <  4:240948/129‑1 (MQ=255)
                                                                                                          ctatgGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCa                                                             >  3:240948/4‑132 (MQ=255)
                                                                                                               gCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGa                                                                         <  2:1628225/115‑1 (MQ=255)
                                                                                                                ccTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGtt                                                                                           >  3:1253671/1‑96 (MQ=255)
                                                                                                                ccTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGtt                                                                                           <  4:1253671/96‑1 (MQ=255)
                                                                                                                 cTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATcc                                                                >  1:1625878/1‑122 (MQ=255)
                                                                                                                 cTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGc                                       <  4:837291/147‑1 (MQ=255)
                                                                                                                 cTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTa                                     >  3:1574295/1‑149 (MQ=255)
                                                                                                                 cTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTa                                     >  4:1317889/1‑149 (MQ=255)
                                                                                                                  tGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAg                                    >  2:874732/1‑149 (MQ=255)
                                                                                                                  tGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAg                                    <  3:1317889/149‑1 (MQ=255)
                                                                                                                        cccccGGTCAGCGAACCTACCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCa                              >  2:1573410/1‑149 (MQ=255)
                                                                                                                         ccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAAc                                                 >  2:856822/1‑129 (MQ=255)
                                                                                                                         ccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAAc                                                 <  1:856822/129‑1 (MQ=255)
                                                                                                                         ccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGc                                       >  1:1005573/1‑139 (MQ=255)
                                                                                                                         ccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGc                                       <  3:449607/139‑1 (MQ=255)
                                                                                                                         ccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGc                                       <  2:1005573/139‑1 (MQ=255)
                                                                                                                         ccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAg                             >  2:1852392/1‑149 (MQ=255)
                                                                                                                         ccccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAg                             >  4:1684582/1‑149 (MQ=255)
                                                                                                                          cccGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAgg                            >  1:1582488/1‑149 (MQ=255)
                                                                                                                             ggTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCt                         <  2:412997/149‑1 (MQ=255)
                                                                                                                               tCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGtt                                           <  3:255419/129‑1 (MQ=255)
                                                                                                                               tCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGtt                                           >  4:255419/1‑129 (MQ=255)
                                                                                                                                    gAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCg                                                                          >  1:1054872/1‑93 (MQ=255)
                                                                                                                                    gAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCg                                                                          <  2:1054872/93‑1 (MQ=255)
                                                                                                                                    gAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCt                                              <  1:1569468/121‑1 (MQ=255)
                                                                                                                                        tttCCGGTTTTGAAAACCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGa              >  1:650147/2‑149 (MQ=255)
                                                                                                                                        cTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGc                                       >  2:6928/1‑124 (MQ=255)
                                                                                                                                        cTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGc                                       <  1:6928/124‑1 (MQ=255)
                                                                                                                                        cTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGa              <  3:964958/149‑1 (MQ=255)
                                                                                                                                         ttCCGGTTTTGAAAACCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAg             <  1:1745957/149‑1 (MQ=255)
                                                                                                                                         ttCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGt                       <  2:563684/139‑1 (MQ=255)
                                                                                                                                          ttccggTTTTGAAAACCCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAg             >  2:1745957/5‑149 (MQ=255)
                                                                                                                                             gggTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACggg                                                                                                         <  1:1324582/53‑1 (MQ=255)
                                                                                                                                             gggTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGt         >  1:412747/1‑149 (MQ=255)
                                                                                                                                                 tttCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATg     <  4:1574295/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAGATGGCATGTAAATTAAGCCGTTTCGGCCCCGGAATCAGCCGCATAGCCTGTTCCAGATCGGCACGTAGCTCACTGGCATTACGCGCTTTGCCCGGATAATTGCCTGTGGCCTGAATCCCCCCGGTCAGCGAACCTTCCGGGTTTTCAAAACCGGAAACATCATCGCCCTGCCAGCAGTGCATTGAAACGGGTAAACGATCAAGTTGGCGCAGCGCCTCCTCGACATCAATCCCCACCGCCGCGAAACGCTGTTTCGCTAGCTCCCAGGCCTGTTCCAGTTGAGTGGTCATGCGC  >  NC_000913/4095690‑4095986

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