Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,094,539 T→C V225A (GTT→GCT)  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,5390TC98.6% 453.0 / ‑2.8 144V225A (GTT→GCT) hisFimidazole glycerol phosphate synthase holoenzyme
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (76/66);  total (77/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

Only 100 of 152 total aligned reads displayed.

AAGTGCGTGAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCA  >  NC_000913/2094449‑2094631
                                                                                          |                                                                                            
 aGTGCGTGAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGa                                                                                           <  2:2170014/93‑1 (MQ=255)
  gTGCGTGAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGAc                                                                                          <  2:1837168/93‑1 (MQ=255)
   tGCGTGAAGTTTGCCACGTCTCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACg                                                                                         <  1:4048782/93‑1 (MQ=255)
     cGTGAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGgc                                                                                       <  1:3160002/93‑1 (MQ=255)
      gTGAAGTTTGCCACGTCCCGCTTATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGgcg                                                                                      >  1:2152856/1‑93 (MQ=255)
      gTGAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGgcg                                                                                      >  1:1271586/1‑93 (MQ=255)
      gTGAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGgcg                                                                                      >  1:1415889/1‑93 (MQ=255)
       tGAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGgcgc                                                                                     >  1:52202/1‑93 (MQ=255)
       tGAAGATTGCCACGTCCCGCTGATTTCCTCCGGTGGCGCGGGCACCATGAAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGgcgc                                                                                     <  2:4050681/93‑1 (MQ=255)
        gAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGgcgcg                                                                                    >  2:3648920/1‑93 (MQ=255)
          aGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCt                                                                                  >  2:2295180/1‑93 (MQ=255)
            tttGCCACTTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTgg                                                                                <  2:3272171/93‑1 (MQ=255)
            tttGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTgg                                                                                <  1:2736396/93‑1 (MQ=255)
             ttGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGc                                                                               <  1:2497340/93‑1 (MQ=255)
              tGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCa                                                                              <  1:3903184/93‑1 (MQ=255)
                ccACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGc                                                                            <  2:1935438/93‑1 (MQ=255)
                  aCGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCtt                                                                          >  2:561027/1‑93 (MQ=255)
                   cgccccGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTc                                                                         <  2:1974727/90‑1 (MQ=255)
                   cGTCTCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTc                                                                         <  1:1736458/93‑1 (MQ=255)
                   cGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTc                                                                         <  2:2089812/93‑1 (MQ=255)
                   cGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTc                                                                         >  1:1236663/1‑93 (MQ=255)
                      cccGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGt                                                                      >  1:2116406/1‑93 (MQ=255)
                       ccGCTGATTGCCTCCGGTGTCGCGGGCCCCATGGAACTCTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTa                                                                     <  2:1483006/93‑1 (MQ=255)
                       ccGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTa                                                                     <  1:3599696/93‑1 (MQ=255)
                       ccGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTa                                                                     <  1:4059390/93‑1 (MQ=255)
                       ccGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTa                                                                     >  2:4261350/1‑93 (MQ=255)
                        cGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTAt                                                                    <  2:571188/93‑1 (MQ=255)
                        cGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTAt                                                                    >  2:2819637/1‑93 (MQ=255)
                        cGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTAt                                                                    <  1:161595/93‑1 (MQ=255)
                           tGATTGCCTCCGGTGGCGCGGTGACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTcc                                                                 >  1:850390/1‑93 (MQ=255)
                           tGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTcc                                                                 <  1:2819479/93‑1 (MQ=255)
                            gATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCca                                                                >  2:629153/1‑93 (MQ=255)
                            gATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCca                                                                >  2:2186077/1‑93 (MQ=255)
                            gATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCca                                                                >  1:2162284/1‑93 (MQ=255)
                            gATAGCCTCCGGTGTCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCca                                                                >  2:2360697/1‑93 (MQ=255)
                              ttGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCcaca                                                              >  1:2963519/1‑93 (MQ=255)
                                gCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCacaaa                                                            >  1:4369946/1‑93 (MQ=255)
                                gCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCacaaa                                                            <  2:2704044/93‑1 (MQ=255)
                                 ccTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGCCGGCGCGCTGGCAGCTTCCGTATTCCacaaac                                                           >  2:3440703/1‑93 (MQ=255)
                                  cTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCacaaaca                                                          >  2:1175476/1‑93 (MQ=255)
                                   tCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCacaaacaa                                                         >  1:4129206/1‑93 (MQ=255)
                                    ccGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCAcaaacaaa                                                        >  2:1067429/1‑93 (MQ=255)
                                      ggTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAaata                                                      <  1:920969/93‑1 (MQ=255)
                                       gTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAaataa                                                     >  1:4276897/1‑93 (MQ=255)
                                       gTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAaataa                                                     >  1:4554055/1‑93 (MQ=255)
                                       gTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAaataa                                                     >  2:220330/1‑93 (MQ=255)
                                        tGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAaataat                                                    <  1:1774627/93‑1 (MQ=255)
                                         ggCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATc                                                   >  1:2859369/1‑93 (MQ=255)
                                           cgcgGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCaa                                                 <  1:3723918/93‑1 (MQ=255)
                                            gcgGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAat                                                <  1:3349501/93‑1 (MQ=255)
                                            gcgGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAat                                                <  1:4649052/93‑1 (MQ=255)
                                            gcgGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCGTCCGTATTCCACAAACGAATAATCAat                                                >  1:1065517/1‑93 (MQ=255)
                                             cgGGCACCATGTAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAata                                               >  1:180266/1‑93 (MQ=255)
                                             cgGGCACCATGTAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAata                                               >  1:180267/1‑93 (MQ=255)
                                               ggCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATAtt                                             >  1:1006283/1‑93 (MQ=255)
                                                  cccATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGt                                          <  1:3721934/92‑1 (MQ=255)
                                                  aCCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGt                                          >  2:4417701/1‑93 (MQ=255)
                                                  aCCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGt                                          <  1:3224829/93‑1 (MQ=255)
                                                     aTGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGaa                                       >  1:3033872/1‑93 (MQ=255)
                                                      tGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACACATAATCAATATTGGTGAAt                                      >  1:1380212/1‑93 (MQ=255)
                                                      tGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAAt                                      >  1:4619311/1‑93 (MQ=255)
                                                      tGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAAt                                      <  1:2681917/93‑1 (MQ=255)
                                                      tGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAAt                                      >  2:1968123/1‑93 (MQ=255)
                                                      tGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAAt                                      <  1:369764/93‑1 (MQ=255)
                                                      tGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAAt                                      <  1:3246891/93‑1 (MQ=255)
                                                       ggAACACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAAtt                                     <  2:3104769/93‑1 (MQ=255)
                                                        gAACACTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTa                                    <  2:4522953/93‑1 (MQ=255)
                                                         aaCACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTaa                                   <  1:4509047/93‑1 (MQ=255)
                                                         aaCACTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTaa                                   >  1:3869212/1‑93 (MQ=255)
                                                           cacTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTaaaa                                 <  1:996643/93‑1 (MQ=255)
                                                           cacTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTaaaa                                 <  2:2627317/93‑1 (MQ=255)
                                                           cacTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTaaaa                                 <  2:1777729/93‑1 (MQ=255)
                                                             cTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGc                               <  2:486636/93‑1 (MQ=255)
                                                             cTTCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGc                               <  2:486637/93‑1 (MQ=255)
                                                               tCCTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGt                             <  1:1691667/93‑1 (MQ=255)
                                                                tcTCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTa                            <  1:3685259/92‑1 (MQ=255)
                                                                  tCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTAcc                          >  1:2794170/1‑93 (MQ=255)
                                                                  tCGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTAcc                          <  2:3289680/93‑1 (MQ=255)
                                                                   cGAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCt                         <  2:897674/93‑1 (MQ=255)
                                                                    gAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTg                        >  2:3005889/1‑93 (MQ=255)
                                                                    gAAGCCTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTg                        <  2:4097404/93‑1 (MQ=255)
                                                                        ccTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCaa                    >  2:3375862/1‑93 (MQ=255)
                                                                         cTTCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAac                   >  1:2809120/1‑93 (MQ=255)
                                                                          ttCCGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAaca                  >  2:1659205/1‑93 (MQ=255)
                                                                            ccGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAATGCGTACCTGGCAAcaca                >  1:4216261/1‑93 (MQ=255)
                                                                            ccGCGATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAAcaca                <  2:4045868/93‑1 (MQ=255)
                                                                             cgcgATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAg               >  2:4611867/1‑93 (MQ=255)
                                                                             cgcgATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAg               >  1:4322694/1‑93 (MQ=255)
                                                                              gcgATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAgg              >  2:550763/1‑93 (MQ=255)
                                                                              gcgATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAgg              >  2:3651598/1‑93 (MQ=255)
                                                                              gcgATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAgg              >  1:2143799/1‑93 (MQ=255)
                                                                               cgATGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAggg             <  2:134571/93‑1 (MQ=255)
                                                                                 aTGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACATAGGTCg           >  2:3308700/1‑93 (MQ=255)
                                                                                  tGCCGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGt          <  2:1910391/93‑1 (MQ=255)
                                                                                    ccGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTgg        >  1:2642397/1‑93 (MQ=255)
                                                                                     cGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGga       >  2:3475343/1‑93 (MQ=255)
                                                                                     cGACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGga       <  1:435488/93‑1 (MQ=255)
                                                                                      gACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGgag      >  2:4177253/1‑93 (MQ=255)
                                                                                      gACGCTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGgag      >  1:2511689/1‑93 (MQ=255)
                                                                                          ttGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCa  <  1:3067807/93‑1 (MQ=255)
                                                                                          |                                                                                            
AAGTGCGTGAAGTTTGCCACGTCCCGCTGATTGCCTCCGGTGGCGCGGGCACCATGGAACACTTCCTCGAAGCCTTCCGCGATGCCGACGTTGACGGCGCGCTGGCAGCTTCCGTATTCCACAAACAAATAATCAATATTGGTGAATTAAAAGCGTACCTGGCAACACAGGGCGTGGAGATCA  >  NC_000913/2094449‑2094631

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