Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,316,595 C→T A406A (GCG→GCA trpC ← N‑(5‑phosphoribosyl)anthranilate isomerase and indole‑3‑glycerolphosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,316,5950CT100.0% 318.3 / NA 89A406A (GCG→GCAtrpCN‑(5‑phosphoribosyl)anthranilate isomerase and indole‑3‑glycerolphosphate synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (43/46);  total (43/46)

GGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCTAAAACA  >  NC_000913/1316509‑1316682
                                                                                      |                                                                                        
ggTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCcagca                                                                                   >  1:3576395/1‑93 (MQ=255)
  ttGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGa                                                                                 <  2:534531/93‑1 (MQ=255)
  ttGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGa                                                                                 <  1:4117705/93‑1 (MQ=255)
  ttGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCACTGCCAGCAGa                                                                                 >  1:3257951/1‑93 (MQ=255)
   tGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGaa                                                                                <  2:2386480/93‑1 (MQ=255)
   tGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGaa                                                                                <  2:3875354/93‑1 (MQ=255)
    gCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAAc                                                                               >  1:2005330/1‑93 (MQ=255)
      gACTCTGGAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGt                                                                             <  2:20520/93‑1 (MQ=255)
      gACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGt                                                                             <  1:1217292/93‑1 (MQ=255)
      gACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGt                                                                             <  1:3207396/93‑1 (MQ=255)
        ctctACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTg                                                                           >  1:4535501/1‑93 (MQ=255)
         tctACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGc                                                                          >  2:206845/1‑93 (MQ=255)
         tctACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGc                                                                          >  2:3149792/1‑93 (MQ=255)
           tACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCa                                                                        >  2:2641419/1‑93 (MQ=255)
           tACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCa                                                                        >  1:4629656/1‑93 (MQ=255)
             cagcagAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAg                                                                      <  2:2969068/93‑1 (MQ=255)
              agcagAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGc                                                                     <  2:1869798/93‑1 (MQ=255)
               gcagAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCg                                                                    >  1:1650543/1‑93 (MQ=255)
                 agAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAt                                                                  <  1:2586683/93‑1 (MQ=255)
                 agAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAt                                                                  >  2:2663406/1‑93 (MQ=255)
                  gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                                 >  1:4674296/1‑93 (MQ=255)
                  gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                                 >  2:2005968/1‑93 (MQ=255)
                  gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                                 >  1:290970/1‑93 (MQ=255)
                  gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                                 <  1:1383140/93‑1 (MQ=255)
                  gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                                 >  1:18521/1‑93 (MQ=255)
                    aTTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGa                                                               >  2:3025566/1‑93 (MQ=255)
                    aTTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGa                                                               >  2:1260171/1‑93 (MQ=255)
                    aTTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGa                                                               >  1:4654634/1‑93 (MQ=255)
                     ttAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGAc                                                              <  2:234463/93‑1 (MQ=255)
                       aaaaTCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCa                                                            <  2:4254229/93‑1 (MQ=255)
                         aaTCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCAtt                                                          >  1:270334/1‑93 (MQ=255)
                         aaTCAAGTCCGCCGCAGCCGGTTTGTGCCGCTCCCACGCAGTTGTCTGCGCCTAGGCCCCTGCCAGCAGAACGGTGCCAAGCGATTGACCAtt                                                          <  2:4576012/93‑1 (MQ=255)
                           tCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTa                                                        <  2:1334400/93‑1 (MQ=255)
                            cAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTaa                                                       <  1:3936581/93‑1 (MQ=255)
                            cAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTaa                                                       <  2:2939879/93‑1 (MQ=255)
                            cAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTaa                                                       <  2:1791320/93‑1 (MQ=255)
                            cAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTaa                                                       >  1:1029264/1‑93 (MQ=255)
                             aaGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAAt                                                      <  2:1129978/93‑1 (MQ=255)
                              aGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATa                                                     <  1:2988856/93‑1 (MQ=255)
                              aGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATa                                                     >  1:2484550/1‑93 (MQ=255)
                              aGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATa                                                     <  2:2017373/93‑1 (MQ=255)
                                tCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGt                                                   <  1:3537843/93‑1 (MQ=255)
                                  cGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGa                                                 >  1:1329017/1‑93 (MQ=255)
                                   ggCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGAc                                                >  1:3165144/1‑93 (MQ=255)
                                   ggCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGAc                                                <  2:1232392/93‑1 (MQ=255)
                                    gcgcAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGAcc                                               <  1:1899159/93‑1 (MQ=255)
                                    gcgcAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGAcc                                               >  2:4278606/1‑93 (MQ=255)
                                      gcAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAg                                             <  2:2156306/93‑1 (MQ=255)
                                       cAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGt                                            >  2:3450553/1‑93 (MQ=255)
                                       cAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGt                                            <  2:3121883/93‑1 (MQ=255)
                                         gCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCa                                          <  1:4137952/93‑1 (MQ=255)
                                          ccGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCaa                                         >  1:1528918/1‑93 (MQ=255)
                                           cGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCaaa                                        >  1:3420819/1‑93 (MQ=255)
                                            ggTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCaaaa                                       >  2:1632002/1‑93 (MQ=255)
                                            ggTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCaaaa                                       >  2:4068766/1‑93 (MQ=255)
                                            ggTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCaaaa                                       <  2:1847414/93‑1 (MQ=255)
                                             gTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGCTTGACCATTTAATAGTGACCAGTCAAAAc                                      >  2:3008541/1‑93 (MQ=255)
                                             gTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAAc                                      >  1:4334240/1‑93 (MQ=255)
                                             gTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAAc                                      <  1:2198690/93‑1 (MQ=255)
                                             gTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAAc                                      <  2:2762337/93‑1 (MQ=255)
                                              tttGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACg                                     >  2:2384072/1‑93 (MQ=255)
                                               ttGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGCCCAGTCAAAACGt                                    <  1:3398688/93‑1 (MQ=255)
                                                tgtgCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGtt                                   <  1:3422608/93‑1 (MQ=255)
                                                  tgCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGc                                 >  2:490640/1‑93 (MQ=255)
                                                        ttCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCtt                           <  1:653819/93‑1 (MQ=255)
                                                        ttCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCtt                           <  1:1329297/93‑1 (MQ=255)
                                                          ccACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTcc                         >  2:2203923/1‑93 (MQ=255)
                                                             cGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCAcc                     >  1:3854719/1‑93 (MQ=255)
                                                              gCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCAccc                    >  2:1001555/1‑93 (MQ=255)
                                                               cAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCt                   >  2:157370/1‑93 (MQ=255)
                                                               cAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGGGACCAGTCAAAACGTTGCCCGCTTCCACCCt                   >  2:1372631/1‑93 (MQ=255)
                                                                aGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTg                  >  1:1011930/1‑93 (MQ=255)
                                                                aGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTg                  <  1:1601144/93‑1 (MQ=255)
                                                                 gTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTgg                 <  1:2028674/93‑1 (MQ=255)
                                                                 gTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTgg                 >  1:2494100/1‑93 (MQ=255)
                                                                 gTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTgg                 <  1:3229451/93‑1 (MQ=255)
                                                                  ttATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGc                <  1:4363139/93‑1 (MQ=255)
                                                                  ttATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGc                <  2:3932370/93‑1 (MQ=255)
                                                                  ttATCTGCACCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGc                <  1:2766561/93‑1 (MQ=255)
                                                                   tATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGcc               >  1:279856/1‑93 (MQ=255)
                                                                    aTCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCg              >  2:3322081/1‑93 (MQ=255)
                                                                      cTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGtt            <  1:3342728/93‑1 (MQ=255)
                                                                       tGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTg           <  2:3142669/93‑1 (MQ=255)
                                                                        gcgcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGCTGt          >  2:2810678/1‑93 (MQ=255)
                                                                         cgcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTc         <  1:1923343/93‑1 (MQ=255)
                                                                         cgcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTc         <  2:2005330/93‑1 (MQ=255)
                                                                         cgcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTc         <  1:4051439/93‑1 (MQ=255)
                                                                         cgcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTc         <  1:954461/93‑1 (MQ=255)
                                                                          gcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCt        >  1:171648/1‑93 (MQ=255)
                                                                               aagcccctGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCTaaaa    <  2:689279/90‑1 (MQ=255)
                                                                                agcccctGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCTAAAAc   <  1:4467811/91‑1 (MQ=255)
                                                                                 gcccctGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCTAAAACa  >  2:4404871/2‑93 (MQ=255)
                                                                                 acccctGCCAGCAGAACGTTTCCAAGCGATTGACCATTTAATAATTACCAGGCAAAAAGTTGCCCGCTTTCACCCCTGCCCGTTGCTTAAACa  >  2:411150/2‑93 (MQ=255)
                                                                                      |                                                                                        
GGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCTAAAACA  >  NC_000913/1316509‑1316682

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