Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,981,031 G→T 73.6% C384* (TGC→TGA)  araE ← arabinose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,981,0310GT73.6% 284.5 / 131.9 179C384* (TGC→TGA) araEarabinose transporter
Reads supporting (aligned to +/- strand):  ref base G (23/24);  new base T (60/70);  total (85/94)
Fisher's exact test for biased strand distribution p-value = 8.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 182 total aligned reads displayed.

GTCAGGAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGTA  >  NC_000913/2980933‑2981130
                                                                                                  |                                                                                                   
gtCAGGAAGGTCGCTCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGc                                                                                                    >  2:3751606/1‑100 (MQ=255)
  cAGGAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATcac                                                                                                  <  1:1962714/100‑1 (MQ=255)
    ggAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTACTACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAg                                                                                                <  1:1452497/100‑1 (MQ=255)
    ggAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAg                                                                                                >  1:3277354/1‑100 (MQ=255)
    ggAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAg                                                                                                >  1:2077620/1‑100 (MQ=255)
    ggAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAg                                                                                                <  1:372150/100‑1 (MQ=255)
    ggAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGACTGAATTTCAGATCACAg                                                                                                <  1:560647/100‑1 (MQ=255)
     gAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAgg                                                                                               <  1:3503998/100‑1 (MQ=255)
        ggTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATc                                                                                            <  1:2017251/100‑1 (MQ=255)
         gTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATcc                                                                                           <  2:4150097/100‑1 (MQ=255)
          tCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCca                                                                                          <  2:3898771/100‑1 (MQ=255)
           cgcgCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATcc                                                                                           <  2:23399/98‑1 (MQ=255)
           cgcgCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCcac                                                                                         <  1:160577/100‑1 (MQ=255)
           cgcgCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCcac                                                                                         >  1:3094666/1‑100 (MQ=255)
             cgcCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAGTTTCAGAGCACAGGATCcaccg                                                                                       >  1:3174714/1‑98 (MQ=255)
             cgcCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCacac                                                                                       >  1:1368931/1‑100 (MQ=255)
             cgcCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCacac                                                                                       >  2:2560044/1‑100 (MQ=255)
              gcCGATAATCATTTTCGACACCCAGTTCGTGTTGGTCGAACTGGTTATACCGAAATCGCGTCTTTTCAGCGGTGAAATGTCAGATCCCAGGATCCAcacc                                                                                      <  2:3399632/100‑1 (MQ=255)
                   tAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTgg                                                                                 >  1:3746269/1‑100 (MQ=255)
                        atatTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTggcgcgg                                                                            <  2:3286904/100‑1 (MQ=255)
                        atatTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTggcgcgg                                                                            <  2:1338649/100‑1 (MQ=255)
                          atTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTggcgcggcg                                                                          <  1:2431913/100‑1 (MQ=255)
                          atTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTggcgcggcg                                                                          >  2:3916731/1‑100 (MQ=255)
                           ttCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTggcgcggcgc                                                                         >  2:643032/1‑100 (MQ=255)
                            tCGACACCCAGTTCGTGGTGGTCGAACAGTTACTACCGACAGCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTggcgcggcgc                                                                         <  1:2770505/99‑1 (MQ=38)
                              gACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCa                                                                      >  1:3617521/1‑100 (MQ=255)
                              gACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCa                                                                      >  2:1397256/1‑100 (MQ=255)
                                 acCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCACTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCg                                                                   <  2:1439055/100‑1 (MQ=255)
                                  cccAGTTCGTGGTGGTCGAACAGGTAATACCGCAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGGTCCACACCACTGGCGCGGCGCTCATCGc                                                                  >  2:1598601/1‑100 (MQ=255)
                                  cccAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGTGCGGCGCTCATCGc                                                                  <  2:938116/100‑1 (MQ=255)
                                  cccAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGc                                                                  <  2:1789285/100‑1 (MQ=255)
                                   ccAGTTCGTGGTGGTCGAACAGGTAATACCGAACTCGCGTCATTTCAGCGGCTGAATTTCAGAACAGAGGATCCACACCACTGGCGCGGCGCTCATCGc                                                                  >  2:2775240/1‑99 (MQ=39)
                                   ccAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCa                                                                 <  1:3561752/100‑1 (MQ=255)
                                   ccAGTTCGTGGTGGTCGAACAGGTAATACCAAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCa                                                                 >  2:2774098/1‑100 (MQ=255)
                                    cAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCAt                                                                <  1:3821421/100‑1 (MQ=255)
                                     aGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATa                                                               <  2:2181255/100‑1 (MQ=255)
                                         cGTGGTGGTCGTACAGGTAATACAGAAATCGCGGCATTTCAGCGGCTGAATTTTCAATCACAGGATCCACACCACTGGCGCGTCGCTCATCTCATAAc                                                             >  2:1054521/1‑98 (MQ=37)
                                          gtggtgGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCgg                                                          <  1:2631735/100‑1 (MQ=255)
                                          gtggtgGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCgg                                                          <  1:643036/100‑1 (MQ=255)
                                           tggtggTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGc                                                         >  1:4096315/1‑100 (MQ=255)
                                           tggtggTCGAACAGGTAATACCGAAATCGCGGCAGTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGc                                                         <  2:1281957/100‑1 (MQ=255)
                                            ggtggtCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCa                                                        >  2:2875404/1‑100 (MQ=255)
                                            ggtggtCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCAGCGCATAACCGGc                                                         <  2:3859898/99‑1 (MQ=255)
                                              tggtCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAAt                                                      <  1:1550104/100‑1 (MQ=255)
                                              tggtCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAAt                                                      <  2:2683911/100‑1 (MQ=255)
                                               ggtCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATa                                                     >  1:119693/1‑100 (MQ=255)
                                               ggtCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATa                                                     >  1:2683927/1‑100 (MQ=255)
                                                gtCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATac                                                    <  1:151116/100‑1 (MQ=255)
                                                 tCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATaca                                                   <  2:2077609/100‑1 (MQ=255)
                                                 tCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATaca                                                   >  2:2837820/1‑100 (MQ=255)
                                                  cGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGTTGAATCTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATacac                                                  <  2:308546/100‑1 (MQ=255)
                                                  cGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATacac                                                  >  2:302073/1‑100 (MQ=255)
                                                   gAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATAcaca                                                 <  2:4126117/100‑1 (MQ=255)
                                                   gAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATAcaca                                                 <  2:1905958/100‑1 (MQ=255)
                                                    aaCAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACAcat                                                >  2:606786/1‑100 (MQ=255)
                                                     aCAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACAcatc                                               >  2:452205/1‑100 (MQ=255)
                                                     aCAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACAcatc                                               <  1:1859205/100‑1 (MQ=255)
                                                       aGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAAATTCAGATCACATGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACAcatcat                                             >  2:1446864/1‑100 (MQ=255)
                                                        ggTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCACATAACCGGCAATACACatcatc                                            >  1:3192342/1‑100 (MQ=255)
                                                         gTAATACAGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCg                                           <  2:1708052/100‑1 (MQ=255)
                                                              aCCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATg                                      >  1:528746/1‑100 (MQ=255)
                                                              aCCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATg                                      <  2:2012750/100‑1 (MQ=255)
                                                              aCCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATg                                      >  1:752086/1‑100 (MQ=255)
                                                              aCCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATg                                      >  1:3375726/1‑100 (MQ=255)
                                                               ccGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGc                                     <  1:3399697/100‑1 (MQ=255)
                                                               ccGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGc                                     >  1:33511/1‑100 (MQ=255)
                                                                cGAAATCGCGGCTTTTCAGCGGCTGAACTGCAGATCCCAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGcc                                    <  2:1785105/100‑1 (MQ=255)
                                                                cGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGcc                                    <  2:235009/100‑1 (MQ=255)
                                                                cGAAATCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGcc                                    <  1:808235/100‑1 (MQ=255)
                                                                cGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGcc                                    <  2:1178488/100‑1 (MQ=255)
                                                                cGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGcc                                    <  1:210478/100‑1 (MQ=255)
                                                                 gAAATCGCGGCATTTCAGCGGCTTAATTTCAGATCACAGGATCCACACCACTGTCGCGGCGCTCATCGCATCACCGGCAATACACATCATCGTc                                         >  2:872147/1‑94 (MQ=255)
                                                                   aaTCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAAc                                 <  2:3471703/100‑1 (MQ=255)
                                                                   aaTCGCGGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCCTCATGCCaa                                  >  2:3554196/1‑99 (MQ=255)
                                                                      cgcgGCCTTTCAGCGGCTGAATGTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACag                               <  2:2886880/99‑1 (MQ=255)
                                                                       gcgGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACagag                             >  1:2145950/1‑100 (MQ=255)
                                                                       gcgGCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACagag                             <  2:518217/100‑1 (MQ=255)
                                                                        cgGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAGCCGGCAATACACATCATCGTCATGCCAACagaga                            >  2:1095373/1‑100 (MQ=255)
                                                                         ggCATTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACagagag                           <  1:2738435/100‑1 (MQ=255)
                                                                           cATTTCAGCAGCTCAATTTCAGAGCACAGGACCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACagaga                            <  2:1638782/97‑1 (MQ=255)
                                                                            aTTTCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCa                        <  2:3617622/100‑1 (MQ=255)
                                                                             tttCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAg                       <  1:1002432/100‑1 (MQ=255)
                                                                              ttCAGCGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAgg                      >  1:1151174/1‑100 (MQ=255)
                                                                                   cGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAg                 >  1:1064274/1‑100 (MQ=255)
                                                                                   cGGCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAg                 <  1:2560059/100‑1 (MQ=255)
                                                                                    ggCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAAAGCCAGGACAAGc                >  1:446670/1‑100 (MQ=255)
                                                                                     gCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACACGcc               >  2:3633326/1‑100 (MQ=255)
                                                                                     gCTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGcc               >  1:3973563/1‑100 (MQ=255)
                                                                                      cTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGcc               >  1:3305732/1‑99 (MQ=255)
                                                                                      cTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGcc               >  2:3079048/1‑99 (MQ=255)
                                                                                      cTGAATTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCa              >  1:1271855/1‑100 (MQ=255)
                                                                                         aaTTTCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTg           >  2:3375586/1‑100 (MQ=255)
                                                                                          aTTTCAGATCACAGGATCCACACCGCTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTgg          <  2:2602007/100‑1 (MQ=255)
                                                                                          aTTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTgg          <  1:2837836/100‑1 (MQ=255)
                                                                                            ttCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTGGaa        <  1:1404767/100‑1 (MQ=255)
                                                                                            ttCAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTGGaa        <  1:2303315/100‑1 (MQ=255)
                                                                                             tCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAg       <  1:1993061/100‑1 (MQ=255)
                                                                                              cAGATCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAGc      >  2:709562/1‑100 (MQ=255)
                                                                                               agatCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAGcc     >  2:3521004/1‑100 (MQ=255)
                                                                                                 agCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGt   >  1:1026748/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GTCAGGAAGGTCGCGCCGATAATCATATTCGACACCCAGTTCGTGGTGGTCGAACAGGTAATACCGAAATCGCGGCATTTCAGCGGCTGAATTTCAGAGCACAGGATCCACACCACTGGCGCGGCGCTCATCGCATAACCGGCAATACACATCATCGTCATGCCAACAGAGAGCCAGGACAAGCCACTGGAAGCCGTA  >  NC_000913/2980933‑2981130

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