Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,316,594 (C)5→4 coding (1219/1362 nt) 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,5900C.98.9% 389.9 / ‑3.1 88coding (1223/1362 nt)trpCN‑(5‑phosphoribosyl)anthranilate isomerase and indole‑3‑glycerolphosphate synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (43/44);  minor base G (0/1);  total (43/45)
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.91e-01

CTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCT  >  NC_000913/1316498‑1316676
                                                                                            |                                                                                      
cTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGc                                                                                        <  2:3613674/93‑1 (MQ=255)
cTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGc                                                                                        <  1:3697727/93‑1 (MQ=255)
 tttGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcc                                                                                       <  1:4683665/93‑1 (MQ=255)
 tttGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcc                                                                                       >  1:4494216/1‑93 (MQ=255)
 tttGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcc                                                                                       >  1:2644803/1‑93 (MQ=255)
  ttGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGccc                                                                                      <  2:4572134/93‑1 (MQ=255)
   tGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcccc                                                                                     <  1:770303/93‑1 (MQ=255)
    gATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcccct                                                                                    >  1:2174184/1‑92 (MQ=255)
     aTGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcccctg                                                                                   >  1:874169/1‑91 (MQ=255)
      tGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcccctgc                                                                                  >  2:2432165/1‑90 (MQ=255)
      tGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGcccctgc                                                                                  <  1:3624408/93‑4 (MQ=255)
           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:3875354/93‑1 (MQ=255)
              tGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGaa                                                                         <  2:2386480/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:18521/1‑93 (MQ=255)
                             gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                          <  1:1383140/93‑1 (MQ=255)
                             gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                          >  1:290970/1‑93 (MQ=255)
                             gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                          >  1:4674296/1‑93 (MQ=255)
                             gAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGAtt                                                          >  2:2005968/1‑93 (MQ=255)
                               aTTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGa                                                        >  1:4654634/1‑93 (MQ=255)
                               aTTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGa                                                        >  2:3025566/1‑93 (MQ=255)
                               aTTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGa                                                        >  2:1260171/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                                                <  2:1791320/93‑1 (MQ=255)
                                       cAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTaa                                                <  1:3936581/93‑1 (MQ=255)
                                       cAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTaa                                                <  2:2939879/93‑1 (MQ=255)
                                       cAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTaa                                                >  1:1029264/1‑93 (MQ=255)
                                        aaGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAAt                                               <  2:1129978/93‑1 (MQ=255)
                                         aGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATa                                              <  2:2017373/93‑1 (MQ=255)
                                         aGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATa                                              <  1:2988856/93‑1 (MQ=255)
                                         aGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATa                                              >  1:2484550/1‑93 (MQ=255)
                                           tCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGt                                            <  1:3537843/93‑1 (MQ=255)
                                             cGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGa                                          >  1:1329017/1‑93 (MQ=255)
                                              ggCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGAc                                         <  2:1232392/93‑1 (MQ=255)
                                              ggCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGAc                                         >  1:3165144/1‑93 (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:3121883/93‑1 (MQ=255)
                                                  cAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGt                                     >  2:3450553/1‑93 (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:1847414/93‑1 (MQ=255)
                                                       ggTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCaaaa                                >  2:4068766/1‑93 (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:3229451/93‑1 (MQ=255)
                                                                            gTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTgg           <  1:2028674/93‑1 (MQ=255)
                                                                            gTTATCTGCGCCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTgg           >  1:2494100/1‑93 (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:4051439/93‑1 (MQ=255)
                                                                                    cgcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTc   <  1:1923343/93‑1 (MQ=255)
                                                                                    cgcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTc   <  2:2005330/93‑1 (MQ=255)
                                                                                    cgcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTc   <  1:954461/93‑1 (MQ=255)
                                                                                     gcCTAAGCCCCTGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCt  >  1:171648/1‑93 (MQ=255)
                                                                                            |                                                                                      
CTTTGATGCCCGGTTGCGACTCTACAGCAGAATTAAAATCAAGTCCGGCGCAGCCGGTTTGTGCCGCTTCCACGCAGTTATCTGCGCCTAAGCCCCCCGCCAGCAGAACGTTGCCAAGCGATTGACCATTTAATAGTGACCAGTCAAAACGTTGCCCGCTTCCACCCTGGCCGTTGTCT  >  NC_000913/1316498‑1316676

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