Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,091,351 C→T G168D (GGC→GAC)  rhaD ← rhamnulose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,091,3510CT99.4% 549.0 / ‑5.7 160G168D (GGC→GAC) rhaDrhamnulose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (67/92);  minor base A (1/0);  total (68/92)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 172 total aligned reads displayed.

GCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091261‑4091440
                                                                                          |                                                                                         
   ggcGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTcc                                                                                      <  2:4064471/91‑1 (MQ=255)
    tCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTccc                                                                                     <  1:2406535/93‑1 (MQ=255)
        cgggCCGATCTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTATTCCCACAg                                                                                 <  1:1014973/91‑1 (MQ=255)
         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCTACAGt                                                                                <  1:3377458/93‑1 (MQ=255)
         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  1:255867/1‑93 (MQ=255)
          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGGCAACGACATCCGGGAATACCACCAGACACTCGGTACTGTCTTCCCACAGtt                                                                               >  2:4392023/1‑93 (MQ=255)
          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               <  2:2328865/93‑1 (MQ=255)
          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               >  1:640392/1‑93 (MQ=255)
          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               >  1:3768377/1‑93 (MQ=255)
          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               >  1:1699199/1‑93 (MQ=255)
           gCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTg                                                                              <  1:2834708/93‑1 (MQ=255)
             cGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGc                                                                            >  1:4274459/1‑93 (MQ=255)
              gaattCGTCCGTGCCGTGCACCTTTCACGGCAAAATGCCTACTCCATCCTGGAATACTACCAGACACTCGGTGCTGTTTTCCCACAGCTGGCg                                                                           >  2:4211342/4‑93 (MQ=255)
              gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           >  1:2244196/1‑93 (MQ=255)
              gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           >  1:1482311/1‑93 (MQ=255)
              gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           >  1:2803320/1‑93 (MQ=255)
               aTTTCGTCCGTGCCGGGCCCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  1:2646583/93‑1 (MQ=255)
               aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  2:4746502/93‑1 (MQ=255)
               aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  2:2626884/93‑1 (MQ=255)
                tttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAg                                                                         >  2:605647/1‑93 (MQ=255)
                   cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                                      <  1:724766/93‑1 (MQ=255)
                    gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                                     <  2:3965258/93‑1 (MQ=255)
                    gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                                     <  2:1698600/93‑1 (MQ=255)
                     tCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                    >  2:4550083/1‑93 (MQ=255)
                      ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                   <  1:4001108/93‑1 (MQ=255)
                       cGTGCCGGGCACCATCCACGTCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  2:2874100/93‑1 (MQ=255)
                       cGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  2:3283952/93‑1 (MQ=255)
                       cGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  2:4121508/93‑1 (MQ=255)
                       cGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  1:240111/93‑1 (MQ=255)
                       cGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  2:2025932/93‑1 (MQ=255)
                       cGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  2:2116479/93‑1 (MQ=255)
                       cGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                  <  1:4570110/93‑1 (MQ=255)
                        gTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAcc                                                                 >  2:2919724/1‑93 (MQ=255)
                           ccGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgcg                                                              <  2:2800118/93‑1 (MQ=255)
                           ccGGGCAACATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgcg                                                              <  1:1555055/93‑1 (MQ=255)
                               gCACCATCCACGGCAAAATGCCACCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  2:2393711/93‑1 (MQ=255)
                               gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          >  1:2255951/1‑93 (MQ=255)
                               gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  1:1688297/93‑1 (MQ=255)
                                cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTc                                                         <  1:3524332/93‑1 (MQ=255)
                                 aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                                        >  2:1560280/1‑93 (MQ=255)
                                   cATCCACGGCAAAATGCCAACGCCATCCGGGAATACAACCAGACACTCGGTGCGTTCTTCCCACAGTTGCCGAGTGAAGACCGCGGGGTCGtt                                                      >  2:48973/1‑93 (MQ=255)
                                    aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTCGCGACTGAAGACCGCGGTGTCGttt                                                     >  1:1893877/1‑93 (MQ=255)
                                    aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCCTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                     >  2:3723502/1‑93 (MQ=255)
                                     tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                                    <  1:1289308/93‑1 (MQ=255)
                                      ccACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTc                                                   <  2:4132098/93‑1 (MQ=255)
                                        ccGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCaa                                                 <  1:3558487/92‑1 (MQ=255)
                                         cGGCAAAATGCCAACGCCATCCGGGAATGCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAg                                                <  1:1229223/93‑1 (MQ=255)
                                         cGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAg                                                <  2:1156654/93‑1 (MQ=255)
                                          ggCAAAATGCCAACGCCATCCGGGAATACCGCCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                               <  1:1982961/93‑1 (MQ=255)
                                           gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTa                                              >  1:2243405/1‑93 (MQ=255)
                                           gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTa                                              >  1:3144074/1‑93 (MQ=255)
                                            cAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAc                                             >  1:4403843/1‑93 (MQ=255)
                                             aaaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                            <  2:1124156/93‑1 (MQ=255)
                                             aaaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                            <  2:295307/93‑1 (MQ=255)
                                              aaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACAt                                           <  1:3309433/93‑1 (MQ=255)
                                               aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          >  1:4529652/1‑93 (MQ=255)
                                               aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          <  2:2773835/93‑1 (MQ=255)
                                                aTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                         <  2:33483/93‑1 (MQ=255)
                                                 tGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAgg                                        <  2:4588208/93‑1 (MQ=255)
                                                  gCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGt                                       >  1:1596460/1‑93 (MQ=255)
                                                  gCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGCGAAGACCGCGGTGTCGTTTTCAAGTACATAGGt                                       <  1:4189221/93‑1 (MQ=255)
                                                     aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    <  1:267541/93‑1 (MQ=255)
                                                     aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    >  1:32397/1‑93 (MQ=255)
                                                     aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGACTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    >  1:4125547/1‑93 (MQ=255)
                                                      aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                                   >  2:2586797/1‑93 (MQ=255)
                                                      aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                                   <  2:3596005/93‑1 (MQ=255)
                                                           aTCCGGGAATCCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTGCATAGGTGAGGGCGAt                              <  1:1809616/93‑1 (MQ=255)
                                                           aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGAt                              >  2:3447241/1‑93 (MQ=255)
                                                            tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             >  2:1792998/1‑93 (MQ=255)
                                                             ccGGGAATCCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            <  1:3809659/93‑1 (MQ=255)
                                                              cGGGAATCCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  2:1898937/93‑1 (MQ=255)
                                                              cGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  2:345826/93‑1 (MQ=255)
                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATATGTGAGGGCGATCAGGt                         <  1:4507386/93‑1 (MQ=255)
                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  2:3762605/93‑1 (MQ=255)
                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  1:2741979/93‑1 (MQ=255)
                                                                 gAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGtt                        <  1:2132057/93‑1 (MQ=255)
                                                                    tcccaccaGACACTCGGTGCTGTCTTCTCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt                     <  1:3460151/91‑1 (MQ=255)
                                                                    tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt                     <  1:4242219/93‑1 (MQ=255)
                                                                     accaccAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg                    <  1:3830409/93‑1 (MQ=255)
                                                                         ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                >  2:3836429/1‑93 (MQ=255)
                                                                         ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                <  2:720223/93‑1 (MQ=255)
                                                                         ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGATGgtggcgg                >  2:2208194/1‑92 (MQ=255)
                                                                          caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtg               <  1:102830/93‑1 (MQ=255)
                                                                          caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtg               <  1:4744442/93‑1 (MQ=255)
                                                                          caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGGGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtg               >  1:796678/1‑93 (MQ=255)
                                                                           aGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtgg              >  1:4017550/1‑93 (MQ=255)
                                                                            gACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtggc             >  2:238544/1‑93 (MQ=255)
                                                                             acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCa            <  1:194341/93‑1 (MQ=255)
                                                                               acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt          >  2:185593/1‑93 (MQ=255)
                                                                                cTCGGTGCTGTCTTCCCACAGTTGGCGCGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGtc         >  2:4461524/1‑92 (MQ=255)
                                                                                  cGGTGCTGTCTTCCCGCAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTAGGTGGCGTGGCAGTGCa       <  2:4224118/93‑1 (MQ=255)
                                                                                  cGGTGCTGTCTTCCCACAGTTTGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa       >  2:311669/1‑93 (MQ=255)
                                                                                   ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCAt      <  2:2440824/93‑1 (MQ=255)
                                                                                    gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATg     >  2:1292729/1‑93 (MQ=255)
                                                                                     tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGCTCAGGTTGGTGGCGTGGCCGTGCCTGa    >  2:1304331/1‑93 (MQ=255)
                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   >  2:3701576/1‑93 (MQ=255)
                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   <  2:4484753/93‑1 (MQ=255)
                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   <  1:643295/93‑1 (MQ=255)
                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   <  2:3145942/93‑1 (MQ=255)
                                                                                       ctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc  <  2:1319156/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091261‑4091440

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