Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,739,740 A→G M141T (ATG→ACG)  yiaJ ← putative regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,739,7400AG99.6% 791.0 / ‑5.1 248M141T (ATG→ACG) yiaJputative regulator
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (111/136);  total (112/136)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 265 total aligned reads displayed.

CGTAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTC  >  NC_000913/3739648‑3739829
                                                                                            |                                                                                         
  tAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTc                                                                                         >  1:929615/1‑93 (MQ=255)
   aGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTcc                                                                                        <  2:2302045/93‑1 (MQ=255)
    gTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTccc                                                                                       >  1:2101982/1‑93 (MQ=255)
       cGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGt                                                                                    >  1:4464097/1‑93 (MQ=255)
       cGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGt                                                                                    <  1:5217055/93‑1 (MQ=255)
        gggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGtt                                                                                   <  1:5923572/93‑1 (MQ=255)
         ggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtg                                                                                  >  2:1468169/1‑93 (MQ=255)
          gTGACTAAACGCCATGCAGATCTTGTCTGTTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATATGCACGGGTTCGCAGCGTCCCGGTtgt                                                                                 <  2:3562083/93‑1 (MQ=255)
          gTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgt                                                                                 <  1:2900503/93‑1 (MQ=255)
           ttaCCAAACGCATTGTAGACCTTGCCCATTGCGGAGCAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                <  2:3907640/91‑1 (MQ=255)
           tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                >  2:3419366/1‑93 (MQ=255)
           tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                >  2:2972803/1‑93 (MQ=255)
           tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                >  1:396486/1‑93 (MQ=255)
           tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                >  1:2794125/1‑93 (MQ=255)
           tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                <  1:4805569/93‑1 (MQ=255)
           tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                <  2:2127625/93‑1 (MQ=255)
            gACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTgg                                                                               <  2:1181996/93‑1 (MQ=255)
             aCCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTggg                                                                              >  1:4188578/1‑93 (MQ=255)
             aCCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTggg                                                                              <  2:2866491/93‑1 (MQ=255)
               cAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGtt                                                                            <  1:1116911/93‑1 (MQ=255)
                aaaCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTc                                                                           >  2:5353912/1‑93 (MQ=255)
                 aaCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTcc                                                                          >  1:2516809/1‑93 (MQ=255)
                 aaCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTcc                                                                          <  1:258398/93‑1 (MQ=255)
                    gCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGc                                                                       <  1:1345424/93‑1 (MQ=255)
                    gCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGc                                                                       <  2:5380327/93‑1 (MQ=255)
                     ccATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCt                                                                      >  1:3712331/1‑93 (MQ=255)
                       aTGTAGATCTTGCCCATTGCGGATCAGTAGAGCGGCATATGCGGGCCAATATAGGCACGGGGTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    >  1:6301894/1‑93 (MQ=255)
                       aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCTGGTTGTGGGTTCCAGCTta                                                                    <  1:5674406/93‑1 (MQ=255)
                       aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    >  1:2476253/1‑93 (MQ=255)
                        tGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtat                                                                   >  2:3083269/1‑93 (MQ=255)
                        tGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtat                                                                   >  2:5327222/1‑93 (MQ=255)
                         gTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtata                                                                  <  1:1039515/93‑1 (MQ=255)
                         gTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtata                                                                  <  1:4258776/93‑1 (MQ=255)
                          tAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaa                                                                 <  2:3870432/93‑1 (MQ=255)
                          tAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaa                                                                 <  2:2444901/93‑1 (MQ=255)
                          tAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaa                                                                 <  1:5287168/93‑1 (MQ=255)
                           aGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaaa                                                                <  2:1741461/93‑1 (MQ=255)
                           aGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaaa                                                                <  1:1210791/93‑1 (MQ=255)
                             aTCTTGCCCATTGCGGAACAGTAGAGCGGCATATGGTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCCGCTTATAAATc                                                              >  2:4835654/1‑93 (MQ=255)
                             aTCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCCGCTTATAAATc                                                              >  2:1393688/1‑93 (MQ=255)
                             aTCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGATCGCAGCGTCCCGGGTGTGGGTTCCAGCATATAAATc                                                              >  2:4000455/1‑93 (MQ=255)
                              tCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCa                                                             >  1:2136512/1‑93 (MQ=255)
                              tCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCa                                                             <  1:5773789/93‑1 (MQ=255)
                               cTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaa                                                            >  1:2323447/1‑93 (MQ=255)
                               cTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaa                                                            <  2:5004496/93‑1 (MQ=255)
                                ttGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaaa                                                           >  2:4594949/1‑93 (MQ=255)
                                ttGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaaa                                                           <  1:5041210/93‑1 (MQ=255)
                                  gCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAAt                                                         <  2:2361103/93‑1 (MQ=255)
                                  gCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAAt                                                         <  1:1021828/93‑1 (MQ=255)
                                    ccATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAg                                                       <  1:3801356/93‑1 (MQ=255)
                                     cATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      <  1:2982006/93‑1 (MQ=255)
                                      aTTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCg                                                     <  1:2304058/93‑1 (MQ=255)
                                       ttGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGt                                                    >  1:989097/1‑93 (MQ=255)
                                          cGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGAt                                                 <  2:1601965/93‑1 (MQ=255)
                                          cGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCAGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGAt                                                 <  2:4726949/93‑1 (MQ=255)
                                           ggAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATc                                                >  2:3640278/1‑93 (MQ=255)
                                           ggAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATc                                                >  2:3741450/1‑93 (MQ=255)
                                            gAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCg                                               >  1:2055844/1‑93 (MQ=255)
                                            gAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCg                                               >  2:1975745/1‑93 (MQ=255)
                                             aaCAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGt                                              <  2:2218258/93‑1 (MQ=255)
                                              aCAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTc                                             >  1:3336959/1‑93 (MQ=255)
                                                aGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCtt                                           <  2:2686131/93‑1 (MQ=255)
                                                 gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          <  1:4816698/93‑1 (MQ=255)
                                                 gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          <  2:1648909/93‑1 (MQ=255)
                                                   agagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcgc                                        <  2:4294969/93‑1 (MQ=255)
                                                        ggCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTg                                   <  2:5868923/93‑1 (MQ=255)
                                                          cATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGga                                 <  2:1758066/93‑1 (MQ=255)
                                                          cATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGATggg                                 >  1:1495317/1‑92 (MQ=255)
                                                            tatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgaga                               >  1:4032242/1‑93 (MQ=255)
                                                            tatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgaga                               <  1:2711150/93‑1 (MQ=255)
                                                            tatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgaga                               >  1:1431925/1‑93 (MQ=255)
                                                             atGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGaa                              >  2:6019192/1‑93 (MQ=255)
                                                              tGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAg                             >  2:4431826/1‑93 (MQ=255)
                                                                 tGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTa                          >  1:1656123/1‑93 (MQ=255)
                                                                  ggCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTaa                         >  2:3819622/1‑93 (MQ=255)
                                                                   gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        <  1:75741/93‑1 (MQ=255)
                                                                    ccAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATg                       <  1:4355636/93‑1 (MQ=255)
                                                                    ccAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATg                       >  1:2203743/1‑93 (MQ=255)
                                                                     cAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGTAGTTAATgg                      >  1:4563292/1‑93 (MQ=255)
                                                                     cAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATgg                      >  1:6081928/1‑93 (MQ=255)
                                                                      aaTATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGt                     <  2:2842863/93‑1 (MQ=255)
                                                                      aaTATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGt                     >  1:3476099/1‑93 (MQ=255)
                                                                         ataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTTATGGTTTc                  <  2:6124934/93‑1 (MQ=255)
                                                                         ataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTc                  <  2:6286770/93‑1 (MQ=255)
                                                                         ataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTc                  >  1:6130196/1‑93 (MQ=255)
                                                                          taGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCa                 <  1:1746738/93‑1 (MQ=255)
                                                                          taGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCa                 >  1:5611674/1‑93 (MQ=255)
                                                                            ggCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCAcc               <  1:5398021/93‑1 (MQ=255)
                                                                             gCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCa              >  1:2281668/1‑93 (MQ=255)
                                                                             gCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCa              <  2:2791382/93‑1 (MQ=255)
                                                                              cACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAg             <  2:1837833/93‑1 (MQ=255)
                                                                               gcGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGt            <  2:4493680/92‑1 (MQ=255)
                                                                                cGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTg           <  1:1590210/93‑1 (MQ=255)
                                                                                 gggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTgg          <  1:5903315/93‑1 (MQ=255)
                                                                                 gggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTgg          >  2:5721863/1‑93 (MQ=255)
                                                                                 gggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTgg          >  1:4112061/1‑93 (MQ=255)
                                                                                     tCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGAt      <  2:4829528/93‑1 (MQ=255)
                                                                                       gcagcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGt    <  2:1269273/93‑1 (MQ=255)
                                                                                       gcagcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGt    <  2:707925/93‑1 (MQ=255)
                                                                                        cagcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGtt   <  1:473711/93‑1 (MQ=255)
                                                                                            |                                                                                         
CGTAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTC  >  NC_000913/3739648‑3739829

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