Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,348,017 C→T A408A (GCG→GCA yjdH ← putative 2‑component sensor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,348,0170CT98.8% 838.8 / ‑4.1 243A408A (GCG→GCAyjdHputative 2‑component sensor protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (124/116);  total (125/118)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

Only 100 of 263 total aligned reads displayed.

ACGTGGTAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTCGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATT  >  NC_000913/4347926‑4348106
                                                                                           |                                                                                         
 cGTGGTAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTtgc                                                                                         <  2:4312691/93‑1 (MQ=255)
 cGTGGTAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTtgc                                                                                         <  1:4405407/93‑1 (MQ=255)
  gTGGTAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTtgcg                                                                                        >  2:434401/1‑93 (MQ=255)
  gTGGTAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTtgcg                                                                                        >  2:3522594/1‑93 (MQ=255)
   tGGTAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTtgcgc                                                                                       <  2:4445208/93‑1 (MQ=255)
     gTAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCgg                                                                                     >  1:1400757/1‑93 (MQ=255)
      tAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGt                                                                                    <  2:3458860/93‑1 (MQ=255)
         tCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTaa                                                                                 <  2:4204838/93‑1 (MQ=255)
         tCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTaa                                                                                 <  1:668346/93‑1 (MQ=255)
         tCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTaa                                                                                 >  1:3794520/1‑93 (MQ=255)
          cAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAAt                                                                                >  1:111107/1‑93 (MQ=255)
          cAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAAt                                                                                <  1:87028/93‑1 (MQ=255)
           aGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATc                                                                               <  2:2057741/93‑1 (MQ=255)
           aGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATc                                                                               >  2:518862/1‑93 (MQ=255)
           aGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGCTAAtt                                                                               >  1:45649/1‑92 (MQ=255)
            gCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCt                                                                              >  1:712109/1‑93 (MQ=255)
            gCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTCGCGCGGTTAATCt                                                                              <  2:3929019/93‑1 (MQ=255)
            gCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGGTGCGCGGTTAATCt                                                                              >  2:942745/1‑93 (MQ=255)
             cGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCtt                                                                             <  2:2662188/93‑1 (MQ=255)
             cGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCtt                                                                             <  2:1292475/93‑1 (MQ=255)
             cGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCtt                                                                             >  1:137897/1‑93 (MQ=255)
              gTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTg                                                                            >  2:500385/1‑93 (MQ=255)
               tCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGc                                                                           >  2:2364073/1‑93 (MQ=255)
               tCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGc                                                                           <  2:1347773/93‑1 (MQ=255)
               tCGCGACCTGGTCCTCACTGCCGCTGTCAGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGc                                                                           >  1:3748283/1‑93 (MQ=255)
                cgcgACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCt                                                                          <  2:3794520/93‑1 (MQ=255)
                cgcgACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCt                                                                          <  1:468384/93‑1 (MQ=255)
                   gACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGAt                                                                       <  1:108465/93‑1 (MQ=255)
                    gccTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGAtt                                                                      <  2:4103170/92‑1 (MQ=255)
                    aCCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGAtt                                                                      <  2:252204/93‑1 (MQ=255)
                     ccTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTa                                                                     >  2:3497736/1‑93 (MQ=255)
                     ccTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTa                                                                     >  1:2867689/1‑93 (MQ=255)
                      cTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTaa                                                                    >  2:1220574/1‑93 (MQ=255)
                       tGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTaaa                                                                   <  2:4518744/93‑1 (MQ=255)
                       tGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAAGCGGTTGCGCGGTTAATCTTGCTGATTaaa                                                                   <  2:1590890/93‑1 (MQ=255)
                        ggTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTaaaa                                                                  >  2:1614308/1‑93 (MQ=255)
                        ggTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTaaaa                                                                  >  2:3858455/1‑93 (MQ=255)
                           ccTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAAc                                                               >  2:3249633/1‑93 (MQ=255)
                             tCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCa                                                             <  2:1108206/93‑1 (MQ=255)
                             tCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCa                                                             >  1:3499155/1‑93 (MQ=255)
                             tCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCa                                                             >  1:2979105/1‑93 (MQ=255)
                               aCTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGc                                                           <  1:551082/93‑1 (MQ=255)
                                cTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCg                                                          >  1:2214741/1‑93 (MQ=255)
                                 tGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGa                                                         <  2:2287030/93‑1 (MQ=255)
                                 tGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGa                                                         <  2:2451788/93‑1 (MQ=255)
                                  gCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGAt                                                        <  1:939724/93‑1 (MQ=255)
                                    cGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATaa                                                      >  2:2894639/1‑93 (MQ=255)
                                     gCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAAc                                                     >  1:3009507/1‑93 (MQ=255)
                                     gCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAAc                                                     <  1:2817886/93‑1 (MQ=255)
                                       tGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCg                                                   <  2:17467/93‑1 (MQ=255)
                                       tGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCg                                                   <  1:1700168/93‑1 (MQ=255)
                                        gTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCgg                                                  <  1:3634353/93‑1 (MQ=255)
                                         tCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGga                                                 <  2:1216000/93‑1 (MQ=255)
                                         tCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGga                                                 <  1:1806014/93‑1 (MQ=255)
                                             gCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGAtt                                             <  1:2763998/93‑1 (MQ=255)
                                              cAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGAttt                                            <  2:1124379/93‑1 (MQ=255)
                                               aGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTg                                           >  2:794693/1‑93 (MQ=255)
                                                 cTGGCTTTCACTGTTTAAAATCATCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGAt                                         <  2:2591552/93‑1 (MQ=255)
                                                  tGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATc                                        >  2:4661527/1‑93 (MQ=255)
                                                  tGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGAGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATc                                        >  2:1449028/1‑93 (MQ=255)
                                                   ggCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCt                                       >  1:1223184/1‑93 (MQ=255)
                                                    gCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCtt                                      <  1:1971864/93‑1 (MQ=255)
                                                    gCTTTCAATGTTTAAAATCAGCGGATGCCGTAAATCGGTTGCGGGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCtt                                      <  1:3520825/93‑1 (MQ=255)
                                                      tttCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTAc                                    <  1:1934046/93‑1 (MQ=255)
                                                      tttCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTAc                                    >  1:3072504/1‑93 (MQ=255)
                                                       ttCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTAcc                                   >  1:1298375/1‑93 (MQ=255)
                                                         cACTGTTTAAAATCAGCGTATGGCCTAAATCGGTCGCGCGGTTAATCGTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTTCCca                                 <  2:4020566/93‑1 (MQ=255)
                                                             gTTTAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCcagcag                             <  1:2553457/93‑1 (MQ=255)
                                                              tttAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCaga                            <  1:618763/93‑1 (MQ=255)
                                                                tAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCagaga                          <  2:1949595/93‑1 (MQ=255)
                                                                tAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCagaga                          <  2:3732596/93‑1 (MQ=255)
                                                                tAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCagaga                          <  2:3573245/93‑1 (MQ=255)
                                                                tAAAATCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCagaga                          >  1:2309442/1‑93 (MQ=255)
                                                                   aaTCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGcc                       <  2:3754755/93‑1 (MQ=255)
                                                                   aaTCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGcc                       >  1:2068245/1‑93 (MQ=255)
                                                                     tCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCaa                     >  1:3611553/1‑93 (MQ=255)
                                                                     tCAGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCaa                     <  2:3398612/93‑1 (MQ=255)
                                                                       aGCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAAtc                   >  2:71816/1‑93 (MQ=255)
                                                                        gCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAAtct                  >  1:1199715/1‑93 (MQ=255)
                                                                        gCGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAAtct                  >  1:751827/1‑93 (MQ=255)
                                                                         cGTATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAAtctc                 <  1:3998418/93‑1 (MQ=255)
                                                                           tATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGGGATTTGATCTTACCCAGCATAAAGCCAATGTCtt               >  2:2310933/1‑93 (MQ=255)
                                                                           tATGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCtt               >  1:208947/1‑93 (MQ=255)
                                                                            aTGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTc              <  1:268008/93‑1 (MQ=255)
                                                                             tGGCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTcc             <  1:258752/93‑1 (MQ=255)
                                                                              ggCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCt            <  1:4487123/93‑1 (MQ=255)
                                                                              ggCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCt            <  2:1020877/93‑1 (MQ=255)
                                                                               gCCTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTg           >  1:3870420/1‑93 (MQ=255)
                                                                                ccTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGa          <  1:4616195/93‑1 (MQ=255)
                                                                                 cTAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCAAGCAGAGAGCCAATCTCTTCCTGAt         >  1:387701/1‑93 (MQ=255)
                                                                                  tAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATa        <  1:205744/93‑1 (MQ=255)
                                                                                  tAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATa        <  1:3237414/93‑1 (MQ=255)
                                                                                  tAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATa        <  1:1547334/93‑1 (MQ=255)
                                                                                  tAAATCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATa        <  1:2040681/93‑1 (MQ=255)
                                                                                  tAAATCGGTTGCGCGGGTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATa        >  1:4059584/1‑93 (MQ=255)
                                                                                    aaTCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGt      <  1:971441/93‑1 (MQ=255)
                                                                                      tCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTa    >  1:3282036/1‑93 (MQ=255)
                                                                                      tCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTa    >  1:4719078/1‑93 (MQ=255)
                                                                                      tCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTa    <  2:665344/93‑1 (MQ=255)
                                                                                      tCGGTTGCGCGGTTAATCTTGCTGATTAAAAAACAAGCGATAACCGGGGATTTGCTCTAACCCAGCAAAAAGCAAACCTCTCCCTCATGGttt    >  2:1265465/1‑92 (MQ=255)
                                                                                           |                                                                                         
ACGTGGTAATCAGCGTCGCGACCTGGTCCTCACTGCCGCTGTCTGGCAGCTGGCTTTCACTGTTTAAAATCAGCGTATGGCCTAAATCGGTCGCGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATT  >  NC_000913/4347926‑4348106

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