Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,094,047 T→C V61A (GTA→GCA)  hisF → imidazole glycerol phosphate synthase holoenzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,094,0470TC83.5% 585.7 / 102.8 260V61A (GTA→GCA) hisFimidazole glycerol phosphate synthase holoenzyme
Reads supporting (aligned to +/- strand):  ref base T (20/23);  new base C (104/113);  total (124/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 261 total aligned reads displayed.

GATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAA  >  NC_000913/2093956‑2094139
                                                                                           |                                                                                            
gATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTa                                                                                             <  2:1084213/93‑1 (MQ=255)
  tatCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGa                                                                                           <  1:6825863/93‑1 (MQ=255)
   atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGAt                                                                                          <  2:5232336/93‑1 (MQ=255)
   atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGAt                                                                                          >  2:2098433/1‑93 (MQ=255)
    tCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATa                                                                                         >  2:6717266/1‑93 (MQ=255)
      gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATaaa                                                                                       >  2:7085877/1‑93 (MQ=255)
       tGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATaaaa                                                                                      >  1:6519060/1‑93 (MQ=255)
         ccGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGc                                                                                    >  2:6392527/1‑93 (MQ=255)
           gCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTg                                                                                  >  2:6306438/1‑93 (MQ=255)
           gCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTg                                                                                  <  2:2570755/93‑1 (MQ=255)
            cTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTgg                                                                                 <  1:6047637/93‑1 (MQ=255)
            cTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTgg                                                                                 <  1:4363218/93‑1 (MQ=255)
             tGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTggg                                                                                <  1:5750735/93‑1 (MQ=255)
             tGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTggg                                                                                <  2:4864109/93‑1 (MQ=255)
             tGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTggg                                                                                >  1:2772904/1‑93 (MQ=255)
               gCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTa                                                                              <  2:5497107/93‑1 (MQ=255)
               gCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTa                                                                              <  2:6078823/93‑1 (MQ=255)
                 aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTAtc                                                                            >  1:591052/1‑93 (MQ=255)
                 aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTAtc                                                                            <  1:4747944/93‑1 (MQ=255)
                 aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGCGTGGCAGATAAAAGCTGGGTAtc                                                                            >  1:2719744/1‑93 (MQ=255)
                   aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTAtctc                                                                          >  2:2825711/1‑93 (MQ=255)
                   aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTAtctc                                                                          >  1:2903890/1‑93 (MQ=255)
                     cgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTcgc                                                                        >  2:5320001/1‑93 (MQ=255)
                     cgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTcgc                                                                        >  1:2389497/1‑93 (MQ=255)
                      gctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTcgcg                                                                       >  2:4449500/1‑93 (MQ=255)
                      gctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTcgcg                                                                       <  1:6278116/93‑1 (MQ=255)
                       ctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGt                                                                      <  2:3956633/93‑1 (MQ=255)
                        tacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTcgcgcg                                                                     >  1:45080/1‑91 (MQ=255)
                        tacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTg                                                                     <  2:2355910/93‑1 (MQ=255)
                          cgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGc                                                                   <  2:1298808/93‑1 (MQ=255)
                          cgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGc                                                                   <  2:2713569/93‑1 (MQ=255)
                          cgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGc                                                                   <  2:3414544/93‑1 (MQ=255)
                           gctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCg                                                                  <  1:3952601/93‑1 (MQ=255)
                            ctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCgg                                                                 >  1:3743143/1‑93 (MQ=255)
                             tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGa                                                                >  1:830439/1‑93 (MQ=255)
                               aagaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAg                                                              >  1:7153660/1‑93 (MQ=255)
                               aagaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAg                                                              <  2:4033238/93‑1 (MQ=255)
                                 gaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTg                                                            >  1:6337943/1‑93 (MQ=255)
                                 gaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTg                                                            >  2:6822467/1‑93 (MQ=255)
                                  aagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGa                                                           <  1:7439723/93‑1 (MQ=255)
                                  aagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGa                                                           <  2:1919006/93‑1 (MQ=255)
                                  aagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGa                                                           <  1:5021466/93‑1 (MQ=255)
                                   agGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGAt                                                          >  1:3308530/1‑93 (MQ=255)
                                    ggCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATc                                                         <  2:5079230/93‑1 (MQ=255)
                                        cTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACa                                                     <  2:4641464/93‑1 (MQ=255)
                                        cTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACa                                                     <  1:993254/93‑1 (MQ=255)
                                         tGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACAt                                                    <  1:2209808/93‑1 (MQ=255)
                                            cGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTcc                                                 <  2:5002841/93‑1 (MQ=255)
                                            cGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTcc                                                 <  2:2913474/93‑1 (MQ=255)
                                             gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCg                                                <  1:1982926/93‑1 (MQ=255)
                                             gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCg                                                >  2:6530032/1‑93 (MQ=255)
                                              aaCTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGt                                               >  2:5863364/1‑93 (MQ=255)
                                                cTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGttt                                             <  2:296215/93‑1 (MQ=255)
                                                cTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGttt                                             <  1:3690837/93‑1 (MQ=255)
                                                  ggTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtg                                           >  1:5155180/1‑93 (MQ=255)
                                                   gtgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtgt                                          >  1:1024074/1‑93 (MQ=255)
                                                   gtgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtgt                                          <  1:4928583/93‑1 (MQ=255)
                                                    tgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtgtg                                         >  2:1250677/1‑93 (MQ=255)
                                                    tgtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtgtg                                         <  2:426825/93‑1 (MQ=255)
                                                     gtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtgtgt                                        <  1:1609758/93‑1 (MQ=255)
                                                     gtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtgtgt                                        <  2:3014766/93‑1 (MQ=255)
                                                      ttCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTtgtgtg                                       <  1:6425870/93‑1 (MQ=255)
                                                       tCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTgg                                      >  2:1457103/1‑93 (MQ=255)
                                                          aCGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCgg                                   <  1:7234248/93‑1 (MQ=255)
                                                           cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCggg                                  <  1:4227617/93‑1 (MQ=255)
                                                           cGATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCggg                                  <  2:3291788/93‑1 (MQ=255)
                                                            gATATCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGt                                 <  1:3665180/93‑1 (MQ=255)
                                                            gATATCACCGCTTCCAGCGATGGCCGCGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGt                                 <  2:4695002/93‑1 (MQ=255)
                                                              tatCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTgg                               >  2:2632748/1‑93 (MQ=255)
                                                              tatCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTgg                               <  1:3764852/93‑1 (MQ=255)
                                                               atCACCGCTTCCAGGGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTggg                              >  2:5188992/1‑93 (MQ=255)
                                                                tCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGa                             <  2:3544952/93‑1 (MQ=255)
                                                                tCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGa                             <  2:1740923/93‑1 (MQ=255)
                                                                tCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGa                             <  1:1684116/93‑1 (MQ=255)
                                                                tCACCGCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAACTGATCGACATTCCGTTTTGTGTGGCGGGTGGGa                             >  2:7362648/1‑93 (MQ=255)
                                                                 cACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGgggggggt                            >  1:6110466/1‑91 (MQ=255)
                                                                     gCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAg                        >  1:7315663/1‑93 (MQ=255)
                                                                     gCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAg                        >  1:5198156/1‑93 (MQ=255)
                                                                     gCTTCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAg                        >  2:6901624/1‑93 (MQ=255)
                                                                       ttCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGtc                      <  2:7510124/93‑1 (MQ=255)
                                                                       ttCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAggc                      >  1:242236/1‑91 (MQ=255)
                                                                       ttCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAAGtc                      >  1:7505940/1‑93 (MQ=255)
                                                                        tCCAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGtct                     >  1:116119/1‑93 (MQ=255)
                                                                          cAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTctct                   >  2:7392094/1‑93 (MQ=255)
                                                                          cAGCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTctct                   <  1:3518630/93‑1 (MQ=255)
                                                                            gCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTgg                 <  2:5439077/93‑1 (MQ=255)
                                                                            gCGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGGGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTgg                 >  1:344805/1‑93 (MQ=255)
                                                                             cGATGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGa                <  2:135636/93‑1 (MQ=255)
                                                                              gATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGaa               <  2:4508357/93‑1 (MQ=255)
                                                                               aTGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAg              >  2:696657/1‑93 (MQ=255)
                                                                                tGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGa             <  2:215439/93‑1 (MQ=255)
                                                                                tGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGa             <  1:4144282/93‑1 (MQ=255)
                                                                                tGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGa             <  1:2552345/93‑1 (MQ=255)
                                                                                tGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGa             <  1:5782324/93‑1 (MQ=255)
                                                                                tGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGa             >  2:1509243/1‑93 (MQ=255)
                                                                                tGGCCGTGTGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGa             <  2:1793418/93‑1 (MQ=255)
                                                                                     gtgtGGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCcg        >  2:570019/1‑93 (MQ=255)
                                                                                     gtgtGACAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCcg        >  1:4699282/1‑93 (MQ=255)
                                                                                       gtgGCAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCcgcg      <  1:1533556/93‑1 (MQ=255)
                                                                                         ggTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGaa    >  1:3157900/1‑93 (MQ=255)
                                                                                           |                                                                                            
GATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAA  >  NC_000913/2093956‑2094139

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